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systems with a short-scale repulsive pair interaction converges to the solution of the cubic non- 
linear Schrodinger equation in a suitable scaling limit. The result is extended to fc-particle density 
' matrices for all positive integer k. 
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1 Introduction 



The fundamental principle of quantum mechanics states that a quantum system of N particles is 
described by a wave function of N variables satisfying a Schrodinger equation. In realistic systems, 
N is so large that a direct solution of the Schrodinger equation for interacting systems is clearly 
an impossible task. Many-body dynamics is thus traditionally approximated by simpler effective 
dynamics where only the time evolution of a few cumulative degrees of freedom is monitored. In the 
simplest case only the one-particle marginal densities are considered. The many-body pair interaction 
is then replaced by an effective non-linear mean-field potential and higher order quantum correlations 
are neglected. 

For Bose systems, there are three important examples: i) the Hartree equation describing Bose 
systems with soft or Coulomb interactions in the mean field limit; ii) the cubic non- linear Schrodinger 
(NLS) equation for Bose systems with short range interactions in suitable scaling limits and iii) the 
Gross-Pitaevskii equation for the condensate of Bose systems with short range interactions (including 
the hard core interaction) in a certain scaling limit. All three equations are non-linear Schrodinger 
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equations and the Gross-Pitaevskii equation is itself a cubic non-linear Schrodinger (NLS) equation. 
The fundamental reason why nonlinear Schrodinger equation describes these three systems is the 
following observation: the effective dynamics of a given particle is governed by the density times 
the two-body scattering process. In particular, all three-body scattering processes are negligible in 
the limit; a fact that should be established rigorously. We now introduce a model covering all these 
cases. 

Before giving the precise definitions, we explain the physical differences among these three effec- 
tive theories. Let Yli<i<j<N U(xi ~~ x j) be the many-body pair interaction, where Xj is the position 
of the j-th particle. The largest lengthscale of the problem is the size of the system and it is typically 
comparable with the variation scale of the particle density g{x). We will set this scale to be 0(1). 
Thus the density of the system, g, is of order TV. The potential U defines two lengthscales: the range 
of U and the scattering length of U (see Appendix for a definition), denoted by tjj and ajj, respec- 
tively. The scattering length determines the effective lengthscale of the two particle correlations. We 
shall set the scattering length ajj ~ 0(N~ 1 ) so that gajj is order one. This is needed to obtain an 
effective one-particle dynamics in our model. 

The case i), the Hartree equation, is obtained when rjj ~ 0(1), i.e. the range of U is comparable 
with spatial variation of the density. In this case, the effective two-body scattering process is of 
mean-field type, and each particle is subject to an effective potential U * g. If au <C r\j <C 1, i.e. the 
range of U is much shorter than the spatial variation of the density but bigger than the scattering 
length, the effective two-body scattering process is its Born approximation and the effective potential 
is g f U . Mathematically, this can be explained by the fact that the weak limit of U(x) in this scaling 
is 8{x) J U. Since g(x) is quadratic in the quantum mechanical wave function, one obtains a cubic 
nonlinear Schrodinger equation for the one-particle orbitals. Finally, the Gross-Pitaevskii equation 
arises when the potential U is so localized that its range is comparable with its scattering length, 
au ~ r u (i n particular, this is the case of the hard-core interaction). In this case, the effective 
two-body scattering is the full two-body scattering process and the effective potential is Sirgajj. 

We now give the precise definition of this model. We consider a system of N interacting bosons 
in d = 3 dimensions. The state space of the iV-boson system is L 2 (M. 3N , dx), the subspace of 
L 2 (M. 3N , dx) containing all functions symmetric with respect to permutations of the N particles. Let 
V be a smooth, compactly supported, non-negative and symmetric ( i.e. V(x) = V(—x)) function. 
Define the rescaling of V by 

V N (x) := N 3f3 V(NPx) , (1.1) 

where [3 is a nonnegative parameter. The Hamiltonian with pair interaction ^Vjv(xj — Xj) is given 
by the nonnegative self-adjoint operator 

N 1 

j=l i<j 

acting on L 2 S (R 3N , dx) . Denote by tpN,t the wave function at time t; it satisfies the Schrodinger 
equation 

idtipN,t = H N ip N , t , (1.3) 

with initial condition Vjv.o- The Schrodinger equation conserves the energy, the L 2 -norm and the 
permutation symmetry of the wave function. 

Instead of describing the system through its wave function ipN,ti we can introduce the correspond- 
ing density matrix jN,t, defined as the orthogonal projection onto ipN,t m the space L 2 (M 3Ar ,dx), 
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i.e., 7iv,t = ^jvt- The Schrodinger equation (jl.3p assumes then the Heisenberg form 

id tlN)t = [H N ,y Ntt ] , [A, B] := AB-BA. (1.4) 
Since || Viv,t || = 1, it follows that Tr77v,t = 1. 

(k) 

For any integer k = 1, . . . , N, we define the k-particle marginal density, j N t , by taking the partial 
trace of jN,t over the last N — k particles. If 7jv,t (x; x') := ^jv,t (x)f/>Ar,t(x') denotes the kernel of 7jv,i> 
then the kernel of 7jy t is given by 

7wl(xfc;x^.)=/ dx A r_ fc 7 A r it (x fc ,x 7V -fc;Xfc,x A r_ A .) . (1.5) 

Here and henceforth we use the notation x = (xi, . . . , £Ejv), x& = (xi, . . . , xjt), Xjv-fc = . . . , £Cjv)j 

and similarly for the primed variables. Due to the permutational symmetry of ipNt m an variables, 
the marginal densities are also symmetric in the sense that 

(fc) / / / I \ (k) I i i i \ it ct\ 

lN,t\ Xl i X2 i • • • ' X li X 2' ■ ■ ■ ) X k) — 7jV,t ^ X cr(l) j x a{2) ) • • • j a; <r(&) i ' X <r(2) ' • • • > X o-(fc) J U- D J 

for any permutation a £ Sk (<?& denotes the set of permutations on k elements). Denote by CT the 
unitary operator 

(Q a iJ))(xx, . . . ,x k ) = ipix^i),. • -,x a(k) ) . (1.7) 

Then (jl.6p is equivalent to 

e CT7 (fe) e CT -i = 7 (fc) (i.8) 

for all a G By definition, density matrices are non- negative trace class operators, 7W > 0, acting 
on L 2 (R 3fc ) and with permutational symmetry (jl.6p . With a slight abuse of notation we will use the 
same notation for the operators and their kernels. 

The two-body potential of this model is U = N^Vn- By scaling, the scattering length a\j = 
0(N~ 1 ). The range of interaction, ru, is of order r\j = 0(N~@). Cases i), ii) and iii) correspond to 
the cases = 0, < (3 < 1 and = 1, respectively. The focus of this paper is to study case ii). Our 
main result, the following Theorem ll.l| states that the time evolution of the one-particle density 
matrix is given by a cubic non-linear Schrodinger equation, provided < < 1/2. The same result 
is expected to hold for all < < 1; the regime > 1/2 is an open problem. The topology and 
spaces used in this theorem, e.g., the weak* topology of C\, will be defined in next section. We state 
the theorem only for the dimension d = 3, but the proof can be extended to d < 2 as well. 

Theorem 1.1. Fix if £ // 1 (M 3 ), with \\tp\\ = 1. Assume the Hamiltonian Hn is defined as in il.2\) 
with < < 1/2. Suppose that the unsealed potential V is smooth, compactly supported, positive 
and symmetric, i.e. it satisfies V(x) = V(—x). Let V>iv(x) := Yij=i l P( x j)! and suppose ipN,t is the 

solution of the Schrodinger equation il.3\) with initial data tpN,t=o = ipN- Let Fjyj = {7yvj}/b=i ^ e 
the family of marginal distributions associated to tpN,t- Then, for every fixed t £ R and integer k > 1 
we have 

7$ -7?° asN^oo (1.9) 
with respect to the weak* topology of C\. Here 



x;xO:=n<M^)W(4), (1-10) 



3 



where (ft is the solution of the nonlinear Schrddinger equation 

idttpt = -Aip t + b \ipt\ 2 pt (1-11) 
with initial data <ft=o = f and bo = j V(x)dx. 

For dimension d = 1, Adami, Golse and Teta recently obtained a result similar to Theorem 11.11 
in [2] (certain partial results, for the one-dimensional case, were already obtained by these authors, 
together with Bardos, in [J). 

Theorem 11.11 holds also for (3 = 0. The nonlinear equation (jl.lip then becomes the Hartree 
equation 

id t u t = -Au t + (V *\u t \ 2 )u t . (1.12) 

In this special case, the result was rigorously proven by Hepp [15] for smooth potentials V{x) and 
by Spohn [22j for bounded potentials. Ginibre and Velo |14j treated integrable potentials, but they 
required the initial state to be coherent. In particular, in the approach of [13], the number of particles 
cannot be fixed. For Coulomb potential, partial results were obtained in [3] and a complete proof 
was given in [11]. If the particles have a relativistic dispersion then (|1.12p has to be replaced by 

id t u t = (1 - A) 1 / 2 ^ + (V * \u t \ 2 )u t . 

This equation was recently derived in [8], starting from many-body dynamics, for the case of a 
Coulomb potential. A concise overview on results and open problems related to the Hartree equation 
in physical context is found in [12] . 

The basic strategy to prove Theorem II .11 is the same as in [11]: We first write down the BBGKY 
hierarchy (|2.ip for the marginal densities (jl.5p . Then we take the limit 7V — > oo to obtain an 
infinite hierarchy of equations (|2.2p . Since this hierarchy was first mentioned in the context of the 
Gross-Pitaevskii scaling, (3 = 1, we will continue to call it Gross-Pitaevskii hierarchy even when 
the coefficient of the nonlinear term is given by the Born approximation of the scattering length. 
Finally we prove that the Gross-Pitaevskii hierarchy has a unique solution in a suitable space. Since 
tensor products of solutions to the non-linear Schrodinger equation (jl.lip are trivial solutions to the 
hierarchy, this identifies the limit and thus proves Theorem ll.il 

The key steps in this approach are an a-priori estimate, the convergence of the BBGKY hierarchy 
to the infinite hierarchy and the proof of the uniqueness of the infinite hierarchy. The first part was 
already proved in [TJ. The convergence has to be proven in a somewhat stronger sense than in [7J. 
The key point of the present paper is the uniqueness result stated as Theorem 19.11 This theorem 
is indeed the well-posedness of the Gross-Pitaevskii hierarchy. Since solutions of the cubic non- 
linear Schrodinger equation naturally generate solutions of the hierarchy by taking tensor products, 
Theorem 19.11 can be viewed as an extension of the well-posedness theorem of the NLS equation to 
infinite dimensions. Therefore, we have to either extend the fundamental tool in the well-posedness of 
the nonlinear Schrodinger equation, the Stricharz inequality, to infinite dimensions or find a method 
avoiding it (for a review on the Stricharz inequality and the well-posedness of nonlinear Schrodinger 
equations see, e.g., [1], [5], [21]). Apart from this issue, we have to control correlation effects of 
many-particle systems which are absent in the nonlinear equation in Euclidean space. Our method, 
based on the analysis of Feynman graphs, provides a solution to both problems and contains in 
certain sense a version of Stricharz inequality in infinite dimension — albeit we deal only with cubic 
nonlinearity. We will explain this issue in the remarks after Theorem 19.31 
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We emphasize that our uniqueness result is valid for any coupling constant in the non-linear 
Schrodinger equation, in particular it is valid also for the special case of the Gross-Pitaevskii equation. 
Thus the main part of the paper actually is independent of the scaling in the iV-body model, in 
particular it is independent of the choice of 0. The restriction < 1/2 is only used to obtain the 
a-priori bound and therefore the convergence to the infinite hierarchy. 

As mentioned earlier, the coupling constant bo = J V in the non-linear Schrodinger equation in 
Theorem 11.11 is the Born approximation to the physically correct coupling constant: 8ir times the 
scattering length a at of the potential j^Vn- To see this, denote the scattering length of V by ao- 
Then we have 

Urn Na N = H 8 * if ° < " < (1.13) 

iV^oo | a if = 1. 

The limit in the first case is proved in Lemma IA.1I for radial potential; the second one is just a 
rescaling. Notice that at = 1, the coupling constant is the scattering length of the unsealed 
potential, which indicates that the full two-body scattering process needs to be taken into account. 
For stationary problem in this case, the connection between the ground state energy of the Bose 
system and the Gross-Pitaevskii energy functional was rigorously established by Lieb and Yngvason 
[18j . Furthermore, the existence of Bose-Einstein condensation in this limit was proved by Lieb and 
Seiringer [16J. An excellent overview on the recent development in these problems, see [17]. For the 
dynamical problem, we have proved in [9] that the family of the reduced density matrices converges 
to a solution of the Gross-Pitaevskii hierarchy with the correct coupling constant 87rao. However, 
the a priori estimate obtained in [9] was not strong enough to apply the uniqueness theorem in this 
paper, Theorem 19.11 Furthermore, the pair interactions in [9] were cutoff whenever many particles 
are in a very small physically unlikely region. To complete the project for = 1, we still have to 
remove this cutoff and establish the a priori estimate needed for Theorem 19.11 

Acknowledgement. We would like to thank the referees for their helpful comments on how to 
improve the presentation of the results in the manuscript. 

2 The BBGKY Hierarchy 

The marginal density matrices (jl.4j) satisfy the BBGKY hierarchy: 



i=i 

+ N 2 ( V ^( X i- X j)- V N(x'i-x' j ))^l(x k -x' k ) (2.1) 

l<i<j<k 

1 ha 

1 ~ ~m) ^2 / dx ^+i { v N(xj ~ Xfc+l) ~ v N{x'j - Xfc+i)) 7^ 1) (x fc ,x fc+ i;xfc,x fe+ i) 



for k = 1,2, . . . , N. The first term on the right hand side of the hierarchy describes the kinetic 
energy of the first k particles; the second term is associated with the interactions among the first k 
particles, and the last term corresponds to interactions between the first k particles and the other 
N — k particles. 
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Recall that bo = f dx V(x) is the L 1 norm of the non-negative potential V. Since Vn(x) — > 6o<H x ) 
as iV — > oo, the BBGKY hierarchy (|2.ip converges formally, as iV ^ oo, to the following Gross- 
Pitaevskii hierarchy of equations: 



k 



(2.2) 



+ &oE / dx fc+i (<K x i ~~ ~~ ^04 ~ x fe+i)) 7t +1) (xfc,x fe+ i;x' fe ,x fc+ i) 

j'=l 



for any fc > 1. It is easy to see that the family of factorized densities 7 t = IJjLi < Pt(%j) ( Pt(%'j) 1S a 
solution of (|2.2p if and only if ^ is a solution of the cubic nonlinear Schrodinger equation 



id t <pt = -Ay? t + bolvtlVt • ( 2 - 3 ) 
If we can establish the uniqueness of the Gross-Pitaevskii hierarchy, then for every fixed k > 1, 

k 

7$^7t (fc) = n^fo)^). asiV^oo (2.4) 

with respect to some suitable topology. 

As a technical point, we remark that the action of the delta-function on general density matrices 
in (|2.2p is not well defined. However, in our case, the density matrices are in TLk- For such density 
matrices (12. 2h can be defined through an appropriate limiting procedure, see Section [H 

3 Banach Spaces of Density Matrices 

In this section we define some Banach spaces which will be useful in order to take the limit N —* oo 

(k) i 

of the marginal densities T N t - For k > 1, we denote by L k and by K\- the space of trace class and, 
respectively, of compact operators on the fc-particle Hilbert space L 2 (M 3fc ,dx fc ). We have 

(4,||.|| 1 ) = (/C fe ,||.||r , (3.1) 

where || . ||i is the trace norm, and || . || is the operator norm (see, for example, Theorem VI. 26 in 
|19|). The density matrices are the nonnegative elements of C\ with permutational symmetry (|1.8p . 

For 7 ( fc ) G Ch we define the norm 

\h {k) \\ Hk = Tr\S 1 ...S k ^ k) S k ...S 1 \ 
and the corresponding Banach space 

H k = {i [k) e£Ml7 (fc) lk<oo}. 

Moreover we define the space of operators 

A k = {T« = S l . . . S k K^S k . . . Si : G IC k } 

with the norm 



\\T {k) \U k = HSf 1 • • • S k 'T^S^ 1 . . . sr 1 ! 



6 



where | . | is the operator norm. Then, analogously to (|3.ip . we have 

(H k ,\\.\\ Hk ) = (A k ,\\.\\ Ak y (3.2) 

for every k > 1. This induces a weak* topology on 7i k (for a proof of (|3.2p see Lemma 3.1 in 
Since ^4^ is separable, we can fix a dense countable subset of the unit ball of A k - we denote it by 
{j\ }i>i G ^4fc, with ||J^||.4 fc < 1 for all i > 1. Using the operators we define the following 
metric on 7^: for jW^w g ?^ fc we set 

oo 

M 7 (fc W &) ):=E 2-i M fc) (7 (fc) -7 W )| • (3.3) 

Then the topology induced by the metric Pfc(-, •) and the weak* topology are equivalent on the unit 
ball of TCk (see [20] . Theorem 3.16) and hence on any ball of finite radius as well. In other words, a 
uniformly bounded sequence 7^ G TL k converges to 7^ G Hk with respect to the weak* topology, 
if and only if puilN j 7^) ^ as iV 



00. 



For a fixed T > 0, let C([0, T],H. k ) be the space of functions of t G [0, T] with values in which 
are continuous with respect to the metric p k . On C([0,T],Tl k ) we define the metric 

Ml {k) (-)n {k) (■))■■= sup p k (^ k \t),jM(t)). (3.4) 
te[o,T] 

Finally, we define the space 7i as the direct sum over k > 1 of the spaces 7^, that is 

H = ®H k = {r = { 7 (fc) }fe>i : 7 (fe) G H k , Vk > l} , 
fe>i 

and, for a fixed T > 0, we consider the space 

C([O,T],«) = 0C([O,T],W t ), 
fe>i 

equipped with the product of the topologies induced by the metric p^ k ^ on C([0, T],H k )- Let r denote 
this topology. That is, for T N>t = {7jv,t}fc>l and T t = { r -ft^}k>i in C([0, T],H), we have - 1 * T t 
for iV — > 00 if and only if, for every fixed k > 1, 

- / (fe) (fch / (fe) (fch n 

PkUN'vTt ) = sup p k {T N ' v Tt ) 
t€[0,T] 

as — > 00. 

Notation. As in (|1.5|) . we use x = (a?i, . . . , xjy) G R 3iV , x^ = (x\, . . . ,x k ) G R 3fc , x^-A: = 
(x/c+i, . . . , xn) G R 3 ^ - ^, and analogously for the primed variables. We also use the notation 
(x) = (1 + x 2 ) 1 / 2 , for all x G R d . We also set S, = {pj) = (1 - A.,-) 1 / 2 , for all integer j > 1 
(p,- = —iV x . is the momentum of the j-th particle). Moreover Tr^ will indicate for the partial trace 
over the Xj-th variable. The norm notation without subscript, || • ||, will always refer to the L 2 -norm 
for functions and to the operator norm in case of operators. Unless stated otherwise, all integrals are 
over R 3 , or on R 3fc if the measure is dx^,dpfc etc. Universal constants will be denoted by (const.). 
Constants, that may depend on other parameters, will be denoted by C. The dependence is indicated 
in the statements but not always in the proofs. Typically C depends on the initial function ip and 
on the potential V . To compare positive numbers A, B, we also use A < B to indicate that there is 
a universal constant (const.) > with A < (const. )B. The fact that A < B and B < A is denoted 
by A ~ B. 
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4 Outline of the Proof of the Main Theorem 

The proof of Theorem 11,11 is divided into several steps. In the following we fix T > arbitrary. 

Step 1: Regularization of the initial data. Fix k > and x £ Co°(^)! with < x ^ 1) x( s ) = 1) 
for < s < 1, and x( s ) = if s > 2. Define the regularized initial function 

_ X (kH n /N)iP n 



\ X (kH n /N)iI> n \\' 

and we denote by ip^t ^ ne solution of the Schrodinger equation (|1.3|) with initial data tp 1 ^. Denote 

by Tnj = {1n\}^=i family of marginal densities associated with ij)% t . The tilde in the notation 
indicates dependence on the cutoff parameter k. 

This regularization cuts off the high energy part of ipN and it allows us to obtain the strong 
a-priori estimate 

Tr (1 - Ax) . . . (1 - A fe )7}g < C fc (4.1) 

for sufficiently large N and uniformly in t (see Theorem 17. ip . Here the constant C depends on the 
cutoff k > 0. We thus have the compactness stated in the following step. 

Step 2: Compactness o/r^t- For fixed k > 1, it follows from Theorem 17.11 that the sequence 
7^ \ £ C([0, T], TLk) is compact with respect to the topology induced by the metric p k . By a standard 
Cantor diagonalization argument, this implies that Tnj is a compact sequence in C([0, T],7i.) with 

respect to the r-topology. It also follows from Theorem l7.ll that, if To^t = {j^\}k>i £ C{[0,T],H) 

~ ~(k) 

denotes an arbitrary limit point of Tnj, then 7^ is non-negative, symmetric w.r.t. permutations 
in the sense of (II. 8p . and satisfies 

Tr \S 1 ...S k ^] t S k ...S 1 \ <C k (4.2) 

for all k > 1 and t £ [0, T]. We put a tilde in the notation for and C, because a-priori they 
may depend on the cutoff k, which is kept fixed here. Notice that C is independent of k. From Step 
3 and Step 4 below it will follow that is actually independent of k. 

Step 3: Convergence to solutions of the Gross-Pitaevskii hierarchy. Define the free evolution 
operator 

Z4 fc) (i) 7 (fc) = exp LitJ^-^ 7 (fc) exp (^itJ^-A^ . (4.3) 

Theorem 18.11 states that an arbitrary limit point G C([0, T\,H) of the sequence T^t satisfies 
the infinite Gross-Pitaevskii hierarchy in the integral form 



~(k) 

7oo,t 



W W (*)7S!o-^oE / C d S W (fc) (t- S )Tr fe+1 [S( XJ - x k+1 ),^] , (4.4) 



with initial data 

k 



7oo,t=o( x fc; 4) = 7 ( x fc; 4) = n vO^M^-)- ( 4 - 5 ) 
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Notice that (|4.5p implies that the initial data F 



oo,t=0 



is independent of the cutoff k. 



Step 4-' Uniqueness of the solutions of the Gross-Pitaevskii hierarchy. In Theorem 19.11 we prove 
that there is at most one solution T t = {7t }fc>i of (|4.4|) in the space C([0, T],H), such that, for 
all k > 1 and t G [0, T], 7^ is non-negative, symmetric w.r.t. permutation (in the sense (jl.8p ) and 
satisfies (|4.2j) and (|4,5p . Hence the family of factorized densities T t = {7t }fc>i defined in (jl.lOp is 
the unique nonnegative symmetric solution of the Gross-Pitaevskii hierarchy (|4.4p . Thus T/vt — ► F t 

as iV — » 00 in the r-topology and this holds for any fixed k, so the limit is independent of n. Since 

~(k) -t 
is bounded in the TCk norm, uniformly in N, the convergence in the metric Qk and the weak* 

convergence of Tit are equivalent. It therefore follows that for every fixed k > 1, t G [0, T] and k > 0, 
we have 7^! —* 7^ as N —* 00 with respect to the weak* topology of W^- Convergence in weak* 
then trivially follows. 



Step 5: Removal of the cutoff and the conclusion of the proof. It follows from Proposition 15. 1| 
part ii), that 

\\i>N,t-rNA\ = \\i>N-m<CK i /\ 

where the constant C is independent of N and k. This implies that, for every G JCk, we have 



Tr 



(fc) 
TNt 



~(fe) 



(4.6) 



where the constant C depends on j( k \ but is independent of N, k or k. 

For fixed A; > 1 and t, we choose j( fc ) G K.^ and e > 0. Then for any k > 0, we have from (|4,6 
that 



< 



Tr J<*> (7$ 



~{k) 
TNt 



+ 



Tr J« 



(k) 



< C^ 1 / 2 + Tr J« 



7jv,t - It 



(k) 



(4.7) 



Since J*t\ — > 7^ with respect to the weak* topology of Ci, the last term vanishes in the limit 
N —* 00. Since /c > is arbitrary, the r.h.s. of (|4.7p is smaller as £ for N large enough. This 
completes the proof of the theorem. 

Remark. Note that the main body of the proof, Steps 1-4, would have proven weak* convergence 
in H^. It is in the removal of the cutoff, Step 5, that we can prove only under a weaker convergence. 
This situation is similar to the Coulomb interaction paper, [TT], where the removal of the cutoffs 
resulted eventually in a weaker sense of convergence. 



5 Cutoff of the initial wave function 

In this section we show how to regularize the initial wave function ^^(x) = J~J J=1 <p(xj). The aim is 
to find an approximate wave function ip^, depending on a cutoff parameter n > 0, so that, on the 
one hand, the expectation of H% in the state ipfj is of the order N k , and, on the other hand, the 
difference between ip^ and i/jn converges to zero, as k — * 0, uniformly in N. 

Proposition 5.1. Let ip G i? 1 (M 3 ), with ||y||x,2 = 1. Assume is defined as in suppose 
that the unsealed potential V is smooth and positive, and that < j3 < 2/3. We define ?/>Ar(x) = 
UjLi and , for k>0, 

X^iV^v (51) 
9N \\ X (kH n /N)iP n \\ { ' 
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Here x £ Co°0^) *s a cutoff function such that < x 5= lj x( s ) = 1 / or < s < 1 and x( s ) = / or 
s > 2. We denote by *y N , for k = 1, . . . , N , the marginal densities associated with ip^. 

i) For every integer k > 1 and for k > small enough, we have 

W.^)^!^ (5-2) 
where the constant C only depends on the H 1 -norm of if and on the unsealed potential V(x). 

ii) We have 

uniformly in N (C only depends on the H 1 norm of <p and on the unsealed potential V). 
Hi) Let 

k 

ji k \xk]x k ) = ~[l l P(xj) 7 P( x 'j)- (5-3) 
3=1 

Then, for every fixed k > 1 and G JC^, we have 

TrjW - 7 f)^0 (5.4) 

as N — > oo (the convergence is uniform in k > 0, /or k small enough). 
Proof. First we compute 

\\ X (kH n /N)^ n - ip N \\ 2 = (^ N , (1 - x{kH n /N)) 2 ^n} < (Vw, 1(kH n > N)^n) , (5.5) 

where l(s > A) is the characteristic function of [A, oo). Next we use that x( s > 1) < s ) f° r an s > 0. 
Therefore 

\\x(kH n /N)iP n -^ n \\ 2 < ^(ip N ,H N ijj N ) 

= ± (n\\ Vpf + ^Hn 3 ^ n , V(Nf s (x 1 - x 2 ))ip N )\ (5-6) 
< (const. )k {\\<p\\ 2 H i + \\V\\ L i\\<p\\ 2 m ) ' 

where we used that N 3 l 3 V(N /3 (x 1 - x 2 )) < (const.) \\V\\ L i (1 - Ai)(l - A 2 ) (see Lemma El from 
Appendix [A]). Hence 

\\ X (kH n /N)^ n -tP n \\< Ck 1 ! 2 (5.7) 

for a constant C only depending on the H 1 norm of ip and on the unsealed potential V. Using (|5.7|) 
we obtain 

m , H » N , = »»;^f^ H ^ < „„ ' < 2 " Nk 1 



|x(k#w/A0<MI 2 ~ K k 1-2\\ X (kH n /N)^ n -Viv|| ~ K fc 1-Ck 1 / 2 

for all < k < 1/C 2 , which proves i). 

Part ii) follows very easily from (|5.7p because, using the shorthand notation % = x( k Hn /N) and 
using that \\ipn\\ = 1 ; 

^jV ~ i, X ,^1, < II^AT - X^jvll + xVw ~ U X ^ N u = \\i>N ~ X^n\\ + |1 - Hx^MII . . 

Ilxwl IIXWll (5.8) 



< 2||^tv - X^ 



AM 
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Finally, we prove iii). For fixed k > 1, J {k) G JC k and e > we prove that 



TrJ^ (7?-^) 



< e 



(5.£ 



if N is large enough (uniformly in k, for ft sufficiently small). To this end, we choose (p* £ H 2 (M, 3 ) with 
\\<p*\\ = 1, such that \\<p— <p*\\ < e/ (24/c|| ||). Then we define f/!jv,*( x ) = ll?=i V 9 *( x j) rL=fc+l vC^i)) 
and V>tv,* = x( k Hn /N)i()n^/\\x(kHn /N)i(;n^\\. Moreover we set 



~(fe) / / 

7v *V x fc, x fc 



dxTV_ fc ^TV,* ( x fc , ^N-k)^N,* (x' fc , XjV-fc). 



(5.10) 



Note that even though ^>jy * is n °t symmetric with respect to permutation of the N particles, it is 

~(k) 

still symmetric in the first k and the last N — k variables; hence 7)^ is a symmetric density matrix. 
Next we define the Hamiltonian 



H 



N 



A i + jf Yl V N {xi- Xj ), 



j>k+l 



k+l<i<j<N 



with Vn(x) = N 3 Pv(N@x). We denote x = x(^-ffjv /N) and we set Vtv = X^n/Wx^nW and ipN,-- 
)&N,*/\\xi>N,*\\- We also define 



7^(xjfe;x4) := / dxjv.fc ^j\r(Xjfc, X A r_. fc )^Ar(x / fc , Xat_ 



7vl(xfc;x' fc 



(5.11) 



dxAr_ fe VAf,*(xfc, X A r_ A .)^ A r i *(x / fc , Xiv_ fc ) . 



Although and V'iv,* are not symmetric with respect to permutations of the N particles, they 

^(k) 

are still symmetric w.r.t. permutations of the first k and the last N — k particles; hence 7^ and 

Jn\ are density matrices symmetric in all their variables in the sense (|1.6|) . Apart from the physical 

(k) 

densities "f N , we introduced, starting from the wave functions ipx and ipN,*, two more sets of density 
matrices; the densities 7 jy and 7^*, regularized with the cutoff x( k Hn/N), and the densities 7 jy 

and 7^, regularized with the cutoff x = xi^H^/N). Observe that, since the operator acts 
trivially on the first k variables of ipm and ipN*, we have 



®k 



XV 



®N-k 



and V^,* = {<p*T k ® 



I XV 



®N-k\ 



XV 



iN-k 



\X<P 



m-k\ 



(5.12) 



where 



k factors 



Hence 7^ = ^q"' (see (I5.3P ) and 



7Ki = i^><^r* = i^)(v*i®-®i^><^ 



k factors 

for every k > and N > k. We estimate the l.h.s. of (|5.9p by 



j(k) (, 



< 



+ 



+ 



(5.13) 
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The first term on the r.h.s. can be bounded by 

X^N XtpN,: 



^j (k) (i^-tl)\<nJ ik) \\ 

4k\\jW 



< 



1 - CkV2 



Wx^nW IIxV'at,*! 

93 — ^*11 < e /3 



<4||J«|| 



llxV> 



AT 



(5.14) 



uniformly in k, for k small enough (recall that x = x( k Hn/N)). Here we used (j5.7|) . ||x(K-ffjv/-^")|| < 
1, = ||y*|| = 1, and the choice \\(p — ip*\\ < £/(24/c||j( fe ) ||). Analogously, the third term on the 
r.h.s. of (|5.13p is bounded by 

^J {k) (7^-7^)1 = |Tr (bWr-|^r)| < 2fc[| J (fc) || ll^-^ll <e/3. (5.15) 
It remains to bound the second term on the r.h.s. of f)5. 13j) . To this end, we note that 
TrJ«(7Sf,l-7Sl)|<2||J«r X*N,. 



\\x^n,*\\ HxVw,*| 



(5.16) 



IIxV'at,. 



where we used (15. 7p . To estimate the last norm we expand the function x using the Helffer-Sjostrand 
functional calculus (see, for example, [6]). Let x be an almost analytic extension of the smooth 
function x of order two (that is |(9jx( z )I < C| j/ 1 2 , for y = Imz near zero): for example we can take 
X(z = x + iy) := (x( x ) + WX'( X ) + x" ( x )(w) 2 l^)^( x -,y)i where 9 £ C^(R 2 ) and 9(x,y) = 1 for 
z = x + iy in some complex neighborhood of the support of x- Then 



(X - X)^N,* = ~~ J &x&yd z x{z) f — ^| 



1 



(kH n /N) z _ (kH n /N) 



K 

"Ntt 



dx dydzx(z) 



1 



Taking the norm we obtain 



Ck 



II (x - x)Vw,*ll < -rr / dxd y 



Next we note that 



N 



Hn — Hn 



z-{kH n /N) 
\d-zX{z)\ 



[Hn — Hn) 



1 



(5.17) 



z-(kH n /N) 



Vat,* 



(Hn — Hn)- 



z-(kH n /N) 

k 



(5.18) 



i<k,i<j<N 



Therefore, using the symmetry of ipN* an d of Hn w.r.t. permutations of the first k and the last 
N — k particles, we obtain 



(Hn — Hn) 



1 



z-(kH n /N) 



< k 



A1 



1 



+ 



+ 



z- (kHn/N) 
V n (xi - x 2 ) 



N 

k(N - k) 



z - (kHn/N) 



N 



V N (xi - x k+ i) 



IpN,* 



1 



(5.19) 



z-(kH n /N) 
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Using Lemma |A.3| we can bound 

1 



V N (xi - Xj) 



z-{kH n /N) 



< (const.)N 3 ^ 2 \\V 2 \\ L x 



'1 - Ai 



z - (kH n /N) 



(5.20) 



Since Hn does not depend on the variable xi, we can commute the derivatives with respect to x\ 
through the resolvent (z — (kHn /iV)) _1 and we conclude that 



N 



H 



1 



z-{kH n /N) 



\y\ 



(5.21) 



for a constant C which depends on k, but is independent of N and k. Inserting this bound into (|5.18l) 
we obtain \\(x — x)*0at,*|| < CW 3 ' 3 / 2 ™ 1 and thus, since we assumed f3 < 2/3 , we have, by (|5.16p . 



TrJ^ (7^-7S,l)|< e /3 



(5.22) 



for N sufficiently large (uniformly in k). Together with (|5.14p and (|5.15p . this completes the proof 
of part hi). □ 



6 A-Priori Estimate 

The aim of this section is to prove an a-priori bound for the solution ■0^- 1 of the Schrddinger equation 
(|1.3|) with initial data ip^ ( as defined in Proposition 15. ip . Introduce the operator 

S r .= (1-A,) 1/2 . 
The a-priori bound is an estimate of the form 

(r N ,t, (1 - Ai) ... (1 - A k )r N ,t) = I dx |(1 - Ax) 1 / 2 ... (1 - A fe ) 1/2 ^, t | 2 < C k (6.1) 

for all k > 1, uniformly in t E M and in N, for A?" large enough (the constant C depends on k, but is 
independent of N, t and k). With density matrix notation, 16. 11 is equivalent to 

Tr SxS 2 . . . S k i$ t S k . . . S 2 S 1 < C k . 

To prove (|6.ip . we make use of the following energy estimate, which gives an upper bound on the 
mixed derivative operator in terms of higher powers of the Hamiltonian H^. 

Proposition 6.1 (Energy Estimate). Let Hjy be defined as in il.ty) . with V smooth and positive, 
and with < (3 < 3/5. Define 

- 1 

H N := + J2 Vn ^ ~ Xm ) = H N + N. 

j=l i^m 

Fix k G N and < C < 1. Then there is Nq = Ao(Zc) such that 

(V, S 2 S 2 . . . 5 2 0) < ^ ty, HU) (6.2) 

/or aZZ A" > Nq and all ip G L 2 (IR 3Ar ) (i/ie subspace of L 2 (M. 3N ) containing all permutation symmetric 
functions). 
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Proof. The proof of this proposition can be found in [7]. The constant 2 could be replaced by any 
constant bigger than 1 at the expense of increasing Nq. □ 

Using this energy estimate, the conservation of the energy along the time evolution, and the fact 
that at time zero, (V'am ^n^n) — C k N k by the choice of tp 1 ^ (see Proposition [57T]) , we obtain, in the 
next theorem, the bound (|6.ip . 



Theorem 6.2 (A-Priori Estimate). Let Hn be defined as in il.ty) , with V smooth and positive and 
withO < j3 < 3/5. Fix f € i7 1 (M 3 ) with \\ip\\ = 1 and k > sufficiently small, k < KodMli/ 1 , V). Let 
^tv(x) = Yij=i fi x j) andifj 1 ^ = x( k Hn /N^n /\\x( k Hn /N)i{jn\\- Suppose thatip^t * s the solution of 
the Schrddinger equation M.Sfy with initial data tp'ir. Then, for every k > 1, there exists Nq = No(k), 
with 

{r N ^sisi...sir N ,t)<c k (6.3) 

for all N > Nq, and for all t £ R. The constant C depends on the unsealed potential V(x), on the 
H 1 -norm of ip, and on the cutoff k > (actually, C is proportional to 1/k for small k), but it is 

~(k) 

independent of t £ M., N and k. If we denote by j N t the marginal densities associated with the wave 
function ipNt> then \6. 3\) is equivalent to the bound 

Tr\S 1 ...S k ^ t S k ...S 1 \<C k (6.4) 

for all N > Nq. 

Proof. For fixed k > 1, using H% < 2 k (H% + N k ), it follows from Proposition O that 

i^s! ...s 2 k r N ,t) < ^(r N ,t,H k N r N ,t) + 2 k m,t\\ 2 = ^m,H k N r N ) + 2 k (6.5) 

for all N > No(k). Here we used that the energy and the L 2 -norm are conserved along the time 
evolution. From Proposition 15.14 we obtain 

■ ■ ■ s l ^) ^ ir^k^ + 2k (6 - 6) 

which completes the proof of the theorem (assuming that < k < 1/C 2 ). □ 

7 Compactness of 7^ 

~(k) 

In this section we keep the cutoff k > fixed. We prove that, for fixed k > 1, ^y N t defines a compact 
sequence in C([0,T],Hk) (recall that 7!^! is the £>particle marginal density associated with the wave 
function ip^t)- To establish this result we prove the equicontinuity of 7^ in t G [0,T], and then 

~(k) 

we apply the Arzela-Ascoli theorem. To prove the equicontinuity of T N t we use a simple criterium, 
stated at the end of this section, in Lemma 17.21 

~(k) 

Theorem 7.1. Fix k > 1, T > 0, and n > small enough. Let j\ri be the k-particle marginal 
distribution associated with the solution ip^t °f the Schrddinger equation M.'J\) . with regularized 

initial data ip^ (see Proposition 15. 1\ for the definition of ip 1 ^). Then we have 7^ E C([0,T],H k ) 

for all N > N{k,n) large enough. Moreover, the sequence 7^ is compact in C([0, T], TLk)- If 
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7oo\ £ C([0,T],Tik) is an arbitrary limit point of "f^i, t/ien 7^ is non-negative, symmetric with 
respect to permutations (in the sense fl.8\) ) and satisfies 

\\lt\tW = Tr\S 1 ... S k ^] t S k ...S 1 \<C k (7.1) 

for all t £ [0,T] and k > 1 (the constant C is the same as in J6',^| ) ; and depends on the unsealed 
potential V{x), on the H 1 norm of tp and on the cutoff k > 0, but it is independent of k). 

Proof. We prove that the sequence 7^ is equicontinuous in t, for t £ [0, T\. To this end we define 
the following dense subset of Ak ■ 

J k := {J« G /C fe : WSiSjJ^S^Sr 1 ]] < 00, 1 < i < j < k} . 

We will prove that there exists a threshold N(k,rt) such that for every e > and for every J® £ 
there exists 5 > such that 



(*) 



sup 

N>N(k,K) 



< e 



(7.2) 



for all t,s & [0,T] with |t — s| < 5. Combining this with Lemma 17.21 below, we will obtain the 
equicontinuity. 

In order to prove (|7,2p . we use the BBGKY hierarchy (|2,ip . rewritten in the integral form 



~(fc) _~(fc) 
7/v,t T/V,s 



3=1 



£ / dr [-A.,7^] - ± £ / d - " ^tS&I 



T,( k ) 



\ / j=1 Js 



(7.3) 



where we recall the notation Vn(x) = N^V(N@x). Multiplying last equation with and taking 
the trace we get the bound 

Tr J (fe) (tS*1 - 7^) I < £ /' dr |Tr (s^ J«S, - 5, jWfl^) S/^i 



-SiSjJ^Sr^^SiSj^SiSjSr^VNixi - x^Sj 1 
—SjJ^Sj 1 SjSk+i^^ T Sk+iSjSj 1 S k ^ 1 V]y(xj — Xk+ijSfr^Sj 1 



(7.4) 

using that Sk+i commutes with JW. Using that \\S^ 1 Sj l VN{xi — Xj)S^ l Sj l \\ is finite, uniformly 
in N (see Lemma [A. 3p . and the assumption that SiSjJ^S^ l Sj l is bounded, for every i,j < k, we 
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find 



TrjWffil-i/p) <2k\t-s\ S np\\Sr 1 J^S J \\ sup IV {S^Sjl 

V ' 3<k j<k,re[s,t] 

+ k 2 N~ l \t-s\ sup WSiSjJ^S^Sr 1 \\ sup Tr \S i S j ^ T S j S l 

i,j<k i<j<k,r<=[s,t] 



+ 2k(l- 4) sup WS^jWSjW sup TV ISjSk+a^Sk+^l 

iV ' j<k j<k,r£[s,t] 

(7.5) 



and thus, by Theorem! 

™* k) C k \t-s\ (7.6) 

for a constant Ck depending on k and (but independent of t, s and N) and for all N large enough 
(depending on k and on the cutoff n > 0). This implies (|7.2p and, by Lemma [7.21 it implies that the 
sequence 7^ G C([0, T], 7ifc) is equicontinuous in t (with respect to the metric p& defined on 

Since moreover the sequence y^ t ^ s uniformly bounded in 7ik (for sufficiently large, by Theorem 
note that here n > is fixed), it follows by the Arzela-Ascoli Theorem that it is compact. 

To prove that an arbitrary limit point 7^ of the non-negative sequence 7^ is also non-negative, 



we observe that, for an arbitrary tp G L 2 (R 3fc ) with \\ip\\ = 1, the orthogonal projection P v = \tp)(ip\ 
is in Ak and therefore we have 

(ip, 7« <p) = TV P^t\ = lim Tr P^f t = lim (<p, j$ t <p) > , (7.7) 

for an appropriate subsequence Nj with Nj — > 00 as j — * 00. Similarly, the symmetry of 7^ w.r.t. 

permutations is inherited from the symmetry of for finite N. In fact, for an arbitrary G Ak 
and a permutation tt G Sfc, we have 



Tr J«7S = lim = lim TV ^6^?^. = lim Tr 6.-, J«9,f t 



(7.8) 



IOO,t X1 J W T /oo,t 

where we used that, since G -4fc, also 0„—i J^@ n G .Afc, because 

||e,-ij( fc )e^iu fe = ||5 1 - 1 ...5 fe 1 e 7r -ij( fc )e 7r ^- 1 ...5 1 - 1 || 

= lie^-xSf 1 . . . s k 1 j^s^ 1 . . . sr 1 ^ II ( 7 - 9 ) 
= \\s^...s- 1 j^s- 1 ...s^\\ = \\j^\\ Ak . 

Finally, the bound (|7.ip follows because in the weak limit the norm can only decrease. □ 

(k) 

Lemma 7.2. Fix k. A sequence of time- dependent density matrices j N t , N = 1,2, . . ., defined for 
t G [0, T] and satisfying 

sup sup hP t \\n k < C , (7.10) 

N>lte[0,T] 

is equicontinuous in C([0, T], Hk) with respect to the metric pk (defined in h3. 3\) ). if and only if there 
exists a dense subset Jk of Ak such that for any G Jk and for every e > there exists a 5 > 
such that 



sup 

N>1 



for all t,s G [0, T] with \t — s\ < 5. 
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Proof. The proof of this lemma is similar to the proof of Lemma 9.2 in [Sj; the main difference is 
that here we keep k fixed, while in [9] we considered equicontinuity in the direct sum C([0,T],7i) = 
©fc>lC([0,r],H fc ) over all k > 1. □ 



8 Convergence to solutions of the Gross-Pitaevskii hierarchy 

From Theorem 17.11 and from the Cantor diagonalization argument explained in Step 2 of the proof 
of Theorem II. 1^ we know that the the sequence Tnj = {l]^\} k= i nas a ^ least ° ne limit point 

Too,* = {7^t}fc>i m C([0, T},Ti) with respect to the r-topology. In the next theorem, we show that 
any such limit point is a solution of the infinite Gross-Pitaevskii hierarchy (12. 2\ in the integral form 
(|4.4p . The analogous theorem from [7] cannot be directly applied since here we work in R 3 in contrast 
to the compact configuration space of [7j. Moreover, the infinite hierarchy (14. 4p is defined somewhat 
differently than (1.8) from [TJ. 

Theorem 8.1. Assume Hn is defined as in ( EH| ), with < [3 < 1/2. For a fixed k > 0, let ip% t be 
the solution of the Schrddinger equation il.3\) . with initial data V^r ( defined as in Proposition \5. 1}) . 

and let FV,* = {lNUk=l ^ e mar 9^ na ^ densities associated with ip^ t - Suppose To^t = {7^| t }fc>i £ 
C([0, T],Tt) is a limit point of the sequence r^t with respect to the T-topology. Then is a 

solution of the infinite hierarchy 

7S = ^ k \t)lZ - ibo £ f ds " «) Tr k+1 [S( Xj - x fc+1 ),^+ 1} ] , (8.1) 

3=1 J ° 

with initial data 

k 

/or all k > 1. 

The action of the delta-function in the second term on the r.h.s. of (|8.ip is defined through 
a limiting procedure. We define the operator B^ k \ acting on densities ^ k+l ^ with smooth kernel, 
7 ( fc + 1 )(x fc+1 ;x' fc+1 )G5(M 6 ( fe+1 ))by 

k 

S (*) 7 (fc+i) = _ i6o £ Trfe+1 [5{x . _ x fc+1 ), 7 (*+D] . (8.3) 
i=i 

If we interpret this definition formally for arbitrary density matrices, then the infinite hierarchy (18. ip 
can be rewritten in the more compact form 

7S = l4 k) (*)t£!o + f d« (* - s)B^^ . (8.4) 

JO 

The action of B^ on kernels is formally given by 
(5 (fc) 7 (fc+1) )(x fc ;x' fc ) = -ib ^2 J dx k+1 (5(xj - x k+l ) - 5(x f j - x k+1 )) 7 (fc+1) (x fc , x k+1 ; x' k , x k+1 ) . 

3=1 

(8.5) 
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For the more precise definition of B^ k \ we choose a positive smooth function h £ C°°(IR 3 ), with 
compact support and such that J dx h(x) = 1. For a > 0, we put 5 a (x) = a~ s h(a^ 1 x). Then, for 
7 (fc+1) 6 Hfc+i, we put 

fB(fe) 7 (*+i)"\ (x fc ;x' fc ) :=-»6o lim V / dx k+1 dx' k+1 5 a2 (x k+1 - x' k+l ) 

\ J ai,a2-+0~J (8.6) 

x (5 ai (xj - x k +i) - 6 ai (x'j - Xk+i)) -/ ( - k+1 \^ k ,x k+1 ;x.' k ,x' k+1 ). 
Lemma 18.21 below will show that is well defined for any ^ k+1 ^ g l~tk+i- We introduce the norm 
sup (X!) 4 . . . (zfc) Vi> 4 • • • <4> 4 f | J (fc) (x fe ;x' fc )| + |V-. J (fe) (x fe ;x',)| + \V X , j( fc )(x fc ; 4 



||j .— sup ...yxfc/vx!/ ...\^ k / I I v & J / 1 I V & 5 A: / 1 I * 2?'- v ™ ' / 

x fc ,x' fe 

(8-7) 

for any j <k and for any function j( fc )(xfc; x^). 

Lemma 8.2. Suppose that 8 a {x) is a function satisfying < 5 a (x) < Ca _3 l(|x| < a) and 
f S a (x)dx = 1 (for example 6 a (x) = a~ 3 h(x/a), for a bounded probability density h(x) supported in 
{x : \x\ < 1}). Then if ^ k+1 \-x. k+ i; x' fc+1 ) is the kernel of a density matrix on L 2 (M. 3( - k+1 ^), we have, 
for any j < k, 

J dx fc+ idx' fc+1 J (fc) (x fc ;x' fc ) (5 ai (x k+1 - x k+1 )5 a2 (xj - x k +i) - 6(x k+1 - x' k+1 )S(xj - x k+1 )) 

x 7 (fc+1) (x fc+ i;4 + i) 

< (const.)* p<-% ( ai + v^) Tr|^ +l7 (fc+1) SA+il . (8.8) 

Recall here that Si = (1 — A^) 1 / 2 . Exactly the same bound holds if xj is replaced with x'a in 118. 8\) 
by symmetry. 

This lemma is similar to Proposition 8.1 in [9], with the difference that here we work in the 
infinite space M 3 instead of a compact set A as in [9]. For completeness we give a proof of Lemma 
8.21 at the end of Appendix [A] 

It follows from the this lemma that the limit (|8.6p exists for ^(' c+1 ) £ 'Hk+ii i n an appropriate 
weak topology, and that it is independent of the choice of h G C°°(]R 3 ). Here, with a slight abuse 
of the notation, TL k+ i is used both for the space of densities defined in Section [3J and for the space 
of kernels associated with these densities. Hence the operator B^ k \ originally defined on Schwarz 
functions, can be extended to a bounded operator from the whole 7i k+ i and with values in some 
sufficiently large Banach space determined by the conditions on the test function j( k \ Moreover, 
the following bound holds 



J dx fc dxU«(x fc ;x / fc )Tr fe+1 [5(x j -x fc+1 ), 7 (fc+1) ] 



< C fe |||j( fc )||L-Tr|5,5 fc+1 7 (fe+1) 5 fc+1 5 J | (8.9) 



for each term in (|8.3p . therefore a similar bound holds for the operator as well. 

The equality in (|8.4p is then interpreted in the sense that there exists a representation of 

ds4 k \t-s)B( k h£t 1] 



o 



( k \ ( k \ ( k \ 

which lies in TC k and such that (|8.4j) holds. This follows from the fact that both 7^ and Uq (i) 7 oo 

are in 7i k and the equality can be checked in a weak sense. 
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Proof of Theorem \8.1\ Without loss of generality we can assume that rVt converges to with 
respect to the T-topology. This implies that, for every fixed k > 1 and t E [0, T] we have 

~(fc) ~(k) 

^N,t ~ * 7oo,t 

with respect to the weak* topology of That is, for every E we have 



TrJ^ (7^-7^) 



5.10) 



for A?" — > 00. 
Let 



:= n c^j) + ^) ( with ^ = c 1 - a %) i/2 ) • 

In the following we assume that the observable E KLk C ^ is such that 



HS 



< OO, 



Ul) 



where ||A||hs denotes the Hilbert-Schmidt norm of the operator A, that is ||-A||jjs = Tn4*A Note 
that the set of observables satisfying the condition (|8.1ip is a dense subset of Ak- 
It is straightforward to check that 



\S x ...S h J (fc) || < n 7 k J^n 7 k , and || J^S* . . . SJ < O^J^fi 

HS 



HS 



Moreover, for any j < k 



|||J (fc) lb-<(const.) fc nlj^nl 



HS 



where the norm |||.|||j is defined in (|8.7p . This follows from the standard Sobolev inequality 
(const.) ||/||2,2 in three dimensions applied to each variable separately in the form 



1.12) 
5.13) 



(sup (s) 4 (x / ) 4 |V se J(a;,x / )|) < (const.) / dxdx' (1 - A x ) (x) 4 (V x J(x, x')) (x'Y 

^ x,x' ' J 

< (const.) Tr (1 - A)(x) 4 V J (x) 8 J* V* (x) 4 (1 - A) 

< (const.) Tr ft 7 Jft 14 J*ft 7 

with O = (x) + (1 — A) 1 / 2 . Similar estimates are valid for each term in the definition of ||| • for 
j < k. Here we commuted derivatives and the weights (x); the commutators can be estimated using 
Schwarz inequalities. 

Rewriting the BBGKY hierarchy (|2.ip in integral form, and multiplying it with we obtain 



TV( fe ) = Tr« jt k W k \ttfjp 



i 1 



k 
N 



V / d S Tr« J( fc ^( fc )(t - s)Tr fc+1 [V N (x 3 - x k+1 ),^\ 
j=i Jo 



U4) 



where we recall the notation Vn(x) = N 3 ^V(N^x) and 



5.15) 
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with 

k 



j=l i<j 

Here we use the notation Tr^ instead of Tr to explicitly stress that we take the trace over the 
degrees of freedom of k particles. 

The l.h.s. of (|8T4"P clearly converges, as N -> oo, to Tr^ J {k) l^] t (by (18301 and because 
G /Cfc C ^4fc by assumption). As for the first term on the r.h.s. of (|8.14p we have 

Tr« J« (u^(t)j% - U^\t)^%) - (8.16) 
for N — > oo. The definition oilA^ is recalled from (|4.3|) . To prove (|8.16p we note that 
TrW (W«(t)^, -^(t)7«) = J W 7$ 



+ Tr« J(^f(t) (7S-7go 



.17) 



The second term converges to zero, for N — > oo, because, if G A, then also U^\-t)jW G .4 fe , 
and hence 

Tr J«l#0(t) (>£ " 1%) = ^ J«) (# - I®,) - 



as iV — > oo. As for the first term on the r.h.s. of (|8.17p we have 

TrW JW ( W «(t) - Z#>(t)) 7^o = ^ £ f d, Tr« - " ^i)^ WtSJJ, • 

(8.18) 



Kj 

This implies that 

A; 2 t||jW 



TrWjW (^)(t)-^(t))7? < " sup \\V N ( Xi - Xj )S.'S. l \\ sup Tr 



ATl-(3/3/2) 



Tr ISiSWVoS^Sil , 



where we used Tr \SiSjj\ < Tr \SiSj~/SjSi\, the permutation symmetry of and the bound 
|| V}v(a^i — x j)S^ 1 S~ 1 \\ < CA 3/3 / 2 with a constant C that only depends on the unsealed potential 
V{x) (see Lemma fA.3p . Since /3<l/2<2/3, we get, from Theorem 16. 2[ 



Tr(fc) j(fc) ^(fc)( t ) _ U^it)) 7$ I -» (8.19) 
as A — > 00. 

Next we consider the second term on the r.h.s. of (18.14p . More precisely we prove that the 
difference 



3=1 



V / dsTrW - s)Tr k+1 [V N ( Xj - x k+1 ),^ + s 1] ] 

~~i J 

-60E/ d S TV( fc )j( fc )^ fe) (i-s)Tr fe+1 [5(x, - x fe+1 ), (8.20) 

7 = 1 ^ 
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converges to zero, as N — > oo. To this end we write this difference as the following sum of four terms 
k . 

-^E / d S TrWj(%«(t- S )Tr fc+1 [V N { Xj - x k+1 ),^f] 

3=1 ° 

+ J2 f'dsTrW JW - s) - W W (i - *)) Tr fc+1 [^( 

3=1 
k . 

+ V / ds Tr( fe ) j(% (fe) (t - s)Tr fc+1 [(^(^ - s fc +i) - - x k+1 )) ,7^' 

■/ 



x ~(fe + l)l 



k rt 
10 



+ b J2 I dsTr^ J^U ( k) (t - s)Tt k+1 \S( Xj - x k+1 ), 

3=1' " 



~(*+l) _ ~(fc+l) 
IN.s 100,8 



and we prove that each one of these terms converges to zero when TV" — > 00. 

Using that S&+1 commutes with J^ k >U^ k \ the first term can be bounded in absolute value by 



k . 

d^Tr^ 1 )^^-*)^ 

3=1 

2kH\\J^\\ no _ 1 



< 



N 



S kliVN{xj - Xk+^S^W sup TV Sk+xj^a S k+ i < CN~ 

se[o,t] 



(8.22) 



as TV — > 00. Here we used Theorem 16.21 and that WS^^VN^Xj — Xk+i)S^} 1 \\ < CN@, for a constant C 
which only depends on the unsealed potential V(x) (by the second statement in Lemma lA.3j) . The 
constant C on the r.h.s. also depends on the cutoff K. 
The second term in (I8.2ip . can be rewritten as 



I dsTr^ J« (uW(t - s)-U { *\t - sj) Tr fc+1 [V N ( Xj - x k+1 )tf£ 
3=1 

= -ijy- 1 fr f ds r S dr Tr( fc+1 ) J« W W (i - s - r) 
T~i 77L Jo Jo 



5.23) 



j=l ^<m 

x W N {x t - x m )M {k) {T) V N (xj - i fc+1 ),T^+ 13 

Expanding the two commutators, we find that the absolute value of the r.h.s. of the last equation 
can be estimated by 



k k 



CT-TV f ds f "dr (\\J {k) S e S m \\ + \\S e S m J^\\) \\Sj l S^V N {x, - x r 
j= ie< m Jo Jo 
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x \\S k liVN(xj - x k+1 )S k l 1 S j 1 \\ Tr \SjSk+ry^gSk+i\ ■ 



Notice that ||5 , fc _^ 1 Vjv(xj — Xk+i)S k ^Sj 1 || < CN 13 ^ 2 using both inequalities from Lemma [A. 31 and a 
Schwarz estimate. Combining this with (|S32]> . with ||5r 1 5~ 1 Viv(a;f - z m )|| < CiV 3/3 / 2 and 



Tr |5jSVfi7jv,s - ^ £i#27jv,s <S2<Si 



5.25) 
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for all j = 1, ... , k, by the symmetry of , we find 



(fc+l)- 

s 



< CN- 



-1+2/3 



1.26) 



k ,.i 

| / d ^ Tr(fc) J (k) {U {k) (t - s)-U ( k \t - s)) Tr k+1 [V N (x 3 - x k+1 )tf£ 

which converges to zero, as N —* oo, for (3 < 1/2. We remark that this is the only step where the 
more restrictive (3 < 1/2 condition is used, the rest of the proof works for (3 < 3/5. 

Next we consider the third term in (|8.2ip . Using the kernel representation of J^j~t , the absolute 
value of this term can be estimated by Lemma 18.21 (with a\ = and 02 = N~@) as 

b J2 J* ds j / dx fe dx' fc dx fc+1 (U { k) (s - t) J«) (x fe ; x' fc ) 

x \J^ V n( x 3 ~ x k+i) -S(xj -Xk+i) \ ^^(xfc.Xfe+ijXfc.Xfc+i) 
j dx fc dx' fc dx fc+1 (U^ k \s - i)J (fc) ) (x fc ;x' fe ) 



5.27) 



x (—V N (x' j -Xfc+i) -<5(4 -x fe+ i)J 7^ 1) (x fe ,x fc+ i;x / fc ,x fc+ i) 

<CN~M 2 sup TrlSjSwTfi+VSwSjl V / ds ||^ fc) ( S - t) J^l,- , 
j'<Jfe,«e[o,t] J=1 jo 

for a constant C depending on /c. Here we recall the definition of the norm ||| • \\j from (|8.7p . Using 
the estimate (|8.13|) . we have 



HS 



i(s-i)p1/ v .\m-i(s-t) p 2 



with a /c-dependent constant C . Since e ( ' i (xj) m e 
1, . . . , k, m & N, with j?j = — iVj, we obtain that 

#(*) 



(xj + 2(s — t)pj) m , for any j 



HS 



From the assumption (|8.1ip on and the a-priori control on , we obtain that the r.h.s. of 

(|8.27p is bounded by CN~^/ 2 with a constant C depending on k, on i, on JW, and on the cutoff At. 
Hence, the third term on the r.h.s. of (|8.21|) converges to zero, as N — > 00. 

Finally, we consider the fourth term in (|8.2ip . For fixed s £ [0, t], and j < k we have 



Tr( fe +D jW Wo (i - s)8( Xj - x k+1 ) " 7^ 1} 

for N — > 00, because J^Uo(t — s)5(xj — x k+ \) £ A k +i- In fact 

\\S? . . . S^ v J^U (t - s)5(x j - Sfc+i)^! . • • S^\\ 

< WJ^SjUS^S^Sixj -Xk+JS^S. 



3.28) 



-1 q-ll 



.29) 
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is finite. It follows that the integrand in the fourth term in (|8.2ip converges to zero, for every 
s G [0, t], and every j < k. Since the integrand is uniformly bounded using the uniform (in s) apriori 
estimates on and 7^1^, and the uniformity of the .Afc+i-norm of J^Uo(t — s)5(xj — x k+ i), 

it follows that the fourth term converges to zero as well, for N — * oo. This proves that, for every 
t G [0,T], k>l and G K k satisfying (fHTTTD . we have 



Tr W j(k)~(k) t = Tr (k) j(^) (t) ^) Q 

- ib f^ r^TrW J^ul k \t - s)Tr k+1 [5( Xj - x k+1 ),^] 
,-1 Jo 



1.30) 



i=i 

This implies that 

lt]t = 4 k \t)l% ~ ibo Y, f dsU t\t ~ s)Tr k+1 [S( Xj - x k+1 ),^} (8.31) 

j=1 Jo 

~(k) 

if we consider 7^ as elements of a large space of density matrices, the dual space of the Banach 

space consisting of all sufficiently smooth (such that satisfies (|8.1ip ). Next since 7^ t G 7i k 

and Wq (i)7^j £ Ti-ki it follows that also the second term on the r.h.s. of (|8.31|) lies in Tt k (or at 
least it has a representation as element of Tt k ), and that (|8.3ip holds as an equality on TC k . 

Finally we prove (|8.2p . For arbitrary N > k and G JC k , we have 

Tr J<*> (tSo " 7«) = Tr J« (t®, - 7$) + * ^ (^0 - 7? } ) • (8-32) 

From (|8.10p (with t = 0), the first term converges to zero, for N —* 00. The second term converges 
to zero, as TV" — ► 00, by Proposition 15.14 part iii). □ 

9 Uniqueness of the infinite hierarchy 

In this section we show the uniqueness of the solution of the infinite hierarchy (|2.2[) . The following 
theorem is the main result of this section. 

Theorem 9.1. Fix T > and 60 > 0. Suppose Tq = {70 }fc>i is such that 7q is non-negative and 
symmetric with respect to permutations (in the sense of i 1.8]) ) and it satisfies 

\h { k) \\n k = Tr\Sx . . . S k7 { Q k) S k . . . S 1 \ < C k (9.1) 
for all k > 1 with some constant C . Then the infinite hierarchy 

^ =U^\t) 1 {k) -ib^ f dsUi k \t-s)Tr k+1 [S(xj-x k+1 ), 7 ^], (9.2) 

pi Jo 

has at most one solutionTt = {7+ }fe>i £ C([0,T],TI) withTt=o = Tq, such that 7+ is non-negative, 
symmetric with respect to permutations and satisfies the bound 



hi h) \\n k <C k (9.3) 



for all k > 1, and t G [0, T}. 
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Remark: In the proof we set bo = 1 for convenience. The inclusion of 60 modifies all bounds in a 
trivial way, but it plays no role in the argument. 

In order to prove this theorem, we will expand the solution in a Duhamel-type series. Recall from 
Section [8l the formal definition of the operator B^ k \ given by 

k 

B (fc) 7 (fc+i) = -i^Tr k+1 [5( Xj -z fc+ i),7 ( * +1) ]. (M 
j=l 

On kernels in momentum space B^ acts according to 

(B {k h {k+1) )(Pk;p' k ) = / d *+idfe+i (i {k+1) (pi,-,Pj -qk+i + q'k+i,..,Pk,qk+i;Pk,qk+i) 

-7 (fc+1) ( Pfc , ...,p'j + q k +i - q'k+i, ■■■,p'ki q'k+ii) 

k „ Ik \ 

= (-t)X; / dq fe+ idq' fc+1 J] «J(p, - % )<^ - $ 7 ( fc+1 >(q fc +i; q' fe+ i) 

x [5(pj- - q'j)S(pj - (qj + qk+i - q'k+i)) - <% - qj)5{p'j - (q'j + q'k+i ~ ■ 

(9.5) 

These definitions are formal: they can be made precise using Lemma 18.21 as explained in Section [HJ 
In the current paper we will work in momentum space, i.e. we apply (|9.5p repeatedly and we will 
show that all integrals are absolute convergent. 

With these notations we can expand the solution {7^ } of (|9.2p for any n > 1 as 



7? } = ^ k \t)^ + ^ / dsi r d S2 . . . d Sm 

m=l J ° J ® J ° 

x W (fc) (t - Sl )B(% (fe+1) ( Sl - S2 )i?( fc+1 ) . . . J B( fc+m - 1 ^ fc+m) ( Sm ) 7 J fc+m) 



ft /"Si 

+ / dsi / ds 2 ... I ds n 



(9.6) 







x &f (i - Sl )S(%( fc+1) ( Sl - , 2 )i?( fc+1 ) . . .^ fc+n " 1 >( an _ 1 - , n )B<*^> 7 g-H0 . 



The terms in the summation will be called fully expanded as they contain only the initial data. The 
last error term involves the density matrix at an intermediate time 5^. 

In Sections 19. II and 19. 21 below we show how the terms in this expansion can be written as a sum of 
contributions of suitable Feynman graphs. In Section 19.31 we show how to bound the contributions 
of the Feynman graphs. Then, in Section 19.61 we use these bounds to prove Theorem 19.11 Some 
technical estimates, used in Section ^. 3l to bound the contributions of the Feynman graphs, are shown 
in Section [TUJ 

Notation. For the rest of this paper we will mostly work in Fourier (momentum) space. We use 
the convention that variables p, q, r always refer to three dimensional Fourier variables, while x, y, z 

are reserved for configuration space variables. With this convention, the usual hat indicating the 

(2) 

Fourier transform will be omitted. For example, the kernel of a two-particle density matrix 7 q ; in 
position space is 73 (xi, X2; x%, x 2 ); in momentum space it is given by the Fourier transform 

7 f(<Zl,<Z2;<zU 2 ) = / dsidsadarida^Csi,^ , 
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with the slight abuse of notation of omitting the hat on left hand side. 

Furthermore, to avoid (27r)-factors in the Fourier transform, we make the convention that the 
integration measure for the three dimensional momentum variables p, q, r are always divided by (27r) 3 , 
i.e. 

dp := ; Lcb ^ for all three dimensional momentum variables 
1 (2vr) 3 

where d Lob denotes the usual Lebesgue measure. We will also use delta functions in momentum space, 
5(p), and they will correspond to the measure dp above, i.e. 

j f(p)5(p-Q)dp = f(q) 

for smooth functions. Delta functions in position space, <5(x), remain subordinated to the usual 
Lebesgue measure. 

Similar convention is used for the frequency variables (dual variables to the time) that will always 
be denoted by a: 

da := Lcba for all one dimensional frequency variables 
and to the delta functions involving a-variables. 

9.1 Graphs 

Graphical representation is a convenient tool to bookkeep various terms in the perturbation expan- 
sion of many body systems; the graphs encode the collision histories of the particles in a concise 
form. We start by introducing classical graphs, which would appear in a diagrammatic expansion 
of the solution of the BBGKY hierarchy of a classical many particle system. Afterwards we will 
define quantum graphs, which are suitable for the diagrammatic representation of the expansion of a 
quantum system like (|9.6p . The quantum nature of the expansion requires a doubling of each edge 
of the classical graphs, corresponding to the fact that in the density matrix description each particle 
is associated with two different variables (xj and x'-). In this paper we will only use the quantum 
graphs. Classical graphs are included only to facilitate the description of the quantum graphs and 
to derive a combinatorial estimate on their numbers. 

9.1.1 Classical graphs 

We consider rooted binary trees on n (internal) vertices and n + 1 leaves. The root and the leaves are 
not considered vertices; instead we identify them with the unique edge they are adjacent to. These 
edges are called external and we will use the terminology root and leaves for the external edges. For 
n = 0, i.e. when there is no vertex, there is only one single edge, that is the root and the single 
leaf at the same time. However, when counting the external edges, we will count this edge twice. 
This tree will be called trivial. At every vertex three edges meet; the one that is closest to the root 
is called father-edge, the other two are called son- edges of this vertex. At every internal vertex we 
mark one of the son-edges. (In the expansion this son-edge will inherit the father's identity.) For 
illustration, one can draw such a graph so that the marked son-edge goes straight through, and the 
unmarked edge joins from below (see Fig. [I]). The set of marked binary rooted trees of n vertices is 
denoted by Q n . 

Note that the number of marked binary trees with n vertices, C n := \Q n \, is given by the recursion 

n 

C n +i = ^2 CkC n -k, Co = 1 . 

fc=0 
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Figure 1: Example of a rooted, marked, binary tree with n = 7 vertices 



To show this formula, note that removing the vertex closest to the root splits any marked binary 
tree into two smaller ones, with k and n — k vertices. This recursion defines the so-called Catalan 
numbers. They are given by the closed formula 



and can be estimated by C n < 4 n . 

For any G G Q n , the set of vertices is denoted by V(G) and the set of edges is denoted by E{G). 
The root is denoted by R = R(G), the set of n + 1 leaves is denoted by L = L(G). They together 
form the set of external edges, Ext(G) := R(G) U L(G). The set of internal edges is defined as 
Int(G) := E{G) \ Ext(G). We will draw the graphs as in Fig. [TJ i.e., the root is on the left, and 
the leaves are on the right of the graph. Note that the vertices V{G) are partially ordered by their 
succession towards the root: for any v, v' £ V(G) we have v -< v' if v lies on the (unique) route 
from v' to R. For any G S Q n , we denote by 0(G) the set of complete orderings of the vertices 
V(G) that are compatible with the partial order -< of V(G). In general, for a given G E Q n , there 
are several complete orderings which are compatible with the partial order of G. The complete 
ordering can be visualized by drawing the graph in such a way that the horizontal coordinates of 
the vertices correspond to the ordering (see Fig. [2]). Two rooted, binary, marked, and completely 
ordered trees G\ and G2 are said to be equivalent if there exists a one-to-one map between the edges 
and the vertices of G\ and G2, such that all adjacency relations, all marks and all labels of the 
ordered vertices are preserved. The total number of inequivalent rooted, binary, marked trees with 
n completely ordered vertices is n!, i.e. 



This follows from the fact that one can build up such a graph on n vertices by successively adjoining 
new leaves: a new leaf will be joined to an existing leaf in such a way that the existing leaf keeps 
its "father" -identity. In a graph with k — 1 vertices and k leaves there are exactly k possibilities to 
adjoin the new leaf: one can create a new vertex on any of the existing leaves and adjoin the new 
(k + l)-th leaf to it as unmarked. 

For k > 1 and n > 1 consider forests consisting of k rooted, marked, binary trees, Gi, G2, ■ ■ ■ , G/-, 
so that the total number of (internal) vertices is n. We assume that the trees, i.e. their roots are 
labelled, i.e. the permutation of the trees results in inequivalent forests. The set of such forests will 




n! . 



Geg, 
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a) Partially ordered graph 



b) Ordered graph: v < z < w 



c) Ordered graph: v < w < z 



Figure 2: Ordering of graphs in Q. t 



be denoted by G n ,k- I n Fig. El we draw an example of a graph in £7 7)4 . Note that the number of 
inequivalent forests in Q n ^ is given by 

(ni,...,n fc ):^ = i «i=n 

where the summation runs over all ^-tuples of nonnegative integers that add up to n. This number 
can be bounded by 

n + k - 1\ „ „ 3n+A . 



Again, for G S G n ,ki we wm denote by V(G) the set of the vertices of G, by E(G) the set of 
edges, by Int(G), and Ext(G) the sets of internal and, respectively, external edges. Moreover, R(G) 
and L(G) will be used for the set of roots (there are k roots in each forest), and for the set of leaves 
(there are n + k leaves) of G, respectively. The vertices in V(G) are again partially ordered by their 
succession towards the roots: for any v, v' £ V{G) within the same tree we have v -< v' if v lies on 
the (unique) route from v' to the root of the tree. There is no order relation between vertices in 
different trees. For a given G £ G n ,k, we define 0(G) as the set of complete ordering of the n vertices 
of G which are compatible with the partial order of G: the number of non-equivalent forests with k 
trees and n completely ordered vertices is then 

( n 4. h — i V 
k(k + l)(k + 2)...(n + k-l) -- ' 



(fc-1) 



This can again be seen by the successive build-up of the forest; if one starts with k trivial trees with 
no vertices, then the first new leaf can be adjoined in k different ways, the second one in (k + 1) ways 
etc. 



9.1.2 Quantum (Feynman) Graphs 

For any n > 0, k > 1, we now define the set of Feynman graphs !F n k- An element V £ J~ n ,k is 
a union of k labelled pairs (Tj,Tj), where, for every j = l,...,k, Tj and T'- are two disconnected 
oriented marked rooted trees. As in the classical graphs, the roots and the leaves of the trees are 
not considered vertices; we identify them with the edges they are adjacent to. These edges are called 
external. Each edge in Tj and Tj is oriented (indicated by an arrow) and the edge is called inward or 
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Figure 3: Example of forest in Q n ^ with n = 7 and k = 4 



outward according to whether the orientation points towards or away from the root. For each j the 
root of Tj is outward, the root of Tj is inward. Each vertex is adjacent to four edges. At each vertex 
we require that precisely two edges are incoming and two edges are outgoing. Note that the concept 
of incoming/outgoing is relative to the adjacent vertex and it does not coincide with the concept of 
inward/outward that is solely a property of the edge. 

Similarly to the classical graphs, at each vertex we can introduce the notion of father-edge and 
son-edges. The father-edge is on the route from the vertex to the root, the other three edges are the 
son-edges. Moreover, two son-edges have the same orientation as the father-edge and one of them is 
distinguished and will be called the marked son-edge (we will say that the marked son-edge inherits 
the identity of the father-edge). 

If the number of vertices of Tj and Tj is rrij and m f j, respectively, then the number of external 
edges is 2m j + 2 and 2m! j + 2, respectively. Recall that for a trivial tree we count the leaf and the 

root as separate external edges. We assume that the total number of vertices is Ylj=i( m j + m j) = n - 
The vertices of T are denoted by V(T). The edges of T are denoted by E(T), the internal edges are 
Int(r) and the external ones are Ext(r). For e G E(T) and v G V(r), the notation e G v indicates 
that the edge e is adjacent to the vertex v. We denote by R(T) the set of the root edges and by L(T) 
the set of leaf edges. Similarly to the classical graphs, the fact that each component of T has a root 
induces a partial ordering among the vertices of T, this will be denoted by -<. Analogously, one can 
define a partial ordering among the edges of T, which will also be denoted by -<. 

Note that, because of the mark at every vertex, there is a natural pairing of the leaves of T. 
To define the pairing, we introduce the notion of the ancestor a{£) of a leaf I as follows: if I is an 
unmarked son-edge, then the ancestor of I is defined as I itself, a(£) = I. On the other hand, if 
£ is a marked son-edge, or if it is the root, then the ancestor a(£) is defined as the minimal edge 
e (minimal with respect to the partial order -<) on the route from i to the root, for which the set 
{e G E(T) : e < e < £} contains only marked son-edges. If we say that at each vertex the marked 
son-edge inherits the identity of the father-edge, then the ancestor a(£) is defined as the closest edge 
to the root whose identity is inherited by the leaf I. Clearly, the ancestor-map a : L(T) — > E(F) is 
injective (any two leaves have different ancestors) and it maps leaves of Tj or Tj into edges of Tj or 
Tj, respectively (the ancestor of £ lies in the same tree as £). Moreover, for every j = 1, . . . , k, the 
root of Tj (and of T-) is always the ancestor of exactly one leaf in Tj (or T', respectively). 
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Using the concept of ancestor, we define next the pairing of the leaves of T. For fixed j, two 
leaves £ and £' in Tj (or in Tj) are paired (we say they are companion leaves) if a(£) and a(£') are 
the unmarked son-edges of a vertex in Tj (or in Tj, respectively). Moreover, if the ancestor of a leaf 
£ is the root of Tj, then we pair £ with the unique leaf of Tj whose ancestor is the root of Tj. This 
completely defines a pairing of the leaves of Tj U Tj, and thus of the whole graph T. Note that, if Tj 
has rrij vertices, then the number of leaves of Tj is 2m j + 1: 2m j leaves are paired within each other 
and one leaf, the one with the identity of the root, is paired with the leaf of Tj whose ancestor is the 
root of Tj. 

For a given graph T G T n ,k-> we define the labelling map, tt% : R(T) —►{!.,..., k}, where ni(r) = j 
if r is the root of Tj or Tj. For e G E(T), we also introduce the notation r e = 1 if e is outward 
and r e = — 1 if e is inward. A root and the corresponding component will be called trivial if it 
contains no (internal) vertex. For T G J~ n ,k we denote by R\(T) the set of trivial roots, and we set 
R 2 (T) := fl(T)\J2i(r). Let fci(r) = \Ri(T)\ and fc 2 (r) = |i2 2 ( r )l = 2k - ki(T). Moreover, we define 
by E 2 (T) := E(T)\Rx(T), then \E 2 {T)\ = 2k + 3n - kx(T). Let Li(T) be the set of leaves of the 
trivial components; this set is naturally identified with Ri(T). Let L 2 (T) := L(T) \ Li(T) be the set 
of leaves of the non-trivial components, clearly |X,2(r)| = 2n + 2k — k±(T). 

When we draw graphs in T ni k, for every j = 1, . . . , k, we superpose the two trees Tj and Tj; this 
makes the pairing of the leaves clearer. As in the classical graphs, at each vertex we draw the marked 
son-edge so that it goes straight through, while the two unmarked son-edges join from below. 

Next we will show that Feynman graphs can be constructed in a natural way starting from the 
classical ones. Later we will demonstrate that all quantum graphs can be obtained in this way. As 
a by-product we will derive a bound for the number of elements in J- n ^. 

For any G G Q n< k and a = {a v G {±1} : v G V(G)} we define a Feynman graph, T = T(G,a) G 
T n fa as follows: We double each edge of G and equip them with an opposite orientation (arrows). 
At any vertex of G we define the new vertex of the six edges of V involved as follows. For o~ v = 1, 
the outward father-edge is joined with both edges of the unmarked son-edge and with the outward 
edge of the marked son-edge; this creates a vertex of V with four edges. The inward father-edge is 
joined with the inward edge of the marked son-edge and we consider it as a simple continuation, 
removing the (virtual) vertex. Analogously, if a v = —1, the inward father-edge is joined with both 
edges of the unmarked son-edge and with the inward edge of the marked son-edge. The outward 
father-edge is joined with the outward edge of the marked son-edge (and we consider it as a simple 
continuation). In Fig. |4]we illustrate this procedure with an example: given the graph G G Qz t i in 
a), we first double each edge of G and equip the new edges with opposite orientation in b). Then 
we define the new vertices: for example, if (a v , o~ w , a z ) = (1, 1, —1) we obtain the Feynman graph in 
drawn in c). Hence, for G G Q n ^k an d o G {±l} n , T(G,o~) is a Feynman graph with n + k incoming 
and n + k outgoing leaves (according to the arrows) and with k outgoing and k incoming root edges. 
The vertices of T, V(T), have four edges: two incoming and two outgoing. The set V(T) is in a 
natural one-to-one correspondence with V(G). 

We will now show that every element T G J- n ,k can be obtained as T(G, a) for some G G Q n ,k and 
a G {±l} n . This representation is, however, not unique. The ambiguity comes from the fact that 
although the vertices of each component of G are partially ordered (according to their distance to the 
root), similarly for the components of T, but the order v -< v' in the classical graph for v,v' G V(G) 
is lost in the graph T if v and v' are assigned to different components of V. Therefore two different 
G or a may lead to identical Feynman graphs (see Fig. [5j the Feynman graphs a) and b) are the 
same element of .7-3,1, but they are obtained from two different classical graphs in c) and d)). 
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a) Classical graph 



v 



w 



z 



b) Doubling the edges 



c) Quantum graph 



Figure 4: From Q rhk to T n ^ k 



b) 



c) 



• •- 



d) 



• • 



Figure 5: The graphs a) and b) are the same element of .7-3,1, but c) and d) are different elements of 
£3,1 
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We can equip each Feynman graph with an extra ordering to make this construction unique. 
First we introduce the notion of orderable pairs of vertices. If v, v G V(T) and they belong to the 
same tree-pair but not to the same tree, i.e., v G Tj and v G Tj for some j, then the pair (v,v) is 
called orderable if there is a leaf-pair (£, £') such that v is on the route from £ to the root of Tj and 
v is on the route from £' to the root of Tj. 

For any T G T n p. we define the set of orderings 0(T) on the vertices V(T) as follows. An ordering 
-<; belongs to 0(T) if: 

1. whenever v,v G l^(r) belong to the same tree, then v -< u if and only if u lies on the route 
from v to the root of the tree. 

2. If v € V(r) belong to different tree-pairs, then there is no order relation between them. 

3. Any orderable pair (v,v) is ordered by -<. 

It is easy to see that 0(T) is not empty; for example the ordering v -< v for every orderable pair 
v £Tj,v £ Tj within the same tree-pair is compatible with the requirements. 

A Feynman graph T G T n k equipped with an ordering -<£ 0(T) is called ordered Feynman graph. 

It is easy to see that every ordered Feynman graph T G J 7 " k can be uniquely represented as 
r(G, a) for some G G Q n ,k an d cr G {±l} n . Moreover, since for any (unordered) Feynman graph 
r G J-nk the set of orderings 0(T) is not empty, we obtain in particular, that every Feynman graph 
r G T n ^ can be represented as T = T(G,a). 

The number of Feynman graphs therefore satisfy the bound 

\Fn,k\ <2 n \G n ,k\ <2 M . (9.7) 
9.1.3 Amplitudes of Feynman graphs 

Each Feynman graph T represents a map acting on density matrices; it encodes how the initial 
density matrix changes as the system undergoes a specific sequence of collisions. In this section we 
describe the kernel of this map, commonly known as the amplitude of the Feynman graph. 

Given an arbitrary quantity x e defined for e G E(T), and a vertex v G V(r), we will use the 
notation Yleev ^ x e to indicate that x e is summed with a plus sign if the edge e is incoming (w.r.t 
v) while it is summed with a minus sign if e is outgoing. 

For any T G T n ^, we choose a family r\ = {rj e } e eE(r)j with the property r\ e > for all e G E(T), 
and such that, at every vertex v G V(r), 

J2 ±T eVe = . (9.8) 

(Recall that r e = 1 for outward and r e = — 1 for inward edges.) It is easy to check that (|9.8p is 
equivalent to the requirement that the rj associated with any father-edge equals the sum of the rj 
associated to its son-edges. In particular, the values of r\ on each of the 2(n + k) leaves uniquely 
determine rj e for all edges e G E(T). For a given T G J- n ,k, and a given family 77, we define the 
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operator Kr t t,r) through its kernel 



*r A „(q*,qii *n + k, r' n+k ) := j— — £ \[ (-ir^'K^ %«; (e) - r«; (e) ) 

n da e d P e n ^-i e 1(e )) n s^-^e^^^^-^^ 

edE 2 {T) e<=R 2 (r) e&L 2 {T) 



n ^zi^r n ffE^WE*). 



(9.9) 



where = q if e points away from the root and g" e = q' if e points toward the root and similar 
notation is used for the r variables. Here the map tti labelling the roots is fixed, since T is considered 
as a graph with labelled roots. The map iT2 : L(T) — ► {1, . . . , n + A;} labels the leaves of V in such a 
way that the two elements of a leaf pair receives the same label. Since there is no natural order of 
the leaves, we sum over all possible labelling TT2, and we divide the result by the number of labelling 
{n + k)\. Notice that Krj,r] maps operators on L^(]R 3 ( n+fc )) into operators on L^(R 3fc ) by the formula 

(^r,t,r,7 (n+fe) )(qfc;q'fc) = J ^r,t,r,(qfc,qfe;r n+fe ,r; +fe ) 7 ( n+fc )(r n+fc ;r; +fc )dr n+fc dr; +fc 

where 7 ( n+fc ) is an operator on Ll(R 3( - n+k )) given by the kernel 7 < -™ +fc ^(r ri+ jt; r^ +fc ) in Fourier space. 

Note that there are I-R2I + I -^2 1 + \V\ = 4/c + 3ra — 2ki momentum delta-functions involving p e 
variables and only \E2\ = 2k + 2>n — k\ momentum integration variables. Together with the k\ 
direct delta- functions related to the roots in Ri(T), we see that the kernel K^t-n contains 2k delta- 
functions among its 2n + 4k variables. This corresponds to the 2k momentum conservation in each 
of the 2k components of T. It can easily be seen that all the p e momenta can indeed be uniquely 
expressed through the external momenta r n+ fc, r' n+k , q^, q^. In particular, the dp e integrations are all 
well defined and they correspond to substituting the appropriate linear combinations of the external 
momenta into p e . 



We will represent the fully expanded terms in the Duhamel expansion (|9.6p as a sum of contribu- 
tions associated with the Feynman graphs. More precisely, we will show in the next subsection that 
the m-th term in the sum on the r.h.s. of (|9.6p can be rewritten as the sum of Kr t t,r]"fo k+m ^ over all 
T G T m ^ independent of the choice of r\. 

On the intuitive level, it is possible to recognize the origin of some of the factors in the formula 
(|9.9p for the kernel K-p^,-q- The appearance of the one-dimensional variables a e and the propagators 
(a e — Pe + ir e r/ e ) _1 , for example, derives from the free evolution e~ ttTeP ^, expressed as 

2 r e —itT e (a e +ir e r) e ) 



e -Ur ePe = ( - Te) I ^ . _ (g l0) 

JR a e -pi + lT e T] e 

The presence of the factor S(Y^ e&v ipe) (conservation of momentum at every vertex), on the other 
hand, is due to the translation invariance of the interaction, and it just reproduces the kernel of 
the operators in momentum space (see (|9.5p ). Finally, one can understand the presence of 
the factors S(Y^ eGv ±a e ) as the result of the integration over the time- variables, after all the free 
evolutions e ±ltp = have been rewritten in terms of resolvents according to (I9,10p , 
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\v 



u 



Figure 6: Integration over a e 



The absolute convergence of the da e integrals in (|9,9p can be proven by induction on n. For 
n = there is no such integration. For n = 1, the da e integrations are of the form 

da^da^da^da^, (5(a e2 ~\~ cx u a ei a^,) ^ 
(«ei - Pi x + ir ei r] ei ){a e2 - p\ 2 + iT e2 r] e2 )(a u -pl + iT u rj u )(a w - pi + ir w r] w ) 

with r ei = Te 2 and r u + t w = 0. Here e 2 corresponds to the marked son-edge of the father-edge e\ 
in T, and u, w correspond to the doubling of the unmarked son-edge (see Fig. [6]). Recall that, by 
definition r/ ei = r/ e2 + r} w + r/ n , and all the r/'s are strictly positive. 
We will use the following inequality with r\ = min(r/i, 772), 



da 9 f da 

< 77 



\a — a — irji\\a — b + ir\ 2 \ Ju |(a — a)/?? — i||(a - 6)/?? + z 



dy Crj 



(9.12) 



i 1 1 S |y — ^| 1_ f \a — b — irj\ 



1-e 



for any a, b £ R, r/i, 772 > and < e < 1 (in the last inequality we applied Lemma flO.ip . 
Using this inequality (|9.11|) can be bounded by (with 77 := m.m(rj e2 ,7] w ,rj u )) 



Crj- 



da ei da e 



Ki -Pei +H We 2 ~Vl 2 +IV\ l«e 2 ~ «ei + ~ Pi + £ 
< ^2 -2e /" dOea 

A K 2 ~Pe 2 + i7 ?l K 2 ~v\ +Pl -Pi + ir]\ 1 ~ 2e { ' ' 

C 3 r]- 3£ C 



-3e 

\p1 2 ~ pIi +pI~Pw + H 1 '^ ~ v 



where we applied a slightly modified version of (19.12h two more times, and where we assumed that 
< e < 1/3. 

For a general n we note that any Feynman graph can be built up from the trivial graph with 2k 
roots by successively adjoining new vertices to leaves. When we adjoin a new vertex to a Feynman 
graph T G ^* n _i fc, we select a leaf, e 6 L(F), split it into two edges, e\ and e = e 2 (e± becomes the 
new leaf and it is equipped with the same orientation as e) and we adjoin two more edges, u and w 
of opposite orientation (see Fig. [7]). 

We create three new denominators, three new a variables and one new delta function among 
them. The additional integration is 

' <U '' ! <1 "" <1 -S(a e + a u -a ei — aj) (9.14) 



(a ei - p\ + iT ei r] ei )(a u -pl + iT u r] u )(a w - pi + ir w r] w ) 
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Figure 7: Insertion of a new edge 



where the 77's are chosen such that r? e2 = rj ei + r) w + r\ u . This integral is absolutely convergent 
uniformly in a e and for any choice of rj ei ,rj u ,r] w > and for any choice of the p- variables. This 
guarantees that the absolute convergence of all the da e integrals in (|9.9p can be checked with a step 
by step reduction. 

Now we prove that the right hand side of (|9.9p is independent of the family rj = {i]e}eeE(r) • Notice 
that, because of the condition that at each vertex the r\ associated with the father-edge equals the 
sum of the 77's associated with the son-edges, the only independent ry's are the ones associated with 
the leaves of T. It is moreover clear that, for every fixed e G L(T), Kr t t,r), as a function of rj s , has 
an analytic extension in the whole half plane Re rj s > 0. It is therefore sufficient to show that Kr,t,r) 
is constant in a small neighborhood of a given value rjg with Re rj s > (while all the other rj e > 0, 
e/e£ L(T), are kept constant). After replacing r/g by ?7e+£, we can shift the a e variables as follows. 
For e G E(T) on the route from e to the unique root connected to e, we shift a e — > a e = a e + ir e £, 
while for all other e G E(T) we leave a e = a e . Here we assume that |Re£| < min e£E ^ Re r/ e to 
avoid deforming the a e integral contour through the pole at a e = pi — ir e r] e . Then we can see that 
the integral remains unchanged. This follows from observing that for all v 

±a e = ±a e 

thanks to the definition of r e and the sign convention of the summation. This proves the independence 
of (|9.9p from the family r\. In particular we can use the simpler notation i^r,t- 

We note, that one may introduce dual variables a e for the trivial components e G Ri(T) as well. 
Using the identity (I9,10p . one then may rewrite the definition of ifr.f; into a more compact form 

Kr,t (q* , q'fe ; r n+k , r' n+k ) = ^ + ^ 

x / / n daed ^ n s (pe-Qt { e)) n ^-rt^ (9.15) 



'eG-E(r) eefl(r) eGL(r) 

Xe-^-CD ^) JJ 1 n <E^)<E^ 

However, this definition will not result in absolute convergent integrals, therefore (|9.10p has to be 
used for each e G -Ri(r) before any estimates. With this remark in mind, we can use the simpler 
formula (|9.15p . 

The amplitudes i^r,t will describe the terms of the summation in (|9.6p . The input of the last 
term in (|9.6|) is somewhat different from the previous terms; it involves the density matrix ^ k+n ^ at 
an intermediate time s n and the last free evolution is absent. We thus introduce a slight modification 
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(9.16) 



of the amplitude ifr,t to represent this term. We define the operator Lr,t, for T G .F^fe, n,k > 1, 
through 

£r,f(qfc,qfc;r n+fc ,i4 +fc ) : = 7—-^ XI X a * [ [ II doLe II d ^ e 

X II ^Oe-^^^Oe'-g^^)) II S ( Pe - r ^2(e))S{Pe' ~ < 2 ( e) ) JJ *(j^±P< 

ee-R(r) egL(r) i.ev(r) eei- 

xexp(-ii X r e (ae + ir e »/e)) II ~ - J +iTr , II 5 (X ±a < 

egR(r) =6B(r) e e ' e ^„ G y(r) e€v 

where M(T) C V(r) is the set of maximal vertices of T, that is the set of all v G V(r) so that there 
exists no v with v >- v (recall that -< was the partial order on V(r) induced by the distance to the 
root). Moreover, for a vertex v G V(T), we denote by S v the set of son-edges of the vertex v (clearly, 
R(T) n S v = 0, for any v G V(r)), and we set 

0y = y~] T e , 

eeS v 

i.e., a v = 1 if of the four edges adjacent to v, three are outward and one is inward, a v = — 1 otherwise. 

In other words, Lyj is defined so that there are no propagators associated with the son-edges of 
a v G M(r). The variables r] = {r]e}eeE(r) are chosen as in (19.80 : although r/ e , for e£5 5 , does not 
appear directly in (|9.16j) . the value of the rj- variable associated to the father-edge of v depends on 
it. Analogously to Kr t t, it can be proven that Lr,t is independent of the choice of rj. In (|9.16p we 
define Lrt using a form analogous to the simpler definition (|9. 15|) : clearly Lrj can be rewritten into 
a form analogous to (19. 9p as well after integrating out all a e , e G i?i(r). 



9.2 Duhamel expansion in terms of Feynman graphs 

In this section, we shall prove that the Duhamel expansion for the solution (|9.6p can be expressed as 
sums over Feynman graphs. Similar representations were used in the physics literature and proofs 
are available for the expansions of the imaginary time many-fermion Green functions via the method 
of (Grassmannian) functional integrals, see e.g., [21]. Our relatively short proof is elementary and 
does not rely on functional integral. It also allows for an explicit remainder term. 

Theorem 9.2. Fix k,n > 1. Then, for any Jq that is symmetric with respect to permutations 
( in the sense of II 1.8\) ), we have 

dsj ds 2 ... ds n U { Q k \t - 81 )BW . . .^ fc+n - 1) (,„_ 1 - s n )B^-M k+n \sn)l$ +n) 

Jo Jo Jo 

K T,t% 

(9.17) 

For n = 0, k > 1 we have 

"o^hf = E K ^o ] (9-18) 
where the summation is only for the trivial graph. 
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Pa« a 



P b a t 



Pd «d 



Pe a e 



Figure 8: The two Feynman graphs in T\\ 

Moreover, for any fixed k,n> 1, if ^ k+n ^ g C([0,T];Tik) is symmetric with respect to permuta- 
tions for all t E [0, T], then we have 

ft fS\ fS n —l 

dsj ds 2 ... ds n 4 k \t - Sl )B^ . ..U^ n - l \s n ^ - s n )B^H k n +n) 

Jo Jo Jo 

/•t 9.19 

£ / dsL r , t -d k+n) 

--7T . JO 



-I 



for all t £ [0,T]. 



Before presenting the proof for the general case, it is instructive to show how the structure of 
the operators K-pj emerges from the Duhamel expansion. We thus consider ()9.17|) in the very simple 
case n = 1, k = 1, i.e., 

f dsU^\t - s)B^U {2 \s) 1 {2) = Kr u a® + K T2 ^f\ (9.20) 
Jo 

where Y\ and T 2 are the two elements of \ drawn in Figure El 
By definition of the map B^\ the l.h.s. of (|9.20p is given by 



-i / dsUl ) 1) (t-s)Tr 2 5(x 1 -x 2 )U^\s)^ 2) +i / ds ll£\t - s)Tv 2 U^ 2) (sh (2) 5(x 1 - x 2 ) . (9.21) 
Jo Jo 

We now show that the first term coincides with the contribution of the Feynman diagram Ti (anal- 
ogously one can prove that the second term equals the contribution of T 2 ). 

Recall that Ky 1 1 £ T\ \ maps operators on L 2 (IR 3 x M 3 ) into operators on L^(R 3 ). In momentum 

(2) (2) 

space, the kernel of 7q has four variables, 7q yPoPd'iPaiPe)- These are the momentum variables on 
the right hand side of Feynman graph T\. The resulting operator, Ky^I^, acts on L 2 (M 3 ), and its 
kernel has two variables, (pb',Pa) represented on the l.h.s. of the Feynman graph. 

(2) 

The pairs of variables, (p c ]Pa) and (pd',Pe), m the argument of 7q correspond to the input and 
output of the first and the second variable of the L 2 (IR 3 x R 3 ) space. Similarly, p^ is the input and p a 



(2) 

is the output variable of the resulting operator KyjTq ■ arrows in the Feynman graph pointing away 
from the roots indicate input- variables, arrows pointing towards the roots indicate output-variables 
of the corresponding density matrix. 

For any kernel ^ 2 \pi,p 2 ',p'i,p 2 ), the free time evolution in momentum space acts as follows 



[U {2 \s)^){p 1 ,p 2 -p^p' 2 ) = e --^+fi) e -(W] 2 +^] 2 )7j 2 )(p 1 ,p 2 ;p' 1 ,^) . 
The multiplication by 5(xi — x 2 ) corresponds to convolution with 5{p\ + p 2 ) in Fourier space, 
(6(xi - x 2 )'y {2) ^j{p 1 ,p 2 ,p 1 ,p 2 ) = / dr 7 (2) (r,pi - r + p 2 ;p' 1 ,p 2 ) , 
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thus after taking the partial trace, we have 

Tr 2 (£(xi - x 2 )7 (2) ) (PWP\)= J drdq 7 (2) (r,pi - r + q;p[,q) 
Applying these elementary steps for the first term of (|9.2ip . we get 

-i / dsU^\t-s)TT 2 5(x 1 -x 2 )l{^(s)4 2) (p b ;p a ) 
Jo J 

= ~ij dae-iQ-'Wjtt-'M J dqdr (lA {2) ' (s)^) (r,p b - r + q;p a ,q) 



: J* dse-^-^le^ti J dp c d Pd dp e 5(p b +p e 



Pc-Pd) 0»)7o ) (PoPdlPa,Pe) ■ 



(9.22) 



In the last step we changed variables, which now correspond to the variables in the Feynman graph. 
In particular, the vertex involves a delta function expressing the Kirchoff law (conservation of mo- 
mentum at the vertex). 

Now we show how the time integrals of the propagators can be expressed in terms of auxiliary 
a-integrals and resolvents. Neglecting the momentum integrations, the last term in (|9.22p contains 
the following propagators: 



n 



d S e ~ i ( t ~ S )Pb e *(*- S )Pa e - is {Pc+P d ) e is (Pa+Pl) _ 



Using the identity (|9,10p for the propagator of each momentum variable, we obtain 



U = i / ds 



da a dabda c dadda e e it ( a °.- i Va)-i(t-s)(a b +ir] b )-is(a c +ir] c )-is(a d +ir] d )+is(a e -ir] e ) 

(«a - Pi ~ ir] a ){oL h -p 2 b + iVb)(a c ~P 2 C + ir}c) (ad -P d + M7d) i a e ~p\~ «7e) 

(In our example, r b = r c = = 1 and T a = r e = — 1 by the definition of r). Notice that the time 
integration can be extended to s £ (— oo,oo), since all 77's are positive and by residue calculation 



da e - is( - a+ir >) 
a — p 2 + irj 







if s < and 77 > 0. 

Using that r\ b = rj c + r]d + rj e and performing the ds integration, we obtain the delta function in 
the a variables: 



n = i 



e it{ aa -i na ) e -it{ ab +i nb ) da a da b da c da d da e 5{a h - a c - a d + a e ) 
(a„ -pl~ ir) a )(ab ~ p\ + ir] b )(a c -p 2 + ir]c)(ad-p d + iVd)(a e -p\- iVe) 



(recall that the 5-function in the a-variables is defined w.r.t. to the measure da = di^,®/^)- 
Combining this formula with (|9.22[) . we arrive at 



/ / dsU^(t - s)Tr 2 8(x! - x 2 )^ 2) ( S ) 7 J 2) ] (p b ; Pa ) 

= j dp c dp d dp e 5(p b +p e - p c - PdWsP (Pc,Pd',Pa,Pe) 



(9.23) 



n 



da, 



j=a,b,c,d,e 



ay - p i + irtfj 



5(a b + a e - a c — ad) 
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which is exactly the action of the kernel K? 1 1 (as defined in (|9,15p ) on 7q ' . 
Now we come to the proof of the general theorem. 

Proof of Theorem\2M We start by proving (EOT)) and (I^TB|) . For k > 1, n > 1, and t G [0,T], 

let 

:= fds, / Sl d S2 ... ds n ut\t-s 1 )B( k K..4 k+n - 1 \s n ^-s n )B^^ 

Jo Jo Jo (9.24) 

For k > 1, n = and t G [0, T], let 

^ufHthP • (9-25) 
We also define / for > 1, n > 0, and i G [0, T] through its kernel given by 

<!*) : = X] y' dr ™+fc dr n+fc^r,i(qfc,q^ (9.26) 



We need to show that / = 6 n [ for all n > 0, A; > 1, and for all t G [0, TJ. The proof is based on 
induction over n. For n = the claim is trivial using the identity (|9.10p . thus from now on we assume 

(k) (k) 

that n > 1. We prove that ^ t = 6 m t for m = n and for all > 1 and f G [0, T], assuming that this 

(fc) (k) 

is true for all m < n. We will first check that Q q = 0^ q, then compare their time derivatives. 

At * = clearly £w = 0. To see 9 n q — 0, we check that the -fCr^ kernel vanishes at t — 0, for 
r G J-" n ,ki n > 1. We use the representation (|9.9p . Since n > 1, we know that ^(r) ^ 0, in particular 
there is at least one propagator factor (ag — p\ + iTgr/g) -1 on the right hand side for some e G ^(T). 
Notice that for i = the exponentially increasing factor n e e_R 2 (r) e??ei * s no ^ P resen f • Using the 
absolute convergence of all the vertex integrals (19. 14ft uniformly for large r] and the estimate f)9.13|) 
for the integration of a maximal vertex, we can let r/ e — ► oo for all e G ^(r) (recall that the 77's on 
the leaves determine all rj values), and conclude that ifr,t=o(qfcj q^S r n+ fc, r^ +fc ), for any fixed values 
of its argument, is bounded by a quantity converging to zero as rj e — > 00, for all e G L2(r). Since 
i^r,t=o is independent of 77, we obtain that it is zero. 

Next we prove that and 9^ t both satisfy the same equations: 

k 

idt&Mk) = E («? - (^) 2 ) €g(q*;qi) + < (S«^53) (q*;qi) 

^«(q*;qi) = Efe 2 - (^) 2 ) ^](q*;qi) + • (s«e2+3) fa^)- 
i=i 

Since, by induction assumption, £^*q = and ^^L*^ = Gn-i \ f° r everv & > 1 and every i G [0,T], 
it follows from (19371) that = 0^ for all k > 1 and i G [0, T]. It remains to prove ([9371) . 

We start by deriving the equation for t . To this end we compute the derivative of / with 
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respect to t. 

'■t 



(9.28) 



Jo Jo 

+ iV dj ds 2 ... " ds n \Aj, U { k \t- Sl )B^ ... 

• =1 JO JO JO L 

Hence, in Fourier space, 

k 

«fc$(q*; <*) = E («? - (^) 2 ) + < (^d-^) (q* ; <&) , (9.29) 



which proves the first equation in (|9.27|) . Next we show the second equation in (|9.27|) . 
Using (|9.9p , we compute the time derivative of 6 n t : 

1 E n (--e)e-^ ( ^« )2 %i; (e) -ri; (e) ) 



re.F„ 

x 



(n + k)\ 

n da e d Pe n 5(p e -gi e 1(e) ) n ^-^ (e) ) 

ee£ 2 (r) ee/? 2 (r) eeL 2 (r) 

x e -^ eeR2(T) Ma e+ ir eVe ) I £ r e (i; (e) ) 2 + £ r e (a e +^ e ) 

\eei?i(r) eeR 2 (r) 

- n ^—^177 n <ewke^)- 

(9.30) 

(Strictly speaking, this calculation is formal, since after the differentiation the da e integral is not 
absolutely convergent. We will remark on this issue at the end.) We write 

E Te(a e + iTeVe) = E r e («e ~ £>e + ZT e 7? e ) + E TeP e ' (9.31) 

eei? 2 (r) eeife(r) ee-R 2 (r) 

Because of the delta- functions rieei?, 2 (r) ^(Pe ~ sf^e))' we a ^ so nave 

E ^ l(e )) 2 + E TePl= E ^*(e)) 2 = i;W-W) 2 )- 

ee.Ri(r) eei? 2 (r) eeiJ(r) i=i 
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From the last two equations we obtain 

k 



re^ n , fe i=i 7 
+ dr n+kdr' n+k ^ n+k \r n+k ;r' n+k ) 

"ot^ n <-*.>^<^ <<.,-*,.>> 

v tt 2 eg_R!(r) 



n da e d Pe n 5(p e -4 (e) ) n ^ 

e£E 2 (r) eei? 2 (r) egL 2 (r) 

x e -i*Ee e « 2 Cr)^(ae+ir eJ7e ) | £ r e (a e - j£ + iT e 7fe 

\eeR 2 (r) 

- n ae -J +i ^, n *(eH<e*. 



e 7T2(e)^ 



(9.32) 



For a given e G -R2(r) let v = v(e) G V(r) be the only vertex such that e S v (by definition of 
i?2(r) there is such vertex). Then the second term on the r.h.s. of (|9.32p can be rewritten as (we 
use that L 2 (T) n R 2 (T) = 0): 



E Tg / dr n+kdr' n+k ^ n+k) (r n+k ;r' n+k ) 
„ 1 \^ T T -a -i^(^r R 2 



E II (— e)e^^<^ %J l(e) -r[; (e) ) 



(n + fc)! 

V y 7T2 eG-Ri(r) 



/" / nd« e dp e e -*M«s+fc*») <5 (E±a e ) 5 ( ) 

^ ^ K eGi) \ee-0 / \eGS / 

n d« e dp e n 6 (pe-qi 1(e) ) n ^-^v 



(9.33) 



=6i5 2 _(r) eGfl 2 (r) eGL 2 (r) 

e^z; e^e 

x g-**E e6 H 2 (r),e^g r «»(a!e+jr e rj e ) 

- n n *(£*«.)<£**)■ 

e^e U^i; 

Next we perform the a§ integration. Using that, by definition of the family rj, rjg — X/e^egs ^ 
and the fact that 

E ±«e = <^=>- r g ag = E TeQ;e (9.34) 
from the definition of r e , we obtain 



(9.35) 
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Roots 



Roots 



Figure 9: The map (r, e) — » T 



Multiplying this contribution with the factor exp(— it ]Cee-R2(r) e^e T e( a e + ii~ e r} e )) from (|9.33p and 
using (|9.34p . we obtain exp(—it^2 e€R T e (a e + iT e r] e )) with R = {e E -R2CO :e/e}U{e£ij:e/e}. 
The set RU R± can be interpreted as the set of roots of a new graph, with k + 1 root-pairs and with 
n — 1 vertices as follows. 

Given T € and e E ^(T), we define a new graph T = T(r,e) E T n -\^+\ as follows: the 
roots of T are given by -R(r) = {e E -R(r) : e / e) U {e £ ii : e / e), that is we remove the edge 
e, and we add all other edges adjacent to the vertex v to the set of roots (recall that v E V(F) is 
the unique vertex to which e is adjacent). The label ii\ of the roots in T is defined so that the two 
unmarked son-edges of e at the vertex v become the (k + l)-th pair of roots of V, while the marked 
son-edge at v inherits the label of e, and all the other roots keep their label. The vertices and the 
edges of the new graph T are ^(T) = V(T) \ {v} and, respectively, E(T) = E(T) \ {e}. Finally, the 
leaves of T are the same as the leaves of T, that is L(T) = L(T). Graphically the map from (r, e) 
(with T E T n ^ and e E i?2(r)) to T E ^"n-i.fe+i corresponds to the cancellation of the edge e and of 
the vertex v and to moving the two new roots (the two unmarked son-edges at v) with their full tree 
graphs to the bottom of the graph (see Fig. [9]). 

Note that the map from T E T n ,k and e E i^CH to T E Fn-i^+i is surjective but not injective: 
in fact for every T E ^n-i.fe+i there are 2k possible choices of T E and e E i?2(r), because the 
last (k + l)-th pair of roots of T can be attached to any of the first 2k roots. This implies that the 
sum over T E T n ^ and e E R2^) in (|9-33[) can be replaced by a sum over T E Fn-i,k+i and a sum 
over the first k pair of roots of T plus a binary choice between the two paralell root edges. 

With all these notations, (j9.33[) can be rewritten as 

/k n k 
dr n+k dr' n+k ^ n+k) (r n+k ; r' n+k ) ^ ^ J dq fc+1 dq' fe+1 5(q e - qe)S(q' e - q t ) 



t+3 



x e 



<%j - Qj ~ Qk+l + tfk+lWj ~ %) ~ %j - % ~ Q'k+l + 9fc+l)<% - qj) 

n da e d Pe n d(pe-ti(e)) n ^-^p 

--' " " m e eE(T) e<E-R(f) egL(r) 



— F 



(ra + k)\ 



(9.36) 



eeE(r) 6 veV(T) eev 
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In the second line we rewrote the integration over the momenta p e , e £ v, in (|9.33p . so that it is 
clear that it describes exactly the action of the operator iB^ (see (|9.5p ). From the last equation, 
together with (|9.32p . we obtain 

id t e$M k ) = £ (<zf - (q' f) ^(q fe ;q' fc ) +z (q fe ;q' fc ) (9.37) 

which proves the second equation in (|9.27|) . This completes the proof of (|9.17|) and (|9,18p . 
Let us now prove (|9.19p . The l.h.s. of (|9.19p can be rewritten as 

d Sl / d* a . . . / d Sn 4 k \t - S1 )B^ . . . ui k+n - X \s n ^ - Sn )B( k +^ 7 t- n} 



o Jo Jo 

t rt—s /-Si 

ds 



o 



ds 1 r\..r- 2 ds n „M k \t-s- sl )BW... O^) 

Jo Jo 

x U { t n - 2 \s n . 2 - Sn^B^-Vut^Hsn-i)} B^H k+n) ■ 
From (|9.17p . the last expression equals 



f ds V K f B(*+»-l) 7 (*+») 
Jo ~ _ ' 



J t k) :-. I d.s «'r,_^"'' + "" IU ; .i A '"" 1 - 1 9.39 i 



Using (|9.15p and (|9.5p . we obtain 

„ t n+fc-1 

4*W) = -«7 ^ E (w+fc _i), E E 



y dr n+fe _idr; +fe _ 1 y J JJ da e dp e 5(p e - q^) JJ 5(p e - r^) 



'eGE(r) ee-R(r) eSL(r) 

"" ,S -* N n ac J M n '(e^we^.) < m °> 

x / d Pn+fe d P ; +fc ]JS(rt - pe)5(r' e - p' e ) ^ n+k \p n+k ; p' n+k ) 

6 ( r j - Pj) 6 ( r j - (Pj + Pn+fc - Pn+k)) ~ <% - Pj) S ( r 'j ~ (Pj + Pn+fc ~ Pn+k)) 

where the sum over j and the last two lines correspond to the action of the operator B^ n+k ^ 1 ^ on 
the density ^1 ■ Note that j actually labels the leaf-pairs of T. Fixing j = 1, . . . , k + n — 1 and 
one of the two terms in the square bracket on the last line of (|9.40p is equivalent to choosing one of 
the leaves of V. For a given T and e £ L(T), we can define a new graph T £ J^ n ,fc by splitting the 
edge e with a new vertex and attaching two new leaf edges to this vertex (see Fig. [T0|) . The same 
graph in T n ,k can clearly be obtained starting from different graphs T £ ^Vi-i,jfc- More precisely, 
if we denote by M(r) the set of maximal vertices of T (u £ M(T) if v £ V(r) and there is no 
v £ V(r) with v ~< v), then we find all possible T by removing a vertex u £ M(T) (and deleting the 
son-edges of u). It is hence clear that, in (I9.40p we can replace the sum over T £ J- n -l,k, the sum 
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r 



r 



r 



Figure 10: The map (T, e) — > T 

over j £ {1, . . . , n + k — 1} and the binary choice of one of the two terms in the last line, by a sum 
over all T £ J~ n .k an d over all v £ M(T). 

In order to rewrite (|9,40p in terms of the new Feynman graph V £ J- n ,k, we observe from (|9,40p . 
that the son-edges of v will have a momentum like all the other edges, and that momentum conser- 
vation holds at v (see the last line of (|9.40p ). but they will not have any a- variable, any r\- variable, 
and any propagator. The labelling iT2 of the leaves of T induces a labelling with {1, . . . , n + k — 1} 
of the leaves of T, with the exception of the two unmarked son-edges of the chosen v £ M(T): these 
two edges are always labelled by the number n + k. Of course, because of the permutation symmetry 
of ^ n+k \ we can restore a full symmetry of the leaf- variables; to this end we replace the sum over 
7T2 £ S , n+ fc_i by a sum over 1T2 £ S n+ k and we replace the factor (n + k — 1)! by (n + k)\. We conclude 
that 



x 



/ n ^ n ^ n n ^-^w 



eeE(r) eeB(r)\s« eeH(r) 

>< n 

e£EE(r)\S s 



eeL(r) 



(9.41) 



where S$ is the set of son-edges of v and we recall that = Yle&Sv Te ' From Q9.16P we obtain (|9.19p . 



Finally we comment on how to make the formal differentiation in (|9.30p rigorous. In the defi- 
nition (|9,9p we could have defined i^r.t^e instead of K^^^ with introducing a regularizing factor 
exp(— ^YleeE2(r) Q e) m the integrals and similarly we would have defined 0^ t£ in (|9.26p . Then the 

(k) 

time derivative of 6 n j. can be computed by differentiating the integrand. Note that before differen- 
tiation, the da e integrals in (|9.9|) are absolutely convergent and the convergence is uniform in e > 
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and t £ [0, T] for any fixed T. Therefore we obtain 



lim 6 

£^0+0 



(A:) 



}(*) 



n,t,e n,t 



(9.42) 



for any fixed n, k and uniformly for t £ [0, T}. In particular, the relation t^ k \ = lim e ^o+o &n0 e = ® 
for n > 1 still holds. The identity (|9.10j) will not hold any more but 



lim 

£->0+0 



00 da P e~ e < e - itTeas 



-oo a e -pi + ITgTj 

•0 _ i;*»> . ._«(*>) 



(-ir e )e- i ' TeP 'e-' ?i 



(9.43) 



still holds, therefore £g t = lim £ _»o+o $o / e - In the integral (|9.35j) we pick up a factor 

\2 



exp 



which converges to 1 in the limit. The integrations in (|9,36p are again absolutely convergent with or 
without the regularizing factors. We therefore conclude the following version of (19. 37ft for e > 0: 



•ft*£k(q*; q' fc ) = £ fe 2 - (^) 2 ) <i(<^ <£) + i { Bik)0 n-il) (q*; q£) + °« (9-44) 

as e — > + 0. Integrating back this system of differential equations, comparing the result with the 
solution to (|9.29f) and using that the difference in the initial conditions vanish as e —* 0, we obtain 
that 

o { n % = + 0(1) 

for any fixed n,k and uniformly on t £ [0,T] for any fixed T > 0. Combining it with (|9.42p . we 
obtain the rigorous proof that the differentiation in (|9.30p is allowed for our purposes. □ 



9.3 Bounds for Amplitudes of Feynman Graphs 

For brevity, we introduce the notation 

for an operator £ IC^ with kernel «/^(qfe;ql.) expressed in momentum space. We also define 
{j( k \ £r,t7 + ) similarly. In the next two theorems we show how to bound (J^ k \ ifr,t7^ n+fc ^) and 
(«/^,-£r,t7 + *^) (for r £ J- n> k and for an observable decaying sufficiently fast in momentum 
space) in terms of the 7^ n _|_fc-norm of 'y( n +*0. By Theorem 19.21 this will allow us to control the 
Duhamel expansion (19. 6p of any solution of the infinite hierarchy (]9.2p . Recall that the contribution 
(j( fc ),K r ,i7 (n+fc) ) shows up in the analysis of the fully expanded terms in (|9.6p . while (J^ , Lr,i7^ n+fc ^) 
shows up in the error term of (19. 6p . 

Theorem 9.3. Fix k > 1. For any n > 1, suppose ^ n+k ^> is non-negative and symmetric with respect 
to permutations (in the sense of (1.8\) ). Assume < t < 1. Choose £ ICk that is symmetric with 
respect to permutations and whose kernel satisfies 

k 
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Then we have 



for every n > 0. 



(J( fc ), i^7 (n+fc) ) < C n Tr (1 - Ax) ... (1 - A n+fe ) 7 (' 



n+k) 



(9.46) 



Theorem 9.4. Fix k > 1. For any n > 1, suppose 7 ( n+fc ) is non-negative and symmetric with 
respect to permutations (in the sense of il.8\) ). Assume < t < 1. Choose £ /C^ symmetric 
with respect to permutations and with a kernel satisfying \9.4-5 ). Then, for every n > 10 + fc/2, w;e 



(JW , Lr, t7 (n+?£) ) < C n (1 - Ax) ... (1 - A„ +fe ) 7 ("+ fc ) . 



(9.47) 



Remark. The integer A; > 1 is fixed. The constant C in (I9.45h depends on k, and so do the 
constants on the r.h.s. of (|9.46p and (|9.47p . The restriction t < 1 plays no significant role in the 
theorems; it simplifies the proof in a trivial manner. It is also possible to obtain a i-power in (j9.46j) 
but we will not need it. Theorem 19.31 will be applied to control the fully expanded terms in the 
Duhamel series (19.6p . which, by Theorem 19.21 can be expressed in terms of the kernels lfr,t- F° r ° ur 
proof of the uniqueness of the solution to the BBGKY hierarchy, it will be sufficient to show that 
these terms are finite. They involve only the initial condition, therefore these terms will be identical 
for any solution with the same initial data. On the other hand, Theorem 19.41 will be applied to 
control the last term in (|9.6p involving a density matrix at an intermediate time s n . In this case it 
is not enough to show the finiteness of the contributions; we also need to prove that they are small. 
It is for this reason that in Theorem 19.41 we need to extract a time dependence and we will use the 
apriori bound ([93]) and that C n t n / 4 -» as n — ► oo if t is small. Note also that the power re/4 in 
(|9.47p is not optimal and we do not aim at the optimal i-dependence, but we remark that this issue 
is related to exploiting an additional smoothing effect of the free Schrodinger evolution. 



Before proving these theorems, we point out that the estimates (|9.46p . (j9.47j) can be viewed as 
Strichartz type inequalities in the many particle setting. Recall that the Strichartz inequality states 
that 



2/ TO 3\ 



\e^f\\;ds<C\\f\\ r 2 , feL 

jo 

for r,p satisfying f + - = § an d 2 < r < oo. This inequality implies that 



[ \\e isA f\\ P ds < Ct 1 ^ 
Jo 



which means that the free evolution smoothes out possible singularities of / at the expense of reducing 
the t-power. 

Another form of the Strichartz inequality asserts that 



fdse^^fs 
Jo 



jrjV < C\\f t \\ , p > , 



ft = ft(x) 



where L\lF x denotes the space L r (M; L P (1R 3 )) and the positive exponents p,r,p',r' satisfy 

12 1 12 7, 

— I = -, r>2; 1 = -, r'<2. 

p 3r 2 p' 3r' 6 



(9.48) 



Once again, this estimate quantifies the smoothing effect of the free evolution operator. The price 
of reducing the t-power is now expressed in terms of the change from the L r norm to the U norm 
in the t variable. 
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The kernels Krt an d Lv,ti though defined in Green function form, actually have representations in 
terms of n-fold time integrations like the formulae appearing in (|9.17p and (|9.19p (with the operator 
j^{m) defined i n (|8,3p ). If we replace the 5-function (which came from the two-body interaction) in 
the definition of £?( m ) by a smooth function, the correct t dependence in (|9,47p would be t n (at least 
for small t). The estimate (|9.47|) states that the 5 interaction is allowed if we give up some power in 
t — in the same spirit as in the Strichartz inequality. 

Each integration step on the left hand side of (|9.17p and (|9.19p actually involves a time integration 
and a space integration via the partial trace in B^ m \ It would therefore be natural to perform an 
iterative estimate involving subsequent one-particle space-time dispersive bounds. Unfortunately, we 
were unable to find an appropriate one-particle scheme to implement this approach. Our method is 
much more complicated and it involves tracking several singularity structures of the density matrices 
in the integration step. 

One reason to use the Feynman diagram representation is to obtain estimates with correct n 
dependence. From the summation in the definition of i?( m ) (see (I9.4p ). the number of terms on the 
l.h.s of (|9.19p is 2 n k(k + 1) • • • (n + k — 1) ~ n\. This factorial can be exactly compensated by the 
multiple time integration, 

dsi / ds 2 ... / ds n = -, (9.49) 



o 



but only if the L 1 -norm is used in time. Higher L r -norms in time result in a partial loss of the 1/n! 
in (|9.49p and thus the summation of these estimates over n will not converge for any t. 

For this reason, we developed a new method, based on the expansions ()9.17j) and (|9.19p in terms 
of Feynman graphs. Our graphical representation, among its other merits, reduces the number of 
terms in the expansion from n\ to C n (see (|9.7p ). This combinatorial reduction stems from combining 
graphs like b) and c) on Fig. El The dispersive properties of e ltA are now captured by the decay 
properties of the integrands in the kernels Kj^ u ^V,t ( see (|9.15p . ()9.16|) ). These multiple integrations 
can be successively performed and we thus obtain the bounds (|9.46p . (19.471) . Our representation 
treats all smoothing effects simultaneously and thus exploits an additional decay which we were not 
able to obtain with one-particle methods. However, we do not know if it is possible to design a 
one-particle inequality similar to the Strichartz inequality to give a short proof for the estimates 
(HEMP and (|9T471) . 



9.4 Proof of Theorem [9731 

To better explain how the dispersive properties of the free evolution are used in our approach, we 
first discuss the main ideas involved in the proof of (|9.46p on a heuristic level; similar ideas apply to 
the proof of (I9.47p . We have to bound all integrals over the three dimensional momentum variables 
p e and over the one-dimensional variables a e appearing in the definition of the kernel Kyj (see 
(I9.15P ). Notice that, because of the singularity of the (5-potential in position space, here we face a 
large momentum problem: we have to make sure that all integrals over p e (and a e ) are convergent 
in the large p e (respectively, large a e ) regime. To this end we will develop an integration scheme 
dictated by the structure of the Feynman graph T. 

We will start by integrating over the variables p e and a e associated with the leaves of T. The 
momenta on the leaves are exactly the variables of the density ^ n+k \ The factor Tr(l — Ai) . . . (1 — 
A ra+ fc)7^ n+fc ^ on the r.h.s. of (j9.46j) implies that each leaf carries a decaying factor |p e |~( 2+A ) (for 
large p e ), for any A < 1/2. In the rigorous proof we will have stronger restrictions on the value of A. 
The idea then is that we integrate over all the p and a variables, starting from the leaves and then 
moving towards the roots. At each step we propagate the momentum decay from the son-edges to 
the father edge of a certain vertex of T by integrating out the variables of the son-edges. 
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P„ «H 



p w a w 



Figure 11: Integration scheme: a typical vertex 



A typical step in the integration scheme is as follows: choose a vertex v G V(r) such that we 
already have demonstrated a decay |p|~( 2+A ) in the momenta p u , pd, p w of the three son edges of 
v (denoted by u,d and w, see Figure [TT]) . This means that all the momentum- and a- variables of 
the edges which are to the right of u, d and w in the graph have already been integrated out. Then 
we first perform the integration over the three ce-variables of the son edges. By power counting, we 
obtain formally 

5(a r = a u + ad - a w ) const 
da u dadda w - ^- tt—. 5— < 



<«u - Pi) {&d ~ P d ) {<*w ~ P 2 W ) (a r -Pl-Pd+Pl) 

(up to logarithmic factors). Then we integrate over the momenta of the son-edges by using their 
decay factor and we obtain, again by simple power counting, 

dp u dp d dp w 5(p r = p u + p d - p w ) const 



\Pu\ 2+x \Pd\ 2+x \Pw\ 2+x (a r - pi - p\ + pi) \Pr\ 2+x 

Here the power counting requires that 3 (2 + A) + 2 — 6 > 2 + A which holds for any A > 0. Thus 
the same decay in the large momentum regime propagated from the son-edges to the father edge. 
This procedure can then be iterated until we reach the roots of V. At that point we can complete 
the integration scheme by using the smoothness (momentum decay) of the observable . 

In this formal computation, the dispersive nature of the free evolution is expressed via the decay 
in p and a of the resolvent (a — p 2 )^ 1 . The decay in both variables is critical to complete the inte- 
gration scheme. The rigorous proof of Theorems 19.31 and 19.41 is more involved than this simple model 
calculation because the singularity structure of (a — p 2 ) -1 is spherical and it cannot be described 
by a simple power counting alone. We will have to consider edges of various types, characterized by 
different decay properties, and we still have to close the iteration scheme. 

Proof of Theorem 1 9. From the definition of Kr t m (|9.9p . we have 

(/ (fc Ur,t7 M) ) = T^| E / dq fc dql,dr n+fc dr; +fe j( fc )(q fc ;q;) 7 (^ (rn+fcir ; +fc) 

x I] (-ir e ) e -^ ( ^« )2 %i; (e) -r[; (e) ) 

eSfli(T)=Li(r) 



X 



e£E 2 (T) e€R 2 (T) eEL 2 (T) 

n ae -J + ,^ n <£^)<2>«)- 



X 

(9.50) 
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Because of the permutation symmetry of ^ n+k \ the integral has the same value for every choice of 
TT2- Hence, instead of averaging, we fix 1T2 G S n+k . We define the sets 

Qi := {e G L(T) : r e = 1} , and Q 2 = {e G L(T) : r e = -1}; (9.51) 

that is Q\ is the set of outward leaves and Q2 is the set of inward leaves. Clearly L(T) = Q1DQ2 and 
\Qi\ = IQ2I = n + k. We use the notation ^ n+k ^ ({(p e ',Pe')} eeQi) to stress the fact that, because of 
the permutation symmetry, the density -y( n+fc ) only depends on the set of pairs (p e ',Pe') of momenta 
associated with the (paired) leaves of T, and not on the order of the pairs. Integrating over the 
variables q^q^ and r n+k ,r' n+k and using all the delta-functions, the absolute value of (|9.50p can be 
estimated by 

Ce ti: eeR{ r^ [ n d Pe n da e n - — i— — - n j_ 

J ±J f ■ L± ±x \a e — p~ + tT e r] e \ ±J * < n N<3 



e£E(T) e&E 2 (T) eeE 2 (T) '~ c x " e ' ~'""=" eeR(V) ^ Pe/ 

X [] 8(^±ae)5C£±Pe)\l {n+k H{(Pe-,Pe>)}eeQ i : 



(9.52) 



where the factor n e gi?(r)(^e) 3 comes from estimating the observable using the assumption 

(Hasp . 

Since the observable is symmetric w.r.t. permutations, and since K^ y t preserves the symme- 
try, to compute the quantity on the l.h.s. of (|9.50p we can replace the density ^ n+k ^> by its restriction 
onto the subspace L 2 (M.^ n+k ^) consisting of all permutation symmetric functions in L 2 (M 3 ( n+fc )). 
Hence, ^ n+k "> can be written as ^ n+k "> = V. Xj\-(pj)(tpj\, with ipj G L 2 s (M. 3 ( n+k ^) such that \\ipj\\ = 1 
for all j, with Xj > for all j (by the non- negativity of ^( n + fc )) ; and with ^ ■ Aj < 00. Hence it is 
enough to prove the bound (|9,46p for ^ n+k ^> being a one-dimensional projection. In the following we 
therefore assume that ^ n+k) (p n +k] p' n+k ) = ^(p n +k)^(p' n+k )- 

We again use the notation ^({PejeeQi) to indicate that ip is a function of the set of the momenta 
associated with leaves in Qi, and not of their order. Moreover we choose r\ e = 1/t, for all e G L(T): 
this implies that r] e > l/t for every e G E(T), and ^ eg ^(r) ^ e = (2«+l)/i. With a weighted Schwarz 
inequality, we can then bound (19. 52ft by 

c e w n n ,„ "g... n n «(e±««we^«) 

17 ees(r) ee£ 2 (r) 1 e l e ^ * 1 ee-R(r) v ' »ev(r) ee« eg-u 

X f n^l^KW)! 2 + ^|l^({Pe}e G Q 2 )| 2 ) 
\IIeeQ 2 ^e YYedQ x Ve J 

< C n J dp n+k p\ .. .p 2 n+k \^{p n+k )\ 2 

- {, «* / n ^ n n <^ n i n *e±o.)«s>»> 
+ sup / n ^ e n T^zVy n tAs ni n 

= : A + B 

(9.53) 

where we used that, since we assumed that t < 1, |a e — j>£ + i/i| _1 < (a e —p 2 )~ l . In the contribution 
j4, resulting from the first term in the parenthesis, we have taken the supremum over all momenta 
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p e associated with the leaves in Q\. We will refer to this estimate as freezing these momenta and the 
corresponding legs e G Qi will be called dead-edges, the rest are called live-edges. In the contribution 
B we froze all momenta of the leaves Qi- 

In order to bound these integrals, we will successively integrate over all a variables and over all 
non-frozen momenta starting from the leaves, until we are left with an integral involving only the 
momenta of the roots. Finally, we integrate out the root momenta: at this point we will also make 
use of the decay factor Y\ e£ j^r)(Pe)~ 3 we gained from the test-function J^ k \ 

The large p e and large a e regimes are critical for the convergence of our integrals. The decay 
of the non- frozen leaf- momenta, |^» e | 2 , alone is not sufficient to render these integrals finite; for 
the intermediate edges even such decay is not available. The propagators provide extra decays, 
(a— p e ) _1 , but they may disappear (even with a possible logarithmic divergence) in the da e -integrals. 
On the other hand, the delta functions of course help since they reduce the effective number of 
integrals. Finally, the test-function provides a strong decay for the root variables (|9,45p . Due to the 
complexity of this structure, it requires a carefully designed successive integration scheme, combined 
with appropriate bounds, to show that these multiple integrals are actually convergent. The precise 
bound will then easily follow along the same lines. Unfortunately, the scheme is complicated by fact 
that in certain estimates (namely when Lemma 110.31 is applied) mild local point singularities arise. 
This is unavoidable even if one uses weights (p e ) 2 instead of pi in the Schwarz inequality in (|9.53j) . 
Therefore some care is needed to avoid accumulation of local divergences. We suggest the reader to 
neglect this issue at the first reading and concentrate only on the large momentum regime of the 
estimates in the proof of (|9.46[> . 

Next we illustrate the successive integration scheme for (I9.53P : we consider the term A, where 
the momenta of the edges in Q\ (recall Q\ is the set of outward leaves) are frozen. The analysis of 
the B is analogous. 

Since the delta functions always relate variables within the connected components of T, the 
integrations can be done independently in each connected component (tree) of T. The order of 
integration is prescribed by the tree structure: we start from the leaf-variables and proceed toward 
the root. 

The key step is what we call integrating out a vertex. It consists in integrating over the a variables 
and the momenta of the live son-edges of this vertex. The integral will be estimated in terms of the 
a-variable and the momentum variable of the father-edge and in terms of the frozen momenta of 
the dead edges from the set {£ £ Qi : v lies on the route from t to its root}. A vertex v will be 
integrated only when all vertices v' with v' >- v have already been integrated out. 

More precisely, we define an increasing sequence of subsets of the vertices V(T), V%(T) C V^(r) C 
. . . C V(r), where V^r) contains all vertices that have been integrated out after the first m inte- 
gration steps, in particular |Vm(r)| = m. In the (m + l)-th integration step we integrate out one 
of the maximal vertices in the set V(r) \ V m (T). The maximality is considered with respect to the 
ordering defined by the restriction of -< onto V(T) \ V m (T). After m = |V(r)| integration steps all 
vertices have been integrated out. 

The process of integrating out each last vertex v (a vertex whose father-edge is a root) in the 
2k — |i?i(r)| non-trivial connected components of T is a little bit different. In this case we integrate 
simultaneously over the a-variable of the son-edges and of the father-edge, and, like in the other 
vertices, we integrate over the momenta of the son-edges. Here we will estimate the integrals in 
terms of the momentum of the father-edge, and of the momenta of the dead leaves in the connected 
component we are considering. As a result, after integrating out all vertices of T, we will be left 
with an integral over the the momenta associated with the roots: the integrand will depend on the 
root-momenta, on the dead momenta and on the observable j( k \ 

Along the procedure we keep track of the available decay factors for each edges. Every edge 
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carries its own propagator, (a e — pi) -1 , and we will focus on the additional decay factors. We see 
from (|9.53p that the every momentum associated with a live leaf carries a decaying factor l/\p e \ 2 . 
We will show that when we integrate out a vertex, a similar polynomial decaying factor can be 
propagated to the father-edge. Because of the propagators (a e — pl)~ l , we will actually be able to 
gain more decay in the momenta of the father-edges, which will be available for the next integration, 
when the momentum of such a father-edge will be integrated out as a son-edge of the next vertex. 
The additional decay is typically in the form of point singularity with a power higher than 2 (i.e. 
of the form \p e — a\~ 2 ~ K , k > 0, with a possible shift a depending only on dead-momenta), but 
sometimes a so-called spherical decay in the form (a e + (3j + (p e — bj) 2 )~ l arises, where the shifts 
bj,Pj depend only on the dead-edge momenta. 

We will therefore distinguish different types of edges, according to the momentum-decay they 
carry. For e £ ^(r); we define the set of vertices 

V e = {v 6 V(Y) : e lies on the route from v to its root}, 

and the set of dead leaves 

D e = {£ £ Qi : e lies on the route from £ to its root}. (9.54) 

We choose A > and e > sufficiently small (the correct conditions will be specified later on). Then 
we have the following type of edges: 

i) d-edges: these are the dead edges, over whose momenta we do not integrate. Dead-edges are 
always leaves; note that one companion of each pair of leaves is dead, the other is live. 

ii) 2-edges: they carry a factor 

i 

(9.55) 



\Pe\ 2 ' 



These edges are exactly the live leaves. 

iii) (2 + A)-edges and (2 + 2A)-edges: they carry a sum of decaying factors 

"(e) j "(e) , 

El prrr, respectively pj— jt- (9.56) 
Ap e -aj\ 2 + x ' y ^ Pe - aj 2 + 2X 

where a,j are linear combinations of the momenta of the dead-edges lying in D e . Here the 
number of terms, v[e) is bounded by v(e) < C'^', for a universal constant C. 

iv) (2 + s + K)-edges with k = 0, A, 2A: they carry a sum of decaying factors of the form: 

"(e) 1 1 

E \ Pe - aj \ 2 +* (a e + (3 3 +(p e -b J ) 2 )^ ' (9 - 57) 

Here aj, bj, are linear combinations of the momenta of the dead-edges in D e , the numbers 
j3j are quadratic functions of the same momenta. The number of terms, v(e), is bounded by 
v(e) < C\ Vb \ for a universal constant C (the symbol s in (2 + s + k) refer to the "spherical 
decay" (a + + (p — a) 2 )~ 1+e ). If e is a root edge, then it can still be of the type (2 + s), 
(2 + s + A) or (2 + s + 2A), but in this case, in (|9.57j) we replace a e with p 2 . 



50 



(e r , a r , r) 



(e p , a p> P) 



(e q , a q , q) 

Figure 12: The vertex of transition 2) 

The summations in (I9,56h and (|9.57p reflect different cases that originate from the fact that the 
three son-edges of a vertex do not play fully symmetric roles. The precise number and the possible 
relations among a,- , bj , 0j of these terms play no role in the procedure since all our estimates will be 
uniform in these shift variables. The reader can therefore safely neglect the complicated structure 
of (19.56h and (|9.57p at the first reading. The only important issue is the type of singularity and the 
power k; these information are carried in the shorthand notation 2, 2 + k, 2 + s + k. 

We will show that every time three edges of one of these types meet as son-edges at a vertex, the 
father-edge will be again of one of these types after integrating out this vertex. We will prove the 
following transitions, to determine the type of the father-edge after integration, given the types of 
the son-edges. 

1) (d,d,2 + k) -> 2 + s + k, forK = 0,A,2A 

2) (d, 2 + «i, 2 + k 2 ) -»■ 2 + 2A, for k 1>2 = 0, A, 2A 

3) (2 + ki,2 + k 2 ,2 + k 3 ) -» 2 + 2A, for fti,2,3 = 0, A, 2A with k\ + k 2 + k 3 > 3A (9.58) 

4) (2 + 2A,2,2) 2 + A 

5) (2 + s + «i,2 + k 2 ,2 + Ka) ->• 2 + 2A for ki )2j3 = 0, A, 2A 

The short notation (A, B,C) — > D means that .A, i?, C type decays on the son-edges yields a -D-type 
decay on the father-edge after integrating out the vertex. The order of A,B,C is irrelevant. For 
example, transition 2) is a short-hand writing of the following estimate: 



v(e p ) u(e q ) 

su p H £ 

^ i=l 0=1 



5{a r = a v + a q — ot q i) da p da q da q '5(r = p + q — q')dpdq 
(ap — p 2 } \p — ai\ 2+K1 (a q — q 2 ) (a q i — (q 1 ) 2 } \q' — bj\ 2+K1 

v(e r ) 

r — cJ 2 + 2 -^ 



(9.59) 



L J \iy> r . 



k=l 



(see Fig. [12] for the notation, where we choose the dashed line with momentum q' to be the dead- 
edge for definiteness). Here c/% are linear combinations of aj, bj variables and of the dead edge 
momentum q' . 

The transitions (19.581) have to be combined with the following set of rules, to see that they indeed 
form a closed system along the successive vertex integration. The rules are as follows: 

a) At every vertex, there are at most two <i-edges. This is clear, because <i-edges are always 
outward pointing leaves, and it is impossible to have a vertex with three son-edges having the 
same orientation. 

b) At every vertex, there are at most two 2-edges: this follows analogously to a), because all 
2-edges are inward pointing leaves. 
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c) There is no vertex with son-edges (2 + A, 2, 2). Note that the (2 + A)-edge can only result as the 
father-edge of a vertex with son-edges (2,2,2 + 2A), and thus (since the 2-edges always point 
inward) it must have inward orientation. But all three son-edges of a vertex cannot have the 
same orientation. 

d) There is no vertex with son-edges (2 + A, 2 + A, 2). As we have seen in b) and c), all 2-edges 
and (2 + A)-edges have inward orientation and all three son-edges of a vertex cannot have the 
same orientation. 

Combining these rules with the transitions in (|9.58p , and observing that the spherical denominator 
(a + j3 + (p — aj) 2 )~ 1+e can be always estimated by one (i.e. +s can always be removed from 2 + s + k 
for free), we see that we have a closed system starting solely from d-edges and 2-edges. Along the 
integration, every edge in V will become one of the types described in i)-iv) above and each vertex 
integration corresponds to one of the steps 1)— 5). 

Before proving the transitions l)-5) rigorously, let us indicate their validity by a simple power 
counting argument. All integrals are locally convergent if A < 1/3, so it is sufficient to focus on their 
large momentum (short distance) behavior. The integration variables p e are momentum variables 
with dimension [length] . Because of the propagators (a e — p 2 )~ l , the a e variables have the same 
dimension as p 2 , i.e. [length]' 2 . Hence the three propagators associated with the three son-edges 
always have the dimension [length]®. 

In the transition 1), we have effectively two a integrations and no momentum integration. In 
fact, there are three a- variables and one live momentum associated with the three son-edges of the 
vertex under consideration; but, because of the ^-functions S(Y^, eev ±a e ) and S(Yleev ^Pe)i we om y 
have to perform two a integration. Each da- integration carries the dimension [length]' 2 . Since the 
momentum decay factor associated with the 2+K-edge has the dimension [length] 2+K , the result of the 
integral has the dimension [length]®~ 2 ' 2+<y2+K ^ ' = [length] i+K . The r.h.s. of 1), on the other hand, has 
the dimension [length] 4+K ~ 2e . Therefore e > guarantees that in the relevant short distance regime, 
the r.h.s is indeed bigger than the l.h.s., as the corresponding exponent on the r.h.s is smaller than 
the exponent on the left. We lose some decay in our estimates to compensate for logarithmic factors. 

In the transition 2), we have effectively two a-integrations and one momentum integration (dp e 
has the dimension [length]' 3 ): hence the left side of 2) has the dimensions [length] 3+Kl+K2 , while the 
r.h.s. has the dimension [length] 2+2X . For A small enough, the exponent on the right is smaller than 
the exponent on the left. 

In 3), 4), and 5) we have two a- and two momentum integrations. The dimension of the left side 
of 3) is [length] 2+Kl+K2+K ' i ; the right side has the dimension [length] 2+2X . This explains where the 
condition K\ + ki + K3 > 3A comes from. The transition 4) is similar. As for 5), the left side has the 
dimension [length] A+Kl+K2+K3 , and the right side has the dimension [length] 2+2X : again, for A small 
enough, the exponent on the left is larger than the exponent on the right. 

Next, we give a rigorous proof of the transitions (|9.58p . Choose a vertex v 6 V(r) and assume 
we already performed the integration over all v 1 >- v. We start by performing the integration over 
the a- variables associated with the son-edges of the vertex v. Let e r denote the father edge of the 
vertex v (see Fig. fT3l) . We distinguish two cases according to whether e r is the root or not. 

Suppose first that e r is not a root (that is v is not the last vertex left). If there is no spherical 
denominator among the son-edges, that is, if the son-edges are all of the type d, 2, 2 + A or 2 + 2A, 
then we estimate the a-integrations, for a fixed and small enough e > 0, by (with the notation of 
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(e r ,a r ,r) (e p , a p , p) 

(e tr , a,, q 1 ) 

(e q , a q . q) 



Figure 13: Integrating out a vertex 

Fig. [H: 



J da p da q da q '5(a r = a p + a q — a q >) — 



1 



p P 2 ) {otq - q 2 ) (ay - (q') 2 ) 
da p da q ^ _ ^ ^ _ ^ ^ _ a ^_ a ^ + ( 9 - 60 ) 

< I 

~ (a r - P 2 -q 2 + {q>) 2 y-z ' 

where we applied Lemma 1 1 . 1 1 twice . once in the a q and once in the a p -integration (after estimating 
(a) -1 < |a|~ 1+<5 , for a sufficiently small 5 > 0). On the other hand, if one of the son-edges (say the 
e p edge in Fig. [T3j) is of the type 2 + s, 2 + s + Aor2 + s + 2A, then we use the bound 

5(a r = a p + a q - a q >) 

da p da q da q i 



< 



(op - p 2 ) {a p + f3+(p- a) 2 y- s {a q - q 2 ) (ay - {q') 2 ) 

da p da q ^ _ ^ ^ + ^ + ^ _ a ^i-e _ q 2j ^ _ ap _ ^ + ( g /)2^ 

da p 

(a p ) (a p +p 2 + (3+{p- a) 2 ) l - e (a r - a p - p 2 - q 2 + {q') 2 ) 



1 fi f 1 1 ^ 1 

5, T7T- — 7T-, \oT~i 7 / da v - — r + 



(l3+p 2 + {p-a) 2 ) l - e J P \(a p ) (a p + + p 2 + (p - a) 2 ) J (a r - a p - p 2 - q 2 + (q') 2 ) 



< I ( I + 1 



~ (/? + (p - a) 2 ) 1 -" V (or " P 2 ~ q 2 + {q') 2 ) 1 - 6 (a r + (3+(p- a) 2 - q 2 + (q') 2 ) 1 ^ 

(9.61) 

where (3 and a, like (3 and a depend only on the frozen momenta associated to the dead leaves in 
D £r (see the definition (|9.54p ). 

Suppose now that e r is a root, that is there is no vertex v with v -< v. In this case we also 
integrate over the a variable associated with the father-edge e r . We use 



/ 



da r da p da q da q i5(a r = a p + a q — a q 



(ay - r 2 )(a p - p 2 ) (a q - q 2 ) {a q > - (q') 2 } 



da r da p da q ^ _ r 2 ){a p - p 2 ) (a q - q 2 ) (a r - a p - a q + (q') 2 ) ( 9 - 62 ) 
1 



(r 2 — p 2 — q 2 + (q') 2 ) 1 6 ' 
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r 



P 



Figure 14: A vertex with two dead edges 

if there is no spherical singularity in the son-edges, and 

5(a r = a p + a q — a q >) 



da r da p da q da q / 



(or - r 2 )(a p - p 2 ) (op + (3 + (p - a) 2 )^{a q - q 2 ) (cy - (q') 2 ) 



I ( I + 



^ +(p _a) 2 )i- £ V(r 2 -p 2 -g 2 + (g') 2 ) 1 - E (r 2 + /3 + (p-a) 2 -g 2 + ( (7 2 ) 1 - £ y ' 

(9.63) 

if there is a spherical denominator in the son-edges. 

Next we have to estimate the momenta integrations. Let us first consider the case 1) in (|9.58p . 
Since we have two dead edges and one momentum delta-function, effectively no integration needs to 
be done. With the notation of Fig. [HI where dashed lines indicate dead-edges, we obtain, using the 
result of MP . 

dg' S(r = p + q-q') 1 1 (g g4) 



\q' - a\ 2+K (a r - p 2 - q 2 + {q') 2 ) l - £ \r - b\ 2+K (a r + (3 + (r - c) 2 ) 1 ^ 

where k can assume the values 0, A, 2A, and /?, b, c depend on a and on the frozen momenta p, q. If 
e r is a root, then, according to (19. 621) . we replace a r by r 2 in (j9.64[) . This proves 1) (note that the 
two dead edges always have the same orientation: with the notation of Fig. [T4"l it is impossible, for 
example, that p and q' are both dead). 

The transition 2) is proven, using the result of (|9.60p (or (|9.62p . if the father-edge is a root) in 
Proposition 110.41 (where k±, K2 can assume the values 0, A and 2A), under the conditions A < 1/6 and 
e < 1/3. Note that the vertices of the type 2) involve two integrals, but because of the (5-function 
from the momentum conservation, one integral is a trivial substitution. 

The transitions 3) and 4) are proven in Proposition 110.61 using the result of the a-integration 
(|9.60p (or (|9.62p if the father-edge is a root) under the condition that A < 1/6 and e < A/2. Finally, 
the transition 5) is shown, using the result of (|9.6ip (or (|9.63p . if the father-edge is a root), in 
Proposition 110.71 (with A replaced by 2A), under the condition that k\ + K2 + k% < 2A and that 
A < 1/10 and e < A. If K\ + «2 + K 3 > 2A, then we can drop the denominator with the spherical 
singularity, and use the transition 4) to prove 5). Therefore, assuming that < A < 1/10 and 
0<e<A/2, we have proven the transitions l)-5) in (19.58p . 

After integrating over all vertices in the graph T we are left with an integral over the momenta 
of the roots in R(T)\Qi (a trivial root can be dead). Recall that we already integrated over the 
a-variable associated to the roots (we performed this integration together with the integration over 
the a- variables of the son-edges of the roots in (|9.62p and (|9.63p ). Estimating all the spherical 
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denominators left by one, we obtain from (|9.53|) that 
A < C n (J ' &v n+kP {...pl +k \^{p 



n+k)\ 2 



u(e) p \ (9.65) 

dp e 



iPeleGQx e eR(r)nQi ^ 7 eGi?(r)\Qi \j=i J 11 Jl ^ 7 



where, for every e, K e = 0, A or 2A, where a e j are linear combinations of the momenta associated 
to dead leaves in D e , and where the number of terms vie) is bounded by vie) < C'^ e L Since the 
integrals are uniformly bounded in the dead-momenta, and since X)eei?(r)\Qi l^el = n > ^ immediately 
follows that 

A<C n (j dp n+k pj . ..p 2 n+k |V(Pn +fe )| 2 ) • (9.66) 



In the same way we can prove that the term B in (19. 53ft satisfies the same bound. Hence we conclude 
that 

J dqfcdq^dr^fcdr^ J {k) (<i k ; q' fc ) K T)t (<ik, 0*; r n+k, r' n+k ) 7 {n+k) (r n+k , r' n+k ) 

< C"Tr (1 - Ax) ... (1 - A n+fe ) 7 (" +fc ) (9.67) 
for every t < 1. □ 

9.5 Proof of Theorem EH 

Despite the obvious analogy, the bound (|9.47p cannot be directly reduced to (|9.46p . The reason is 
that there are three propagators missing at the truncated vertex v that appears in the definition 
(|9.16p . Their missing decays need to be propagated through the whole integration procedure until 
the strong decay of the observable J will compensate for them. For this reason, edges carrying a 
momentum decay of the type i)-iv) introduced in (|9.55p - (|9.57p are not sufficient to prove (|9.47p . and 
we need to introduce additional types of decay. Of course this also requires to consider additional 
vertex integrations, slightly different from the transitions l)-5) in (|9.58p . 

Proof of Theorem \9.4\ As in the proof of Theorem 19.31 it is enough to prove (|9.47p for rank-one 
projectors, 7 (n+fc) (P«+fc; p' n+k ) = 4>(Pn+k)ip(p' n+k ), for a tp £ L^(M 3(n+fc) ); the general case follows 
then from the expansion 

7 (n+fc) (Pn +fc ;pU) = E A ^ J (P^)^(Pn+fc)> wi * h ^ e ^(M 3(n+fc) ), 11^-11 = 1, Vj 

j 

where the eigenvalues of ^( n + fe ) satisfy Xj > for all j and Ylj Xj < oo. 

In order to prove (|9.47p we start with the expression (|9.16p . Recall that if there is only one denom- 
inator containing a e , then the da e integral would not be absolutely convergent, so this integration 
has to be performed before taking the absolute value. This was the reason behind distinguishing the 
roots of the trivial components, Ri(T), in (19.9H . For the kernel Lr t, the same problem arises for the 
component containing v G M(T) (see ([9.16P ). if this component has only one vertex. This justifies 
the following definition. 

For a given T € T nk and v £ M(T), we define the set of edges E2(T,v) as follows: if there exists 
e G R(T) such that e G v, then E 2 (T,v) := E 2 (T)\{e}. Otherwise E 2 (T,v) := E 2 (T). In the former 
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case, starting from (|9.16p . we perform the integration over the «g using (|9.10p before we take the 
absolute value; in other words, we treat e as the trivial roots in Ri(T). This is necessary because, 
since the son-edges of v do not carry a propagator, there is only one denominator containing a g and 
the da e integral is not absolutely convergent. After performing the integration over all a e associated 
to e E" E2(T,v), we take the absolute value of the integrand, and we conclude that (|9.47|) is bounded, 
similarly to (|9.53j) . by (recall that we take here ^ n+k \p n +k; P' n+ k) = ^(Pn+fc)^(Pn+fc)) 



C n J dp n+k ( Pl ) 2 . . . (p n+k ) 2 \^(p n+k )\' 



sup yi [ n dpe n i — z^zii n 

>< n h n n *£±°o n *(E±^ (9 68) 
+ su P yj / n <fe n n 1 



{p e } eeQ2 _ gM(r) J eeB(r) \Q 2 ~ Fe -P e + tl egi?(r) \Pe 



x n ? n n n 

Recall that Sy denotes the set of son-edges of v and that Q\ and Q2 denote the set of outward and, 
respectively, inward leaves (j9.51|) . Note that the weights in the Schwarz inequality are somewhat 
different from the ones used in (|9.53|) : we take the decay factor {p)~ 2 instead of \p\~ 2 in the son-edges 
of the vertex v (the reason will be clear later on). Moreover, in contrast to (I9.53h . we keep track of 
the 1/t factors in the propagators to detect the short time behavior. To do so, we select rj e = t^ 1 
for all leaf-edges, e E L(T), and notice that for an arbitrary e E E(T) the value of r/ e is t" 1 times the 
number of edges in the subgraph of descendents of e. In particular 

Ve=t- 1 -(3n) 

ee-R(r) 

since the total number of non-root edges is 3n. Therefore every propagator in (19.161) carries at least 
the regularization i -1 (that is r/ e > t~ 1 , for all e E E(T)); on the other side, the exponential prefactor 
is at most e 3n = C n after taking absolute value. This justifies (|9.68p . 

To scale out the t variables, we rescale p e , for all e E E(T), and we rescale a e , for e E £b(T, v)\Sv 

as 

p e — * t~zp e for e E E(T) and a e — ► a e t~ x for e E ^(r, v)\S^. 

Then 

1 c -f^e 1 ti 1 e 1 e£v egi) eG-u eGu 
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a) 



b) 



Figure 15: Integrating out the vertex v 



From (|9"l?gj) we can bound (|9~27) by (recall t < 1) 

C n t^- 3 [ d Pn+k ( Pl ) 2 . . . (Pn+kf \^(Pn+k)\' 



X 



su p e / n d ^ n ^ d " e „2\ n 



(9.69) 



""F / / I I u Fe II / 2\11 I 

{Pe} eeQl eeM(r) •/ eeE(r) \ Ql ee E 2 (T,v)\S, ~ Pe) eeR(T) (t~~ 2 Pe) 3 

>< n ^ n tAi n n 

+ sup e / n n / d " e 2\ n ~ r — 

>< n A n 7A2 n *(E^) n *(E^ 

eeQi\5 s Fe eeQiDSs We/ »ev(r) eei) t)GV(r) eGi) 

= : A L + £ L . 

Note that here we rescale also the frozen momenta with t -1 / 2 : of course this is allowed, since we 
take the supremum over them. For e G Q 2 n S» in the term Al (and for e € Qi n 5^ in £?l), we 
used that (t _1 / 2 p e )^ 2 < ( Pe )~ 2 (because < t < 1). We show how to control ^4^, where the outward 
leaves are dead; the proof for Bl is then analogous. 

To estimate Al we proceed very similarly as in our analysis of the contribution A in (|9.53p . But 
here we first have to get rid of the vertex v, whose son-edges do not have propagators. We distinguish 
two cases. 

In the first case, we assume that v is a vertex involving two dead-edges (see Fig. [T5l part a)). 
Then the momentum integration at v (there is no a-integration here) gives 

~ iP + <1 ~ q ' )] = (r-p-q)* - \r- P -g\^ ^ 

for A > small enough. Edges carrying such a decay will be called (1 + 2A)-edges, similarly to the 
notation introduced in i)-iv) above (see (I9,55p - (l9.57p ). Note that in the last inequality, we dropped 
part of the decay for large momenta: this simplifies a little bit the classification of the possible types 
of edges in T (see the points v)-vii) below). Here we used the fact that in (|9.69p we chose the decay 
(p e )~ 2 instead of \p e \~ 2 for the son-edges of the vertex v to avoid an irrelevant complication in the 
short momentum regime. 

In the second case, if v only involves one dead-edge (there is always at least one dead edge 
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adjacent to v), then we obtain (see Fig. [151 part b)) 



/ 



-^-5(r -{p + q- q')) < . (9.71) 

\P/ \Q/ \r + q \ 



In this case, the father-edge of the ^-vertex will be called a type 1-edge. 

After dealing with the vertex v, we are faced with the problem of integrating out the other n — 1 
vertices of T. This is very similar to the problem we encountered when we proved the bound f|9.46[) : 
every non-trivial edge carries now a propagator (a e — pi) -1 exactly as in (|9.53p . Here an edge is 
called trivial if it is a trivial root or if it is a root adjacent to v. The difference is that now we start 
with more types of edges (in the analysis of Kr,t, we started solely with d- and 2-edges; here we also 
have a 1- or a 1 + 2A-edge): this implies that, integrating out the other n — 1 vertices, we will create 
new types of edges, which were not defined in i)-iv) (see fj9.55|) - (j9.5T[) in the proof of Theorem I9.3|) . 
Therefore, we have to supplement the definitions i)-iv) with the following new types of edges. 

v) 1-edges: they carry a decaying factor: 

1 



\Pe ~ a\ 

where a is a linear combination of the dead momenta in D e . As we will see, there can be only 
one 1-edge along the integration procedure; it is the father-edge of v, if it is not a root and if 
only one of the son-edges of v is a dead edge (case b) in Fig. [T5|) . In this case the two live 
son-edges have inward orientation, so does the father-edge of v, therefore the 1-edge is always 
inward. 

vi) (1 + A)-edges and (1 + 2A)-edges: they carry a sum of decaying factors 

"(e) ^ v{e) 

S3 |„ -o-ll+A ' res P ectivel y S3 i„ _„.|i+2A ( 9J2 ) 
3=1 J 3=1 J 

where a,- are linear combinations of the momenta of the dead-edges lying in D e . Here the 
number of terms, v{e) is bounded by v(e) < C 1 '^', for a universal constant C. 

vii) (1 + s + K)-edges with k = 0, A, 2A: they carry a sum of decaying factors of the form: 

E \ Pe - a ^ (ae + Pj + fa-btfy-e ' (9 - 73) 

Here a,j, bj are linear combinations of the momenta of the dead-edges in D e , the numbers 
f3j are quadratic functions of the same momenta. The number of terms, v(e), is bounded by 
v(e) < C'^', for a universal constant C (the symbol s in (1 + s + k) refer to the "spherical 
decay" (a + /3 + (p — a) 2 )~ 1+e ). If e is a root edge, then it can still be of the type (1 + s), 
(1 + s + A) or (l + s + 2A), but in this case, in (|9.73p we replace a e with p 2 e . 

Note that in our classification of the edges of T, we disregard the edges in S^: their only effect is to 
produce a (1 + 2A)- or a 1-edge. Notice also that the new types of edges appear only on the route 
from v to the root: all other edges are still either 2-, (2 + A)-, (2 + 2A)-, (2 + s)-, (2 + s + A)-, or 
(2 + s + 2A)-edges. In order to take care of the new type of edges we have to add to (I9.58P the 
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following new transitions: 

6) (d, d, 1 + k) -> 1 + s + k for k = 0, A, 2A 

7) (d, 1 + «i, 2 + « 2 ) -> 1 + 2A for ki >2 = 0, A, 2A 

8) (1 + «i,2 + «2,2 + k 3 ) -»• 1 + 2A for ki j2 ,3 = 0, A,2A with ki + k 2 + k 3 > 3A 

9) (1 + 2A,2,2) -» 1 + A (9.74) 

10) (1,2 + 2A,2) 1 + A 

11) (1 + s + ki, 2 + k 2 , 2 + k 3 ) -» 1 + 2A for « li2)3 = 0, A, 2A 

12) (1 + m, 2 + 8 + «2, 2 + « 3 ) -» 1 + 2A for « li2)3 = 0, A, 2A 

Moreover we have to add the following rules to the set a)-d) introduced above. 

e) At every integration step, among the son-edges there can only be one edge of the type v)-vii): 
this is clear because these edges can only be found on the route from v to the root. 

f) There is no vertex with son-edges of the type (1, 2, 2), (1, 2 + A, 2), or (1, 2 + A, 2 + A). This 
follows because there cannot be three son-edges with the same orientation, and because, from 
the previous observations (item b), c) and v)), any 1-edge, 2-edge and (2 + A)-edge has always 
inward orientation. 

g) There is no vertex with son-edges of the type (1 + A, 2, 2) or (1 + A, 2 + A, 2). This follows from 
the observation that the 1 + A edge, which can only result from a transition 8) or 9), has always 
inward orientation (because 1- and 2-edges have inward orientation). 

Taking into account the fact that the spherical denominator (a + j3 + (p — a) 2 )~ 1+e can be always 
estimated by one (i.e. +s can always be removed from the decay characterization of any edge), it is 
clear that the transitions 1)-12) together with the rules a)-g) define a closed system, so that every 
edge that arises in the successive vertex- integration of T (with the exception of the son-edges of v) 
is of either one of the types i)-vii) described above. 

Let us now prove the transitions 6)-12). The a-integration can be performed as in (|9.60j) . (|9.6ip 
(or in (|9.62p and (|9.63p . if we consider vertices adjacent to root-edges). As for the momenta integra- 
tion we proceed as follows. The transition 6) can be shown similarly to the transition 1) (see (]9,64j) ). 
The transition 7) follows, using the result of (|9.60p (or of (|9.62p . if the father-edge is a root) to bound 
the a-integration, from the second part of Proposition 110.41 under the condition that < A < 1/6 
and e < 1/3. The transitions 8), 9) and 10) follow, again with the help of (|9.60p or (|9.62p . from 
Proposition 110.61 under the assumption that < A < 1/6 and e < A/2. For ki + k 2 + k 3 < 2A, the 
transitions 11) and 12) follow, using the result of (|9.6ip (or (|9.63p . if the father-edge is a root), by 
Proposition 110.71 under the condition that < A < 1/10 and e < A (here we use this Proposition 
with 2A instead of A). If k\ + k 2 + k 3 > 2A, then we can drop the spherical denominator, and 11) 
and 12) follow from 8). Assuming that < A < 1/10 and < e < A/2, this completes the proof of 

Using the transitions (I9.58|) and ()9.74p we can iteratively integrate over all vertices in T, until we 
are left with an integral over the non-frozen root-momenta. From (|9.69p we obtain 

A L < C n f^ ( / dp n+k (pi) 2 . . . ( Pn+k ) 2 \4>(p n+k )\ 2 ) sup J ]'[ — r 




(9.75) 
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where we denote by e the unique root edge connected with the vertex v. As in (|9.65p . K e = 0, A or 
2A, the linear combinations of the momenta of the dead edges in D e , and v{e) < C'^L 

Note that the (1 + «)-type is inherited within the tree containing v, in particular the decay of the 
corresponding root-edge is weaker than that of all other root-edges. 
Rescaling the momenta, we observe that 

/dp i—K i dp 2 — k 

W^W^W - ct ^ ' and TV \ P -a\^(t^w - ct ^- (9 - 76) 

Since we assumed t < 1, since |i?(r)\Qi| > k, \M(T)\ < 2(n + k), and since K e < 2A, we obtain 

A L < (n + £;)C7 n i^- Afc TV (1 - A x ) ... (1 - A n+fc ) 7^+ fc ) . (9.77) 
Since we assumed n > 10 + fc/2, and A < 1/10, we find 

A L < C n i? Tr (1 - Ax) ... (1 - A n+fe ) 7 <»+*) . 
Since the same is true for Bl (see (19.69P ). this completes the proof of Theorem 19.41 □ 

9.6 Proof of Theorem ED 

Proof of Theorem \9.1[ Suppose that Ti :t = {7it}fc>i and T%,t = {72*}fc>i are * wo solutions in 
C([0, T\,TL) of the infinite hierarchy (|9.2p . such that, for j = 1,2, 7^ is non-negative, symmetric 

w.r.t. permutations and satisfies Ht^ \\n k < C k , for all k > 1 and t £ [0, T], and such that 7^ = 7^0 > 
for all A; > 1. We want to prove that Tit = ^2,t, for every t G [0,T]. To this end we will prove that, 
for every fixed k > 1, 7^ = 7^ for every t G [0, T] (as elements of Wfc). By a simple approximation 
argument it is then sufficient to prove that 

Tr J« (7$ - 7 g?) = (9.78) 

for all in a dense subset of the dual space oiTik- Since we assumed 7^ and 7^ to be symmetric 
w.r.t. permutations, it is enough to consider permutation symmetric observables j( k \ We will show 
(|9.78p for all permutation symmetric with kernel J^ k '(pkiP'k) C m momentum space) satisfying 



i^w;)i<cn--i 



For fixed k > 1 we can expand 7^ in a Duhamel-type expansion as in (|9.6p . With Theorem 19.21 we 
can identify each term in the expansion ([9.6p as the sum of contributions of Feynman graphs. We 
obtain, for any n, that 



n-l 



m=i re^ m , k rej„, t - 70 

for j = 1,2. Multiplying with the observable and taking the trace we obtain 



n-l „t 



(n+k)\ 



(9.80) 
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for j = 1,2. From Theorem 19.31 it follows that the terms in the sum over m are bounded in 
absolute value by C m \\^j^ +k ^ \\n m+k > i n particular they are finite. Since ^q™^ = 720 f° r ever y 
m > 1, when we take the difference between Tr «7^7i* and Tr J^^ti the free evolution terms 
{J^\U^^q) and all the terms in the sum over m disappear and it only remains to bound the 
contributions from the last term in (|9.80p . From Theorem 19.31 and since \ J- n ,k\ < C n+k , we obtain, 
under the assumption that t < 1 and n > 10 + fc/2, 

TV J« ( 7 g | < Cf» jT d S (t- S )i (||7& +fe) ||^ + Il7£ +fc) ||« fe+n ) < CM , (9.81) 

where we used that, by assumption, sup s6 j 0j T] Il7j" + ^ llw n +fc — C' n+fc for j = 1, 2. Hence, if we choose 
i < min(l, (1/2C) 4 ) we conclude that 

TrjW( 7 g- 7 ^) <2"». (9.82) 



Since n > 1 is arbitrary, this clearly proves (|9.78p for every t < min(l, (1/2C) 4 ). The proof can then 
be iterated to show that 7^ = 7^ for all t € [0, T]. □ 

10 Integrating Out a Vertex 

In this section we prove some estimates used in the proofs of Theorem 19. 31 and Theorem [93] to control 
the momentum integration in the transitions l)-5) in (|9.58p and 6)-12) in (|9.74p . 

10.1 Preliminary Estimates 

We begin by proving some useful lemmas: they contain the prototypes of integration we have to deal 
with when integrating out a vertex. 

Lemma 10.1. For every e, A, n with < e < A < 1 and < r/ < A — e there exists a constant C\ iEt11 
such that 

d " < T^Sh (10.1) 



for all ael. 



Proof. If \a\ < 1, then (a) ~ 1, (a — (3) ~ (/?) and (llO.ip is trivial. For \a\ > 1 we split the integral 
into two parts. For |/3| < |a|/2, we use 



1 < 1 



(a-p) 1 - 6 ~ (a-/3) 1 - A + J ?(a) A - £ - ? ? 
and we obtain 

d/3 C /■ d/3 



/" d l < C [ 



*W^{L»W + / d "((^ + <^ n (ia2> 

< c 

~ (a) A - e -^ 
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where we used the Schwarz inequality and the fact that, if \/3\ > 1/2, \/3\ ~ {(5). For \(3\ > \cx\/2, on 
the other hand, we use 

1 < I < 1 



\(3\ x ~ \a\ x -z~v\p\z+v ~ (a)^-'?!/?!^ 
and conclude similarly. □ 

Lemma 10.2. For every e, (5, 7 with < e < 1, 5 < (1/2) — e, 5 > —1/2, and < 7 < min(l — e; 1 + 

25; 1 — 25 — 2e), and /or every 77 > sufficiently small (depending on e, 6,j), there exists a constant 

C = Cs,ea,V w tth 



dp C 

\2-28 {a _ {p _ a)2) l-s ^ — ~ _ a2) i-l-s- ^ ■■ (J(U) 



-e-rj 



for all a £l and a G 



Proof. We consider first the case \a\ < 1. Then (a — a 2 ) ~ (a) and (a) ~ 1, so it is sufficient to prove 
the estimate (110. 3p when a 2 and (a) are removed from the r.h.s. of (I10.3|) . We now distinguish two 
cases, depending on the size of \a\. 

If \a\ < 10, then the integral I is comparable with 

dp f dp f ( 1 1 



1 £ / 119-M 7- ToZoJ £ / T-P^ox + / dp 



2-25 ( p _ a )2-2 £ |p|2-2« 1 J ^^)4-M-3te ^ _ a )4-2<5-2 £ 

which is uniformly bounded in a and a. Here we applied a Schwarz inequality, and we used that, by 
assumption, 5 > —1/2 and 5 + e < 1/2. Since in this case (a) ~ 1, this proves f)10.3|) . 

If \a\ > 10, then we split the dp-integration into three different regimes. In the regime p 2 < \a\/2 
we have (a — (p — a) 2 ) ~ (a) and, putting y = p 2 , 



I 

J v 



dp 1 f H/2 ^_ 1 ^ 1 



d?/ 771 7 ^5 — — i — ~r (10.4) 



W<M/2 N 2 " 2<5 (a- (p-a) 2 )i- ~ (a)i-* 7 |y|H ~ (a)^- s 

because 5 > —1/2. In the regime \a\/2 < p 2 < 2\a\ we have 



dp If dp 



:|/2< p 2<2|q| |p| 2 26 (a - (p- a) 2 ) 1 e ~ (a) 1 5 7| a |/2<p2< 2 | a | (a-(p-a) 2 ) 1 6 

<^—f dq dQ 

(a) 1 5 J\a\/2<(q+a)2<2\a\ (a-q 2 ) 1 6 (a) 1 S J\a\/3<q2<3\a\ {a - q 2 ) ] 

1 f 3H dy^y 1 /^H dy 1 



~ /|«|/3 (« " J/) 1 " 6 ~ (a)H 7|a|/3 l« " l/l 1 " 6 ~ (a)!— 5 

Finally, for p 2 > 2\a\, we have (a — (p — a) 2 ) ~ p 2 and hence 

dp < r dp < 1 



>>2|a| |p| 2 - 25 («-(p-a) 2 )l--7 p2 > 2H |p|4-2 5 -2 £ ~ (a) I 



-e-<5 



(10.6) 



Combining (|10.4|) . (|10.5|) . and (|10.6|) . we obtain (|10.3|) for arbitrary 7 > and rj > 0. 

Now we turn to the case \a\ > 1. By rotational symmetry, we can assume that a = (|a|,0, 0). 
After a change of variables and introducing g : = p\ + p 2 we find 

^J/^Jo \g + p 2 \ l - s (a-a 2 -g-p 2 + 2\a\ Pl y-c- 
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We define the new variables: 

u:=g + Pi, v := a — a —u + 2\a\p\. 
The map D 3 (u, v ) — ► (p\, g) £M x [0, oo) is one-to-one if we choose 

„ f / n ™9 /a — a 2 — u — u\ 2 ~i 

: 2H ) I' 

Computing the Jacobian of this transformation, we obtain 

1 f dudv 



I< 



\a\ J D lul 1-5 ^) 1 - 5 ' 
Using the definition of the domain D, we get, for < 7 < 1, 

I< — [ ^ < — / — (10 7) 



a\ Jn 1 1 1+2_,5 , > 1 

1 ^ 2+2 



-e 



Q — CL — U — V 



2 a 



|m| 2 + 2 \a — a 2 — u — v] 1 ^ ( v y-£ 



Applying the bound (llO.ip twice, to integrate first over v and then over u, and using the assumptions 
that < 7 < min(l — e; 1 + 25; 1 — 25 — 2e) and that 77 is small enough, we find 

I<JL[ ^ < Cs ^ e (10 8) 

~ M 7 JR (a - a 2 - n) 1 -^-! ~ (a^(a - a 2 ) i-*-«-e-iJ ' 

Here we used that \a\ > 1, to replace |a| by (a). □ 

Lemma 10.3. For any e, 5, 77 mf/i < e < 25 < 1, and < 77 < 25 — e, there exists a constant C§„e 
such that 

dp Cg,^ 

,|2+25 ( a _ p . a )l-e - ( a )2S-e- V \ 

for every a E M 3 , a£l. 



J= / ^ < (10 c 

1 I l_J9.4-9„>> /_. „ — /„AM-k-« I „ 11-2(5 ^ 



Proof. By rotational symmetry, we can assume a = (|a|,0,0). Introducing the variable Q = p\ 
we find 

J <r I dPl dQ < f dPl ' f ^ 



{a-pM) 1 - 6 Jo \PI + e\ 1+S ~ 7r (a- Pl \a\y-?\ Pl \ 2S \a\ l ~ 25 J R (a - y) l ^\y\ 25 
and we conclude by (jlO.ip . □ 

10.2 Momentum Integration 

Using the results of the last subsection, we provide here bounds for the integration over the momenta 
carried by the son-edges of a given vertex. The bounds in the next proposition are used to integrate 
out vertices with one dead edge, i.e. vertices of type 2) in (|9.58j) and type 7) in (|9.74p ). These 
vertices involve integration over two momenta; however, because of the momentum delta-function, 
effectively we only need to integrate over one momentum (see Fig. [T6|) . 
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Figure 16: The vertices of Prop. 110.41 



Proposition 10.4. Suppose < A < 1/6 and < e < 1/3. Let k\,K2 > 0, with k\ + K2 < 4A. Then 
there exists a constant C = C(A, e, k%, K2) such that 

f dp 1 ^ 1 

su p / tz — : P „ — : — ri^jz — ztttz — : — i^vw < - — r-^- (10.10) 



« 



\p — a\\ 2+Kl \r — p — q — a2\ 2+K2 (a — p 2 ± (r — p — q) 2 ) 1 6 \r — bj\ 2+2X 



and 



/dp 1 ^ 1 

|p-ai| 1+Kl |r-p-g-a 2 | 2+K2 (a-p 2 ± (r - p - q) 2 ) l - £ ~ |r-6j| 1+2A ( 10 - 1:L ) 

for any r,q,a\,a2 £ M 3 . -Here 6-,- are linear combinations of 0,1,0,2 fl^d 0/ i/ie frozen momentum q, 
and [i is a universal integer constant. 

Remark. The different signs in the propagator on the l.h.s. of (jlO.lOp are needed because of 
the two possible orientation of the dead edge. The positive sign in (jlO.lOp corresponds to the vertex 
on the left side of Fig. [16] since in this case the two live son-edges, carrying the propagators p 2 and 
[p + q — r) 2 , have opposite orientation. For the vertex on the right, first replace q with —q' in (jlO.lOp 
and then use the negative sign, corresponding to the paralell orientations of the live son-edges. The 
bound (llO.lip is used when one of the son-edges is a 1 + K-edge, for k = 0, A or 2A. The different signs 
in the propagator again take care of the possible orientations of the dead edge and of the 1 + K-edge. 

Proof of Proposition \10-4] We will make use of the following inequality to separate denominators. 

Lemma 10.5. For arbitrary a, (3 > and < 7 < min(a,/3), there exists a constant C a ^ n such 
that 



\a\ a \b-a\P~ |6|7 \\a\ a +P-~/ + \b - a\ a +^ ) ' (1 °' 12) 
Proof of Lemma \1(J.5[ Note that, since |6| < \b — a\ + \a\, we have |6| 7 < C 7 (|6 — a| 7 + |a| 7 ) and thus 

W <C 7 f^ + ^V (10,3) 



\a\~f\b - a|T ~ ' V|a| 7 \b - a| 7 

Hence 

1 1 1 < <h_(J_ + 



\a\ a \b-a\P \a\~t\b-a\~f \a\ a ~i\b - a\^ ~ |6| 7 \>| 7 \b - J \a\ a ~i\b - a\P-y 



<&.(, , ' + 1 



;io.i4) 



|6|t \\a\ a \b- a\^ \a\ a ^\b - a\P , 
and (|10.12p follows by a Schwarz inequality. □ 
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Returning to the proof of Proposition 110.41 we introduce a parameter 9 which can assume the 
values 1 and 2. In this way we can prove (jlO.lOp and (jlO.lip in parallel; we use 9 = 2 for the proof of 
(|10.1U|) . 9=1 for the proof of (j 1U. 1 1 1) . According to the sign in the denominator, we have to bound 
one of the two integrals 

dp 1 

p — a\\ e+K1 \r — p — q — a 2 \ 2+K2 (a — 2p ■ (r — q)) 1 ^ 6 

dp 1 ( 10 - 15 ) 



(I) := sup / 
a J 

(II) := sup / 

a J 



p - ai\ 0+K i\r - p - q - a 2 \ 2+K2 (a - 2(p - ^) 2 ) 1 - £ 



2 

Here we shifted the a variable by a p-independent number: this is of course allowed because we take 
the supremum over a. 

We consider first the integral (I). Using (|10.12p we obtain, for arbitrary — 1 < 7 < min(Ki, k 2 ), 

(I)< \ 

1 ' ~ \r-q-ai - a 2 \ 9 +^ 

x sup / dp ( - ———— + ' N 



\p - ai| 2+K i +K 2-T \r - q - p - a 2 \ 2+Kl+K2 ~' r J {a - 2p ■ (r - q)) 
1 [dp 1 



^ tz — : — : nor; sup 



\r - q - 01 - a 2 \ e +~< a J \p\ 2 +^+^~"/ (a - 2p ■ (r - g)) 1 ^ ' 
Assuming that 

£<ki + k 2 -7<1 (10.16) 
we can find rj > such that < rj < n± + k 2 — 7 — e and we can apply Lemma 110.31 to find 

1 



(I) 



,r - q - 01 - a 2 | e +T \r - g|i+7-«i-« 
< I 

~ \ r — g|6»+l+27-Ki-K2 \ r _ g _ fll _ a2 |6»+l+27-Ki-K2 

Choosing 27 = K\ + k 2 — 1 + 2A we can bound (I) by the r.h.s. of (jlO.lOp or of (jlO.lip (according to 
whether 9 = 2 or 9 = 1), with fj, = 2, 61 = q and 6 2 = g + a\ + a 2 : we only have to check that this 
choice of 7 is compatible with the condition —1 < 7 < min(Ki,K 2 ) and with ()1Q. 16j) . This follows 
from the assumptions that K\ + k 2 < 4A, A < 1/6 and < e < 1/3. 

Next we consider the term (II) in f 1 1 . 1 5 j) . Using (|10.12|) and changes of variables, we conclude 
that 

1 



(II) 



< 



a\ — a 2 \ 

x SU P / d P I -TiTZTT^ + 

a 



< 



| p _ ai |2+ K1+K2 | r _ p _ q _ a2 \2+K 1+ K 2 J ( a _ 2 ( p _ rr2)2)l-e 
! " (10.18) 



|r — g — ai — a 2 | 



dp / 1 1 
x sup / r-^-T^ , r— T ._„, ov1 - + 



|p| 2+K1+K2 \ (a- 2(p + 01 - ^S) 2 ) 1 ^ (a - 2(p + a 2 - E=2Wi- B 



2 y / \" "VP 1 2 

Applying Lemma 110.21 with 7 = 2A and 5 = — (k± + k 2 )/2, we obtain for < rj < 1/2 — A — e, 

(II) ~ 7- — — — -j* (7- — : 1 n„ m + n — : * n„ nn ) su p ' 



ai - a 2 \ e \\r - q - 2ai\ 2X \r - q - 2a 2 \ 2X J a / n \ 1+K i +K a 

\ ' (10.19) 



ai-a 2 | e+2A |r - q - 2a\\ e+2X \r - q - 2a 2 \ 0+2X ' 
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q =p+q-r 



Figure 17: The vertex of Prop. 110,61 

Here we used that A < 1/6, e < 1/3 and K\ + K2 < 4A. This completes the proof of (llO.lOp and 
flTI| . □ 

In the next proposition we show how to integrate out vertices where all the son-edges are live 
(2 + «)— or (1 + K;)-edges, with k = 0, A, or 2A. These will include the vertex integrations of type 
3) and 4) in (19. 58ft and type 8), 9) and 10) in (|9.74j) . After the a- integrations, these vertices involve 
integration over three momenta p, q, q' . Using the delta-function S(r — (p + q — q')), we are left with 
two effective integrations which need to be controlled (see Fig. fT7|) . 



Proposition 10.6. Suppose Ki,K2,K3 > with < k% + K2 + K3 < 1- Let < k < Ki + K2 + K 3, <ind 
e < (k\ + K2 + K3 — k)/2. Then there is a constant C = C(k\, ^2,^3, k, e) such that 

dp do 

sup ' 



\p — ai| 2+K1 \q — a2| 2+K2 \r — p — q — a3| 2+K3 



1 1 

-, 5 5 7 — - <CV 1 —rrr- (10.20) 

{a — p z — q z + [r — p — q) z ) 1 £ z — ^|r — 0j| z+K 



and 
sup 



dpdg 



\p — ai| 1+Kl \q — a2\ 2+K2 \r — p — q — a3| 2+K3 



~, 5 5 7 < CV -j — p- (10.21) 

(a -p 2 - q 2 + (r -p - g) 2 ) 1_£ ~[\ r ~bj\ 1+K 



and 



sup 



dpdg 



\p — ai| 2+Kl \q — a2\ 2+K2 \r — p — q — a3| 1+K3 



x 1 9 5 7 T9Tw < C 7" 1 T-nT7 ( 10 - 22 ) 

(a — p z — g z + (r — p — q) ) £ ^!^l r — "il 

/or all r, 01,02,03 £ R 3 . Here the bj are linear combinations of 01,02,03, and pi is a universal 
constant. 

Remark. The bound (|10.2ip is used when the edge carrying momentum p is a (1 + K)-edge, 
for k = 0, A, 2A. The bound (|10.22p . on the other hand, is used when the edge with momentum 
q' = r — p — g is a (1 + «)-edge (if the edge with momentum q is a (1 + K)-edge, then, after 
exchanging p <-> q, we can use (|10.2ip ). 
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Proof. In this proof we will assume that a\ = a 2 = 03 = 0: the generalization to 01,02,03 7^ can 
be obtained with similar shifts as we did in the proof of Proposition 110.41 Using 

1 /ll 1 



|p| K1 |g| K2 |r — p — q\ K3 \\p\ e \<l\ 9 \f — p — q\ l 

with 6 = K\ + K2 + Ks, and the symmetry of the integrand w.r.t. the exchange p «-> q, the l.h.s. of 
(fTO20j) is bounded by 



2 sup 



dpdq 1 



a 7 \p\ 2+e \q\ 2 \r — p — q\ 2 (a — p 2 — q 2 + (r — p — q) 2 ) 1 e 
dpdq 1 



+ sup 



a 



|2|g|2| r — p — q\ 2 +0 (a — p 1 — q 2 + (r — p — g) 2 ) 1 e 



~ SUP / |g| 2 |r-g| 2 j dp (\p\ 2 + 9 + \r - q - p\ 2 + 9 ) (a - p 2 - q 2 + (r - p - g) 2 ) 1 ^ (1 °' 23) 



^ sup 



dq f dp 1 



+ sup 

«2 



|o| 2 |r - g| 2 J \p\ 2+e (ai - 2r ■ q - 2p ■ (r - g)) 1_£ 
dq f dp 1 



\q\ 2 \r - q\ 2 J \p\ 2+e (a 2 - 2(q - §) 2 + 2p ■ (r - q)) 



l-e 



where we applied the inequality (jl0,12l) and then we shifted the variable q — > (r — p — q) to obtain 
the last term. Applying Lemma 110.31 we conclude that the l.h.s. of (|10.20|) can be estimated by 



sup 



dpdg 1 



2\l-e 



|p|2+Ki | |2+k 2 | r _ p — q,|2+K 3 _ p2 _ g2 _|_ (j, _ p _ 

f dq 1 

-TJw^Fm^ (10 - 24) 



+ sup 



dq 1 



, | g |2| r _ g |3-9( a2 _2(g-§)2)e-^' 
for any sufficiently small rj > 0. To bound the first term we use again (|10. 12j) and Lemma 110.31 
dg 1 



sup 

ai 



|g| 2 |r-g| 3 ~ e (ati - 2r-q) e -^ 

~]r^™?l dq (^ + | r _^) {ai -2r.q)°-e-V ^) 
< 1 1 1 

^ - — — S11D 7, < - 



because n < 6, < 2e < 9 — k and rj > is arbitrarily small. Using ()1Q. 12|) and Lemma 110.21 (with 
7 = k), the second term on the r.h.s. of (|10.24p can be controlled by: 

dg 1 



"af J \q\ 2 \r - (a 2 - 2(g - |) 2 )^ 

~ W S «?J dP + F^F*) («2-2(g-|) 2 )^ ( 10 - 26 ) 

1 1.1 



< — = — SU P Hi < 

"3 (a 3 ) 2 



1 !--'« SUp '/_ K -c-2ri ~ | r | 2 + K 
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because 9 > k and 2e < 9 — k, and n is arbitrarily small. This completes the proof of (|10.20p . 
Next we prove (|10.2ip . With 9 = k,\ + k 2 , we have, similarly to (|1Q.24|) . 
dpdq 1 



sup 



2\l-e 



\p\ \q\ 2+Kl \r — p — g| 2+K 2 ( a _ p2 _ g2 _|_ (j, _ p _ 

~ S « P / \p\\r-p\ 2 J dq (w T ° + \r-q-p\ 2+e ) (a - p 2 - q 2 + (r - p - q) 2 ) 1 ^ (10 ' 27) 

/" dp 1 [dp 1 

~ SUP / U 1^ _ r.13-0 Irv, _ 9 r, . r \e~e-ri + SU P 



\p\\r-p\ 3 - e <ai -2 P t)«- £ -i |p| |r - p| 3 ~ e <a 2 - 2 (p - §) 2 ) e -^ ' 

for any small 77 > 0. Using Lemma Il0.3[ the first term can be handled as follows 

dp 1 



sup 

01 



|p||r — p| 3 9 (ai—2p-r) t 



-e-ri 



~ (j^T + jT^r) <ai-2p 1 r)^ (10 ' 28) 

S rr su p — 5 — r ~ rr 

because 9 > 0, e < 9/2, and ry is arbitrarily small. As for the second term on the r.h.s. of (|10.27p . 
from Lemma |10.2| we find (using again that < e < 9/2 and choosing 7 = 0) 



sup 



dp 1 



a ?J \p\\r-p\^(a 2 -2(p-^)^-v 



< r- SUp / dp r-^g + rN2\0- e - n ( 10 - 29 ) 



3-0 | r -p|3-<7 (a 2 -2(p- §) 2 ) e - £ - r ' 



< 1 1 1 

£ rr su p — s — r ~ rr • 

Exactly the same proof also works for (|10.22p . □ 

Finally, in the next proposition, we show how to integrate out vertices with three alive son-edges, 
one of which carries a spherical denominator (that is one of the son-edges is a (2 + s + «)— or a 
(1 + s + K)~edge): these are vertices of the type 5) in (|9.58p and of the type 11) or 12) in (19.741) . 

Proposition 10.7. Suppose < A < 1/5 and < e < A/2. Let k\, k 2 , k 3 > 0, with K\ + k 2 + ft 3 < A. 
Then there exists a constant C = C(X,e,kx,K2,Ks) such that 



aJ,cJ \P ~ ci| 2+Kl |g - c 2 \ 2+K2 \r -p-q- c 3 | 2+K3 {/3+{p- c) 2 ) 1 ' 



dpdq 1 



1 



< 



c E urzhm (ia30) 



(a — p 2 ± q 2 ± (r — p — q) 2 ) 1 6 \r — 6j| 2+A 



and 

dpdq 



sup 

Q,/3,C 



|p - ci| 1+K1 |g - c 2 \ 2+K2 \r - p - q - c 3 \ 2+K ^ {(3 + (p - c) 2 ) 1 ^ 

— < g E |._[n + A (10-31) 



1 



(a — p 2 ± g 2 ± (r — p — g) 2 ) 1 £ ~^l r — ^j'l 
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and 



sup 



dpdq 



a,/3,cJ \P — ci| 2+Kl |<7 — C2| 1+K2 |r — p — q — C3| 2+K3 (/3+ {p — c) 2 ) 1_e 

1 M 1 

x 1 9,9,/ vKT-, < C Y^ 1 TnTT ( 10 - 32 ) 

(a-p 2 ±q 2 ± (r-p-q) 2 } 1 ' 6 j-^ \r - bj\ 1+x 

for every 01,02,0% £ M 3 . Here the bj are linear combinations of ci, 02,03 and fi is a universal integer 
constant. The bounds hold for all four possible choices of the two signs. 

Remarks. The bound (|10.3ip is used if one of the son-edges is a (1 + n)-edge and one of the 
other two son-edges is a (2 + s + K)-edge, with k = 0, A, or 2A. The bound (|1U.32|) . on the other hand, 
is used if one of the son-edges is a (1 + s + K)-edge. The different signs in the propagators are needed 
depending on the orientation of the edge carrying the spherical denominator (in our notation this is 
the edge with momentum p), and on which one of the two terms on the r.h.s. of (|9.6ip arising from 
the a-integration we are considering. Note that since one of the son edges always has an opposite 

orientation with respect to the other two, the sign combination ( ) in the propagators does not 

arise in our applications. Although the estimates remain true also in this case, we prove them only 
for the other combinations (++), (H — ) and ( ). 

Proof. We prove the proposition in the case c\ = C2 = C3 = 0: to generalize the proof for c\, C2, C3 7^ 
one can proceed as we did in Proposition 110.41 (since we assume c\ = C2 = C3 = 0, in our proof we 
will only need one term on the r.h.s. of (|1Q.30|) . (|10.3ip and (|10,32p . that is we can take [i = 1 
and b\ = 0; however, when ci, 02,03 do not vanish, one need fi > 1, as in Prop. 110. 4p . Without 
loss of generality we can also assume that e > (the l.h.s. of (|10.30p . (| 10.31 j) and (|10.32p is clearly 
increasing in e). We begin by proving (jl0.30j) and (|10.3ip : to this end we introduce the parameter 
9 which can assume the values 1,2 (to prove (|10.30p we use 9 = 2, to prove (|10.3ip we use 6 = 1). 
The possible combinations of the two signs lead to the three contributions 

dpdq 1 1 



(I) := sup 
(II) := sup 



|p|fl+Ki \ q \2+ K2 \ r -p- g|2+« 3 (p+(p- c) 2 ) 1 ^ (a - \{p + r) 2 + (q + E^) 2 ) 1 ^ 
dpdq 1 1 



a ,p,cJ \P\ 9+Kl \q\ 2+K2 \r-p- q\ 2+ ^ (13 + - c) 2 ) 1 ^ (a-(p-^) 2 -q-(p- r)) 1 " 
dpdq 1 1 



(III) . Slip j | p | e+Kl |^|2 +K2 | r _ p _ q \2+K, 3 (p + (jj_ c )2)l-e ( a _ p . r + q . (p _ r ))l- 

(10.33) 

The first formula corresponds to the sign choice (++), the second and third ones are the choices 

(H — ) and ( — h) (as remarked above, we do not consider, in this proof, the case ( ), because we 

do not need it in our applications). Recall that 9 = 2 is needed for the proof of (I10.30P and 9 = 1 
for the proof of (jl0.31 j) . To derive (|1Q.33|) . we shifted the variable a by some number independent 
of p and q: this is clearly allowed, because we take the supremum over a. 

We start by estimating the contribution (I). With k = K\ + K2 + K3, and using the assumption 
k<A, we find, from (jl0.12D . 

1 ^1/1 1 



< 



\2+K 2 \ r _ n _ n \2+K 3 ~ _ „|2-2A-f +e ^ |^|2+2A+f-e | r _ p _ ^|2+2A+f -e 



\p\ 6 + K i\q\ 2 + K *\r -p- q\ 2+K 3 ~ \p\ 9 +*\ r - p\' 



+ 



1 

~, 777, To n\~T~K~T~7 \ , o i o \ i « 77 



\ p \e\ r _ p |2-2A+f+ £ \ v | 9 |2+2A+f- £ \ r _ p _ q \2+2\+±- l 

(10.34) 



G9 



This implies, with a simple shift of the q variable, that 

dp 



(I) < sup f 



ajc J \p\ 9 ^\r - p| 2 ~ 2A -f +£ {P+(p- c) 2 ) 1 ^ 
dq 1 



+ sup 



|g| 2+2A+ f- £ (a - \{p + r) 2 + (q + 2?) 2 ) 1 " 6 
dp 1 



(10.35) 



;,Xc7 |p| e |r -p| 2 - 2A +f + £ (0 + (p - c) 2 ) 1 - 
dg 1 



x 



| ? |2+2A+f- e (a _ 1^ + r)2 + ((? + 2=21)2)1-5 



To bound the g-integrals, we apply Lemma [10. 21 with 7 = 1 — 2A — ^ — e, 25 = — 2A— f+e and r) = e; 
this is allowed, because 1 — 2A — | — e > (since e < A/2 and A < 1/5), because 1 — 2A — (k/2) — e < 
1 — 2A — (k/2) + e (since we assumed e > 0), and because 1 — 2A — (k/2) — e < 1 + 2A + (re/2) — 3e. 
We obtain (using that k < A < 1/5) 



(I) < sup / 



dp 1 1 



aJcJ \P\ 8+K \r - p|3~4A- K (p+(p- C )2)l- e <Q _ p . r) 2A+f-2e 

/• dp 1 1 

+ 3 P ci |p| 9 |r - p|3"4A {(3+(p- C )2)l- £ (a _ - . r) 2A+f -2e < 1(, - ;!( >> 



1 [dp 

< — rg SUp 



r\ e a/J M 3 " 4A (f3 + (p- c) 2 ) 1 - (a-p- r) 2A +! " 



-2<r 



where we used ()10. 12[) and a simple shift of the p-variable (and also of the variables a and c, over 
which we take the supremum). To estimate this term we use a Holder inequality: 



dp 1 1 

17cJ |p|3" 4A p + (p - c) 2 )i- (a - p • r) 2A +f - 2e 

dp 1 



l-2A-f+2e 



^ sup 

a,/3,c 



+ (/3+(p~c) 2 ) / (10-37) 



bl 

dp 




2A+f-2e 



A+Q/2)-2e / _ . \ 
Z+ 2A+(k/2)-2 £ \" ' / 



where we used that (l-e)/(l-2A-(«/2)+2e) > 1 (this follows from e < A/2) and that (/?-(p-c) 2 ) > 
1. The first integral is bounded uniformly in (3 and c: this follows from Lemma 110.21 (with 7 = 0), 
because -1 < (1 - 5A - (k/2) + 2e)/(l - 2A - (k/2) + 2e) < 1 (this follows from k < A < 1/5). As 
for the second integral, since < (A + (k/2) — 2e)/(2\ + (k/2) — 2e) < 1 (from e < A/2), we can 
apply Lemma ri0.3t 

f dp 1 1 < J_ nf) „^ 

Sy H 3 - 4A (/3 + (p-c) 2 ) 1 -(5-2p-r) 2A +f- 2£ ~ |r| A ' l±U ^ J 

The powers of |p| in (|10.37|) were chosen exactly in order to get this decay. With (|10.36p . we conclude 
that (I) < \r\- - x . 
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Next we consider the term (II) in f|10.33[) . Instead of (|10.34p we use here the similar bound 

(10.39) 



1 < I ( 1 + ! 



\p\9+Ki\q\2+K2\ r - p - q\2+K 3 ~ | p |0+K| r _p|2-2A-§ ^| g |2+2A+f | r _ p _ g |2+2A+f 



1/1 I 



| p |0| r _ p |2-2A+f ^| 9 |2+2A+f | r _ p _ ff |2+2A+* 



With a shift of the q variable, we obtain 

dp 1 f dq 



(II) < sup f 



ZJ \p\ e +-\r-p\ 2 ~ 2X -^ (P+(p-c)*)i-*J |g| 2 + 2A +f (a-ip-^-q-ip-r))^ 
dp 1 f dq 1 



+ S P J H^|r-p| 2 " 2A -f (/3 + (p-c) 2 )^7 |g| 2 + 2A +! ( a -p-r-g-(p-r))i 
dp 1 f dq 1 

+ sup ' 



a 



f dp 1 / 

J m<V-h 2 - 2A+ ^ {P + ( P - cyy-c J ' 



+ sup 



ftj M*|r-p| 2 - 2A +§ {(3 + {p-cYY-z J |g| 2 + 2A +f (a-(p-|) 2 -g-b-r))i- 
dp 1 f dq 1 



_p|2-2A+f ( / 5 + ( p _ c )2)l- £ y | g |2+2A+f ( Q _ p . r _ g . ( p _ r )) 1-e 

(10.40) 



Applying Lemma 110.31 (with r/ = e) to bound the ^-integral in the four terms, we find 

dp 1 1 



(II) < sup f 

a,0,c J 



ZJ |p| 9 +1r-p|3-4A-« {( 3 +{p - c) 2 ) l-e {a _ {p _r )2) 



2A+^-2e 



2' 



dp 1 1 

SU - P / h\ e + K \r -p|3-4A-« (p+(p- C )2)l-e ^ _ ~ . r )2A+f-2 £ 



a,/3,c 

sup 



dp 1 1 



(10.41) 



ZJ |pHr-p|3"4A ( / g + ( p _ c) 2 ) l-, (a _ (p _ i)2) 2A +f -2 £ 



sup 



2/ 

dp 1 1 



ZJ \ P \ e \r-p\ 3 - 4X </3+(p-c) 2 )i-<a-p. r ) 2 A+f-2 £ - 

Using that k < A < 1/5, it follows from (jlO. 12[) and from a shift of the p variable (and of the a and 
c variable as well), that 



I ' I a.p.c J 



dp 1 1 

L'' 

1 (10.42) 



ZJ N 3 ~ 4A (/3+(p-c) 2 ) 1 - £ («-(p-|) 2 ) 2A+ f' 



•2e 



1 f dp 1 1 

+ Tifl SU P 



rf Q /c7 H 3 " 4A (/5+(p-c) 2 ) 1 - £ (a-p. r ) 2A +f- 2£ " 

The second term is identical to the r.h.s. of (|1U.36|) . and can be estimated in the same way. As for 
the first term on the r.h.s. of the last equation, we apply a Holder inequality (we use here the same 
exponents as in (|10.37p but here we divide the powers of |p| in a different way): 



sup 



dp 



ZJ M 3 ~ 4A (P + (P ~ C) 2 ) 1 ^ (a - (p - I)2)2A+f -2 £ 



1-2A-- £+2e 



^ sup 



dp 1 \ [ f dp 



2A+f-2e 



a,Z\J l^+ T^f+Ii (1+(P-C) 2 ) ) \J |p| 2 (a-(p-|) 2 ), 

(10.43) 
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Since -1 < (1 — 4A)/(1 — 2A- («/2) + 2e) < 1 (as follows from A < 1/5), the first integral is bounded 
uniformly in (3 and c by Lemma 110.21 To estimate the second integral we use again Lemma I10.2[ 
with 7 = A/(2A + (k/2) - 2e) (this is allowed because A/(2A + (k/2) - 2e) < 1 for s < A/2). We 
conclude that 

/"dp 1 1 < 1 1 

"Z/ M 3 ~ 4A (P~(P- c) 2 V- £ (a - (p - I)2)2A+f-2 £ ~ (^A ^ |^A • (°- ) 

From (fTIM it follows that (II) < |r|- e - A . 

Next we consider the term (III) in (110.33 j) . With (jl0.39j> . we find 

dp 1 f dq 1 



(III) < sup / 

aS,cJ 



p,J \p\e+-\ r -p\ 2 ~ 2X -^ {P + {p ~ c) 2 ) l ~ e J \q\ 2 + 2 ^ (a-p-r + q-(p-r)) 

f dp 1 f dq 1 

+ ST i P / | m |(9+«U_„|2-2A-$ /fl 4- Cti - ^2\l-s / |„|2+2A+# 



a/J |p|^|r-p| 2 - 2A -f (P + ip-c)*)^ J |g| 2 + 2A +! (a-(p-|) 2 + g-(p-r))i- 
dp 1 f dq 1 



+ sup 



+ sup 



Ui \p\ e \r-p\ 2 - 2X+ % (P+ip-c) 2 ) 1 -* J |g| 2 + 2A +f (a-p-r + q-(p-r)y-z 
dp 1 f dq 1 



a/ci \p\e\ r -p\ 2 - 2 ^ ((3+(p-c) 2 y-c J \ q \ 2 + 2 ^ (a-(p-^ 2 + q-(p-r)y-e- 

(10.45) 

These terms can be estimated as we did with the four terms on the r.h.s. of (|10.40p (the different 
sign in front of the factor q ■ (p — r) plays no role in our bounds). 

As for (I10.32p . similarly to (I10.33P we find the three contributions 



, n f dpdq 1 1 

1 >' a /J H 2 + K M9l 1+K2 k-p-<z| 2 + K M/5 + (p-c) 2 ) 1 -M«-l(p + 2 + (9 + £ ? : ) 2 ) 1 - £ 

(IF) := sup f 
(IIF) := sup j 

a,0,c J 



4\y t- ' j t \q -r — 
dp dq 1 1 

\p\2+Ki\q\l+K2\r - p - q\ 2 + K 3 {[3 + (p - c) 2 ) 1 ^ (a - (p - §) 2 - g • (p - r)) 1 " 
dp dg 1 1 



Ac-' |p| 2+Kl |g| 1+K2 |r -p - g| 2 +«3 (/J + (p - c ) 2 ) l - £ (a - p • r + q ■ (p - r)) 1 ^ ' 

(10.46) 

The analysis of (F)-(IIF) is then very similar to the one of (I)— (III) : the only difference is that instead 
of using (|10.34p and (|10.39|) . we employ 

1 < I ( : + 



|2+ Kl | g |l+ K2 | r _ p _ q \2+K 3 ~ | p |2+ K | r _ p |l-2A-f+e ^ | g |2+2A+f -e | r _ p _ ^|2+2A+f - £ 

1 / 1 1 



|p|2| r _p|l-2A+f+e \ v | g |2+2A+f- £ | r _ p _ g |2+2A+f- £ ^ 

(10.47) 



to bound the term (F), and 

f*a ~ |„|2+k|^ _ „|l-2A-§ \ |„|2+2A+f T |„_^_^|2+2A+^ ) 

(10.48) 



| p |2 +Kl | 9 |2+ K2 | r _ p _ 9 |2+ K3 ~ |p|2+«| r _ p |l-2A-f \^|2+2A+f ^ _ p _ g |2+2A+f 



1/1 1 



|p|2| r _p|l-2A+f V| g | 2 + 2A +f | r _p _ g|2+2A+f 

to bound (IF) and (IIF) (we use (|10.47p when the propagator is quadratic in q, ([10.48P when it is 
linear). □ 
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A Some Technical Bounds 



In this Appendix we collect some simple information that are used throughout the paper. 

Lemma A.l. Let on be the scattering length of j^Vn with Vjy given in Sl.l\) and let &o := J* V/v = 
f V. We assume additionally that V is radially symmetric. Then for any < /5 < 1 

lim Na N = . (A.l) 

Proof. Let R be the radius of the support of V, i.e. V(x) = for \x\ > R. An upper bound for ajv 
can be obtained by the inequality (see, for example |23j) 

n Vn = it ■ (A - 2) 

To derive a lower bound for a at, we recall that the scattering length can be computed through the 
integral 

8na N = j ^V N f , (A.3) 
where f(x) is the radial symmetric solution of the zero energy equation 

with f(x) — > 1 as \x\ — > 00. It is easy to show that 

„ , J 1-S for \x\ >a N 

fix) > < 1 1 , • (A.4) 

J v 1 ~ \ for \x\ < a N y ' 

This can be proven as follows. We write f{x) = g{\x\)/\x\. Then g satisfies the one-dimensional 
equation 

Simple arguments show that g(r) > for all r > 0. Moreover, for r > RN~@, we have g(r) = r — aj\f 
(it is easy to see that this definition of the scattering length a agrees with (IA.3j) ). For r < RN~@ we 
have 

'■RN-P 



g{RN-P) - g(r) = / dsg'(s) 

J r 

= J ds Ig'(RN-P)- j drg"(r)\ (A.5) 



r-RN-f 3 rRN-P JU3/3-1 

(RN-P-r)- J ds J dr — — V{N P T)g(r) 



because g'(RN~^) = 1. Since g(RN~@) = RN' 13 - a N , we obtain 

i-RN-f 3 i-RN-f 3 

r-g(r)=a N - ds dr N^^V^N^g^r) . (A.6) 
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From the non-negativity of the potential and from g > 0, it follows that r — g{r) < a at for all 
r < RN~@ ', and hence f(x) = <?(|x|)/|x| > 1 — ajs?/\x\. Eq. (|A.4p now follows because /(x) > for 
all x e R 3 . 

Inserting (|A.4p into (|A.3P , and using the bound (|A.2j) , we conclude that 



ttNci n > I dxN w V{N p x) ( 1 



&o /" dx 



> / dx N 3/3 V(N x) \= / ^ N 3p V{N x 



x\>b /8nN J \x\ 

>b - f dx N^V(N^x) - ^N?- 1 [ p- V(x) 

J\x\<b /87rN 87r J \ x \ 



(A.7) 



>b -*-( b -^) 3 \\v\\~-^ moo [ £ 

3 V 87r / 87r J\x\<R Fl 



l\x\<R \ x \ 

This, together with (|A.2p . implies that 

b - dN 3/3 - 3 - C 2 N^ 1 < 8ttNcl n < b , (A.8) 
for two ^-independent constants C\,C%. Hence, for < (3 < 1, we obtain (jA.ip . □ 

In the next lemma we prove that solutions of the nonlinear Schrodinger equation which are in 
the space H 1 ^?) at time t = 0, have H 1 norm uniformly bounded in time. 

Lemma A. 2. Suppose ip G i/ 1 (R 3 ), and let ipt be the solution of the nonlinear Schrodinger equation 

id t (ft = -Atpt + b \(p t \ 2 (pt (A.9) 

with bo > 0. Then 

(<p t , (1 - A)<p t ) = [dx (|^ t (x)| 2 + |V^(x)| 2 ) < C (A.10) 



for all t € R. Hence, with 7 t ^' ) (x/ c ; x' fc ) = J^jLi ^(^jVt^j); we ^ aue 



Tr|5i...5 fc7 f ) 5 fc ...5i| <C fc (A.ll) 



XYie constant C only depends on bo and on the H l -norm of ip. 
Proof. The L 2 -norm of <p>% is conserved in time. Also the energy 

E(<p) = J dx|V^(x)| 2 + | J dx|^(x)| 4 

is conserved. By the Sobolev inequality, we have 

Jdx\v<p(x)\ 2 < E(ip) < cywjj, , 

for a constant C only depending on 6 - Hence 

J |Vv9 t (x)| 2 < E(<pt) = E(<p) < C\\<pf H1 . 



□ 
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The following lemma is useful to bound the pair interaction Vn(x) in terms of the kinetic energy. 
Lemma A.3. Let V e L l {^?) n L 3 / 2 (M 3 ). Put V N {x) = N^V {N ^ x) . Then 

V N (xi - x 2 ) < (const.)\\V\\ L i (1 - Ai)(l - A 2 ) 

and 

V N (X! - x 2 ) < (consL)N p \\V\\ Li /2(l - Ai) 
hold with universal constants. 

Proof. For a proof of these results, see Lemma 5.2 in [9]. □ 

Finally, we give a proof of Lemma [8.2l that is a slight modification of the proof of Proposition 8.1 
in [§]. 

Proof of Lemma\KM By the positivity of j( k+1 \ it is enough to prove (|8.8)1 for the special case 
-y( fc+1 )(x/ c+1 , x' fc+1 ) = /(xfc +1 )/(x^, +1 ). We can then bound the l.h.s. of f|8.8j) by the sum 

J dx fc+ idx^. +1 J (fe) (x fc ;x' fc ) (<5 ai (4+i ~ a*+i) ~ <K4+i - a?fc+i)) ^Oj - x k+1 ) f (x k+1 )J (x k+1 ) 

+ J dx fc+ idx' fe+1 j( k Xx k ;-x' k )6(x , k+1 - x k+1 ) (5 a2 (xj - x k+1 ) - 5(xj - x k+1 )) /(x fe+1 )7(x' fc+1 ) . 

(A.12) 

The first term can be bounded by 

dx fc+1 dx' fc+1 J (fc) (x fc ;x' fc )(5 Ql (4 +1 - x k+l ) - 5(x' k+1 - x k+1 ))5 a2 ( Xj - x k+l ) f (x k+1 )J (x' k+1 ) 

;x' fc )| 5 a2 (xj - xk+i) |/(x fc+ i)| 
x / d 4+i d «i ( x k+l ~ 4+1 ) [/(**» «k+i) - /(**> 4+1 )] • 

(A.13) 

We use the estimate 8 ai {x) < 4 • 1b (a?) where 5 := {x : \x\ < a{\. A standard Poincare-type 
inequality (see, e.g. Lemma 7.16 in [13]) yields that 



/ dx' k+1 5 ai (x k+1 - x' k+l ) [f(x' k ,x k+1 ) - /(x' fc ,4+i)] <C 
J J\y\<oti 



|V fc+ i/(x^,a; fe+1 + y)\ 

~T2 d y 



\y\<ai \V\ 

(A.14) 

for any xi and x k+ \. Inserting this inequality on the r.h.s. of (|A. 13|) and applying a Schwarz 
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inequality we get 

J dx k+1 dx' k+1 J (fe) (x fe ;x' fc ) (S ai (x' k+1 - ajfc+l) ~ S ( x 'k+1 ~ x k+l)) Sa 2 {xj - x k+1 ) f {x k+1 )J (x k+l ) 

<C f dx fc+1 dx^ ai) 5 a2 ( Xj - x k+1 )\ J (fc) (x fc ;x^)| 

x (|/(x fe ,x fc+ i)| 2 + |Vfe+i/(xj b ,x fc+ i +y)| 2 ) 
<Cai^sup y dx' fc | J (fc) (x fc ;x^.)|^ ^ dx fe+ i 5 a2 (xj - x k+1 )\f(x k , x k+ i)\ 2 

+ c| ; SU P y dxfc^^ - x k+ i)\ J (fc) (x fc ;x^,)|^ 

J dx' fe drr fc+1 dy — ^-p — |Vjfc + i/(x' fe , + y)| 2 . 



X 



(A.15) 



In the first term we apply Lemma IA.3I and in the second term we shift the x k+ \ variable, and then 
we compute the y-integral. Moreover we use that 

sup/dxUJ (fc) (x fc ;x' fc )|<C fc |||j( fc )|||, (A.16) 

and 

sup / dx fc )\J^(x k ;x' k )\<C k \\J^l (A.17) 

for a universal constant C (recall the definition of the norm I J W I j from (HZD). Thus, we find 

J dx fc+ idx' fc+1 J {k \x k ;x' k ) (5 ai (x k+l - x k+1 ) - 5(4+1 ~ &k+l)) S a2 (xj - x k+1 ) f (x k+1 )J (x' k+1 ) 
< C k ai \l Jl% TV (1 — A,)(l - A fe+1 ) 7 ( fc+1 ) . 
In order to control the second term on the r.h.s. of (|A.12p . we use that 

J dx fe+ idx' fe+1 j( k \x k ;x k )8(x' k+1 - x k+1 ) (5 a2 (xj - x k+1 ) - 5{x.j - x k+1 )) f(x k+1 )J(x' k+1 ) 

= J dx k dx' k dx k+1 {5 a2 (xj -Xfc+i) - 5(xj - x k+1 )) J (/c) (x fc ; x' k ) f(x k , x k+1 )f(x k , x k+1 ) 

= - j dxi . . . dxj-i&Xj+i . . . dx k+l dx' k J(x' k , x k+1 ) 

x ^J (fc) (x fc ;x' fc )/(x fc ,x fc+1 ) - J dxjd a2 (xj - x k+1 ) J (fc) (x fe ; x' k )f(x k , x k+1 ) 



(A.18) 



where we introduced the notation x^ = (x\, ... , Xj-i,x k +i,Xj+x, . . . , x k ) (that is x k is the same as 
x k , but with Xj replaced by x k +\). Using again the generalization of the Poincare inequality that 
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led to (|A.14|) . we obtain 

J dx fc +idx' fc+1 J^(x k ;x' k )6(x' k+1 - x k+1 ) (S a2 (xj - x k+1 ) - 5{xj - x k+1 )) f (x k+ x) f (x' k+1 ) 

<C f d Xk dx k dx k+1 1{lx { ~ * k+l] f 2 a2) V 3 [J^(x k ;x' k )f(x k ,x k+1 )] 1/(4,^+1)1 
J \ x j "Pfc+il 

I a a i a - x k+ i\ < a 2 ) 
< C I dx fc dx fc dx fc+ i J - 12 



x (jVjJ^x^x'JI \f{x k ,x k+1 )\ |/(x' fe ,x fc +i)| + |j( fc )(x fc ;x^)| | V i /(x fc , x k+1 )\ |/(x' fc , x k+1 )\ 

iV^^Cxfc;^)! ( K |/(x fc ,x fc+1 ))| 2 + K- 1 |/(x , fc ,^ + i)| 2 ) 



l(\xj - x k +i\ < a 2 ) 



< C / dx fc dx' fc dx fc+1 ^ ^ 



+ 



C [ dx.dx^+x 1 ^ f 2 a2) |J (fc) (x fc ;xl.)| («|V i /(x fcj x fc+1 )| 2 + «- 1 |/(4^fc+i 



(A.19) 

In the terms proportional to k we drop the restriction l(|xj — x k+ \\ < a 2 ) and we apply the Hardy's 
inequality to the x^+i-integration. In the terms containing K — 1 , on the other hand, we perform the 
Xj integration (after estimating | J^\ and \\7jJ^\ by their supremum). We get 

J dx fc+ idx' fc+1 j( k) (x k ;x' k )5(x k+1 - x k+1 ) (5 a2 (xj - x k+1 ) - 5(xj - x k+1 )) f (x k+1 )J (x' k+1 ) 

< Ck (sup J dx' k (\J^(x k ;x' k )\ + iV.J^Xfcjx',)!)) J dx k+1 |(1 - A fc+1 ) 1 / 2 (l - A,) 1/2 /(x fe+1 )| 2 

+ CK~ 1 a 2 l sup / dxi . . .dxj_idx J+ i . . .dx fc (| J (fc) (x fc ;x' fc )| + |VjJ (fc) (x fc ;x' fc )|) 

x j dx^dx fc+ i |/(x' fc ,x fc+ i)| 2 

< C k {K + aa/c- 1 )!! J% Tr (1 - A i )(l - A fc+1 ) 7 ( fc+1 ). 



1 /2 

Choosing k = a 2 , we find (jS 



□ 
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